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Preface 


As  our  nation  moves  toward  the  goal  of  a  manned  presence  in  space, 
we  have  come  to  the  realization  that  we  must  make  every  effort  possible  to  reduce 
our  dependancies  on  Earth  resources  if  we  are  to  maintain  this  presence  and  our 
economic  well  being.  The  high  cost  of  transporting  fuels  and  resources  into  space 
and  returning  waste  products  to  Earth  may  make  our  stay  in  space  a  relatively 
short  one.  We  must  begin  to  nurture  and  research  a  space  self-sufficiency 
concept  which  will  allow  us  to  exploit  the  resources  of  space  without  doing  the 
same  to  the  Earth. 

I  have  often  been  intrigued  by  the  man  who  powered  his  home  with  a 
wind  generator  and  the  farmer  who  ran  his  farm  on  the  waste  of  his  pigs  and 
cattle.  Why  not  this  type  of  self-generated  power  in  space?  is  it  economically 
feasible  to  develop  a  system  that  would  turn  our  waste  products  in  space  into  fuel 
energy?  This  research  project  explores  this  question  with  regard  to  the  fuel  energy 
created  being  applied  to  maintain  a  space  station’s  orbit. 

The  energy  generating  concept  proposed  in  this  study  is  anerobic 
digestion.  This  process  has  four  benefits  for  space  application;  1)  it  can  stabilize 
human  waste  products,  2)  it  can  reduce  solid  wastes,  3)  it  can  provide  a  fairly 
clear  effluent  for  water  recovery,  and  4}  it  can  provide  a  fuel  in  the  form  of  a  gas. 
The  method  envisioned  for  employing  this  energy  producing  process  was  to  first 
determine  a  scenario  so  that  input  load  to  the  digester  system  would  be  known  and 
also  the  size  of  the  spacecraft,  which  would  be  "powered"  by  this  energy  source, 
could  be  defined.  The  size  and  shape  of  the  vehicle  is  necessary  in  determining 


the  atmospheric  drag  which  must  be  opposed  to  maintain  the  orbit. 

To  properly  analyze  this  concept  it  was  necessary  to  research  three 
distinctly  unrelated  areas;  1)  simulation  analysis  of  biochemistry,  2)  thermochemical 
analysis  and,  3)  cost  analysis.  In  performing  the  cost  analysis  it  was  necessary  to 
determine  an  alternative  system  for  comparison.  The  alternative  selected  was  the 
transport  of  conventional  propellants  from  the  Earth.  This  alternative  does  not 
consider  a  replacement  of  the  anaerobic  digester  with  some  other  system  which 
must  be  present  to  stabilize  the  waste  products  of  the  space  station.  Nor  does  it 
consider  the  additional  benefits  of  the  anaerobic  digester  listed  in  the  preceding 
paragraphs.  In  this  respect  the  analysis  can  only  be  regarded  as  a  partial  one; 
however,  its  utility  lies  in  its  use  as  a  comparison  tool  which  can  be  applied 
toward  analysis  of  any  other  waste  stabilization  system  that  may  be  selected. 
The  benefits  may  be  taken  as  decision  variables  which  would  be  considered  in  any 
decision  of  a  waste  treatment/fuel  source  decision. 

A  good  deal  of  the  foundation  for  understanding  the  biochemistry 
involved  in  anaerobic  digestion,  and  the  analytical  descriptions  of  the  methane 
generation  process,  came  from  Price  and  Cheremisinoff’s  book,  Biogas  Product  ton 
8  Utilisation.  Their  equations  provided  the  primary  relationships  upon  which  the 
simulation  model  was  built.  During  the  study  it  was  realized  that  the  methane 
produced  by  the  anaerobic  digester  would  have  to  react  with  an  oxidizing  agent  if 
there  was  any  hope  of  the  system  making  a  significant  contribution  to  the  orbit 
maintenance  requirement.  Air  Force  Institute  of  Technology  (AFIT)  professor  Dr. 
William  Elrod’s  assistance  in  the  thermochemical  analysis  of  this  thesis  is  greatly 
appreciated.  He  dedicated  many  hours  to  enhancing  my  understanding  of  a  subject 
area  I  knew  very  little  about.  I  thank  him  also  for  the  guidance  he  provided  as  a 
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reader  in  keeping  this  report  accurate  and  on  track.  The  computer  program  used  in 
the  thermochemical  anlaysis  was  a  modification  of  a  program  by  Capt.  Robert 
Demmick.  His  generosity,  in  letting  me  use  his  program,  was  essential  in 
completion  of  this  project.  Dr.  Lynn  Wolaver,  Dean  for  Research  and  Professional 
Development  at  AFIT,  was  also  beneficial  as  a  reader  in  keeping  the  biochemistry 
aspects  of  this  thesis  accurate  and  concise.  I  truly  appreciate  the  giving  of  his 
time  from  his  busy  schedule  for  the  benefit  of  this  project.  My  advisor,  Dr 
Joseph  Cain,  Associate  Profesor  of  Economics  at  AFIT,  deserves  a  great  deal  of 
recognition  for  his  guidance  throughout  this  project,  especially  in  the  cost  analysis 
portion;  thank-you  for  your  support.  The  experts  also  deserve  recognition  for 
their  time  and  support  in  providing  me  with  the  inputs  I  needed  to  configure  the 
digester  system  and  cost  the  component  parts;  I  thank  David  Hill  and  Richard 
Westerfield,  waste  treatment  engineers  at  the  Dayton  Municipal  Waste  Treatment 
Facility.  I  also  thank  Gary  Lubin  at  The  Henry  P.  Thompson  Company  for  his 
assistance  in  the  estimation  of  many  of  the  digester  system  components  A 
special  thanks  to  Kathy  Cook,  superintendent  of  the  waste  treatment  facility  for 
the  city  of  Fairborn,  for  providing  me  with  many  of  the  EPA  and  Water  Pollution 
Control  Federation  documents  used  as  references  in  this  research  work. 

There  are  two  other  people  who,  though  they  did  not  contribute  to  this 
thesis  directly,  were  fundamental  in  its  development  and  completion.  The 
appreciation  I  have  for  the  support  I  have  received  from  my  wife,  Edna,  and  my 
daughter,  Janae,  cannot  be  expressed  in  words.  I  believe  their  devotion  to  this 
effort  has  been  greater  than  mine.  Without  them  this  thesis  would  not  have  been 
possible.  Thank  you  both. 
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v  Abstract 

This  project  explores  the  economic  feasibility  of  creating  fuel  energy  in 
space  from  human  waste  with  application  toward  space  station  orbit  maintenance. 
The  energy  generating  concept  proposed  in  this  study  is  anaerobic  digestion.  This 
process  has  four  benefits  for  space  application;  1)  it  can  stabilize  human  waste 
products,  2)  it  can  reduce  solid  wastes,  3)  it  can  provide  a  fairly  clear  effluent  for 
water  recovery,  and  4)  it  can  provide  a  fuel  in  the  form  of  a  gas. 

The  analysis  is  dependent  upon  a  predetermined  scenario  defining  the 
input  load  to  the  digester  system  and  the  size  of  the  spacecraft.  The  size,  shape, 
and  altitude  of  the  vehicle  determine  the  atmospheric  drag  which  must  be  opposed 
to  maintain  the  orbit.  The  basic  elements  of  the  study  involve  1)  simulation 
analysis  of  biochemistry,  2)  thermochemical  analysis  and,  3)  cost  analysis  using 
the  Monte  Carlo  method.  An  alternative  system  to  which  the  digester  is  compared 
is  transport  of  conventional  propellants  from  Earth.  This  alternative  does  not 
consider  a  replacement  of  the  anaerobic  digester  with  some  other  system  to 
stabilize  the  waste  products  of  the  space  station,  or  the  additional  benefits  of  the 
anaerobic  digester  listed  above.  In  this  respect  the  analysis  can  only  be  regarded 
as  a  partial  one;  however,  its  utility  lies  in  its  use  as  a  comparison  tool  which  can 
be  applied  toward  analysis  of  any  other  waste  stabilization  system  that  may  be 
selected.  ~The  results  of  this  study  show  a  statistically  significant  advantage  of 
the  digester  system  over  transported  conventional  propellants  due  to  the  high  cost 
of  space  transportation^  Recommendations  are  to  investigate  other  altitude 
scenarios  and  to  compare  the  digester  to  other  waste  stabilization  methods. 
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AN  ECONOMIC  FEASIBILITY  STUDY  ON  THE 
SPACE-BASED  PRODUCTION  OF  METHANE  GAS  FROM 
HUMAN  WASTE  THROUGH  ANAEROBIC  DIGESTION 
FOR  USE  AS  AN  ORBIT  MAINTENANCE  PROPELLANT 


I.  INTRODUCTION 


Background 

Current  procedures  in  disposal  management  of  human  waste  in  the 
space  environment  involve  thermal  drying  or  space  vacuum  drying  of  fecal  solids 
and  return  of  these  solids  to  Earth.  Both  of  these  methods  require  energy 
expenditure.  Heat  is  required  in  the  thermal  drying  process,  and  energy  is 
expended  when  returning  the  residual  waste  mass  to  Earth  through  a  controlled 
reentry.  As  man  extends  his  time  in  space  the  need  to  establish  requirements  for 
handling  metabolic  wastes  within  his  space  habitat  over  long  duration  missions  will 
be  the  focus  of  much  research  and  development  efforts.  These  efforts  will 
concentrate  on  the  most  efficient  means  of  reducing  and  utilizing  the  waste.  The 
high  costs  of  space  transportation  places  additional  emphasis  on  recycling  the 
wastes  of  a  space  habitat.  By  reclaiming  energy  from  the  waste  materials  through 
an  anaerobic  digestion  process,  reduction  of  these  costs  may  be  possible.  (31:6; 
32:82;  33:18;  34:53;  35:922) 

In  future  long  duration  manned  space  missions,  life  support  systems  can 
be  expected  to  supply  by-product  gases,  mainly  hydrogen,  methane  and  carbon 
dioxide.  All  are  candidate  gaseous  propellants  which  can  be  used  for  spac^  raft 
attitude  control.  These  gaseous  outputs  can  be  provided  by  various 


bioregenerative  systems  that  use  microorganisms  to  stabilize  waste  materials 


before  they  build  up  to  toxic  levels.  With  such  systems,  human  and  other  waste 
products  in  the  spacecraft  are  potential  sources  of  future  propellants.  Utility  in 
the  production  of  energy  biogas,  through  anaerobic  decomposition  of  human  waste 
solids,  may  provide  a  viable  energy  source  on  long  duration  space  missions  where 
quantities  of  wastes  must  be  managed  to  provide  the  least  amount  of  energy 
expenditure  for  the  most  amount  of  energy  return.  (45:228,469) 

Though  a  limited  amount  of  scientific  study  has  been  done  on  biogas 
production  from  solid  waste  materials  on  Earth,  research  proposals  have  only  been 
suggested  in  this  area  of  solid  waste  management  applicable  to  space. 

Problem  Statement 

The  technological  area  of  methane  generation  from  human  waste  has 
received  very  little  research  attention  for  spacecraft  application.  Sufficient 
engineering  data  relating  to  the  biological  treatment  of  concentrated  wastes  for 
precise  design  calculations  of  components  to  be  used  in  space  are  not  available. 
Most  of  the  technological  developments  to  date  are  the  products  of  research  and 
applications  relating  to  municipal  waste  processing.  (30:24;  36:266) 

The  use  of  methane  as  a  gaseous  propellant  for  attitude  control  on 
space  craft  is  already  a  well  proven  concept  and  hence  is  not  presented  as  a  main 
topic  in  this  thesis  research.  However,  analysis  of  methane  quality,  quantity, 
collection  method,  and  the  biogas  process  is  presented  in  determination  of  concept 
validity.  A  specific  scenario  based  on  space  platform  mission  provides  the 
boundaries  for  modeling  these  factors.  The  scenario  is  based  on  previous  studies 
performed  by  NASA  which  were  concerned  specifically  with  designing  space 
colonies.  In  determining  the  economic  feasibility  of  the  methane  generation 
concept,  a  comparison  of  acquisition  and  operating  costs  between  propellant  supply 
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via  shuttle  and  space  based  methane  generation  is  necessary.  (30;  34;  45:228) 


Purpose  of  Study 

As  manned  mission  requirements  begin  to  extend  up  to  one  year  and 
beyond,  research  and  development  efforts  will  seek  life  support  systems  and  waste 
subsystems  which  will  reduce  weight  and  volume  requirements  of  the  total  manned 
vehicle  complex.  At  the  same  time,  increasing  fuel  and  transportation  costs  will 
require  more  conservative  and  possibly  less  sophisticated  fuel  sources  if  we  are  to 
remain  guardians  of  the  high  frontier.  The  anaerobic  digester  provides  one  such 
waste  subsystem  which  can  reduce  metabolic  waste  volume  and  provide  a  useful 
by-product  gas  which  may  satisfy  some,  if  not  all,  of  the  orbit  maintenance  fuel 
requirements. 

The  focus  of  this  initiative  is  to  determine  the  cost  advantage,  if  any, 
of  generating  methane  gas  from  human  waste  in  space  (for  use  as  an  orbit 
maintenance  propellant)  compared  to  providing  propellants  from  Earth  via  the 
shuttle.  The  term  orbit  maintenance  is  considered  here  as  it  applies  to  overcoming 
the  Earth’s  atmospheric  drag  on  a  permanently  manned  space  station,  since  this  is 
the  force  that  finally  removes  the  satellite’s  energy  and  causes  it  to  spiral  inwards 
toward  Earth.  This  problem  will  be  addressed  in  four  phases:  1)  development  of 
an  applicable  scenario  to  define  propellant  requirements,  2)  sizing  of  the  methane 
generation  system  in  an  attempt  to  meet  these  requirements,  3)  estimation  of  the 
acquisition  and  operational  cost  of  the  system,  and  4)  comparison  of  the  system 
costs  with  the  logistics  cost  of  transporting  propellants.  (44:  229) 


The  broad  scope  of  this  proposal  requires  a  method  of  using  established 
theory,  simulation  modeling  and  cost  analysis,  together  with  sound  engineering 
principles.  Existing  data  is  used  to  the  greatest  extent  possible  and  all  methane 
generation  modeling  and  volume  and  weight  estimates  incorporate  mathematical 
relationships  found  in  current  literature  applicable  to  Earth  based  systems.  Cost 
estimations  for  this  study  are  based  on  a  comparison  method  or  use  current 
available  costs  on  items  having  equivalent  or  existing  counterparts. 

The  station  will  be  assumed  to  have  pseudogravity,  due  to  the  human 
physiological  need  for  gravity  on  long  duration  space  missions.  Pseudogravity  is 
artificial  gravity  created  by  rotation  of  the  space  station  about  a  central  axis. 
Use  of  analysis  based  on  standard  lg  Earth  gravity  is  appropriate  in  specifying  a 
system  that  will  operate  in  a  lg  pseudogravity  environment.  (34:21) 

The  space  station  will  not  be  totally  self  sufficient  and  will  be 
dependent  on  Earth’s  resources  for  such  things  as  food.  The  location  in  space  of 
the  station,  therefore,  must  be  within  the  range  of  the  shuttle  and  will  be  in  a  low 
earth  orbit  [LEO]. 

Size  and  weight  of  the  space  station  is  determined  by  the  living  and 
working  space  required  for  the  crew  members.  In  estimating  crew  size,  area 
requirements,  and  structural  mass,  data  from  two  programs  in  engineering  design  of 
habitats  for  sustaining  life  in  space  on  a  large  scale  will  be  used.  The  first  design 
program  was  a  10  week  study  held  at  Stanford  University  and  the  Ames  Research 
Center  of  the  National  Aeronautics  and  Space  Administration  in  1975.  The  results 
were  published  in  a  NASA  pamphlet  titled,  Space  Settlements,  A  Destgn  Study. 
The  second  program  was  a  similar  study  that  took  place  in  1977.  The  technical 


papers  resulting  from  that  program  were  published  in  a  NASA  pamphlet  titled, 
Space  Resources  and  Space  Settlemenis.  This  study  develops  its  scenario  from 
NASA’s  work  involved  with  the  initial  LEO  operations  which  provide  a  staging 
base  for  assembly  of  lunar  orbit  payloads,  space  manufacturing,  and  laboratory 
experiments  in  simulated  gravity.  A  modular  habitat  design  is  assumed  to  allow 
for  future  expansion  and  will  accommodate  150-250  people.  (30:66-71;  34:47). 

Determination  of  the  methane  generator  performance  and  size  is 
dependent  on  the  rate  of  methane  expenditure  for  attitude  control  and  can  be 
interpreted  as  the  generation  rate  requirement.  Factors  affecting  methane 
generation  are  primarily  temperature,  moisture,  input  feed  rate,  residence  time  of 
waste  material  in  the  generator,  degree  of  mixing,  and  loading  rate.  Using  design 
criteria  and  formulas  presented  in  Price  and  Cheremisinoff’s  book.  Biogas 
Production  and  Utilisation ,  a  continuous  simulation  model  is  developed  in  SLAM 
[simulation  language  for  alternative  modeling]  with  the  above  factors  representing 
the  input  parameters.  Input  feed  rate  is  dependent  on  the  number  of  crew 
members  selected,  and  human  waste  input  rates  are  based  on  data  obtained  from 
studies  done  by  NASA.  (35:916,  41:97) 

Since  there  is  no  night  or  day  on  a  space  station  in  a  LEO,  it  is 
assumed  that  individuals  aboard  the  station  will  work  on  shift  schedules  extended 
over  24  hour  periods.  Therefore  waste  feed  input  to  the  digester  will  vary  on  a  six 
hour  cycle  following  what  would  be  the  natural  waste  disposal  flow  of  morning, 
noon,  evening  and  midnight  activities  of  the  crews. 

It  has  been  shown  that  the  main  methane-producing  bacteria  are  the 
mesophiles.  These  organisms  are  most  active  in  the  30  -  40  C  range.  Above  this 
temperature  no  significant  increase  in  gas  production  occurs  until  approximately 


55C.  At  this  temperature  the  thermophilic  bacteria  are  favored.  Although  there  is 
a  potential  for  increased  gas  production,  the  thermophlicic  temperature  range  is 
rarely  used  because  it  is  impractical  to  provide  the  heat  necessary  to  favor  the 
bacteria.  This  study,  therefore,  will  consider  operating  temperatures  of  the  digester 
in  the  mesophilic  range.  (41:121) 

One  of  the  main  objectives  of  all  space  based  systems  is  the  reduction 
of  weight  by  integrating  as  many  systems  as  possible.  It  is  assumed  that  the 
anaerobic  digester  proposed  in  this  study  will  be  a  subsystem  of  the  life  support 
system,  and  will  be  expected  to  provide  a  fairly  clear  effluent  for  water  recovery. 
The  digester  shall  therefore  be  a  two  stage  system  in  which  the  first  tank  shall  be 
used  for  digestion  and  be  the  primary  methane  generation  source.  The  second  tank 
will  be  for  storage  and  concentration  of  digested  sludge  and  will  provide  a 
relatively  clear  supernatant  to  a  water  recycling  system.  (41:93) 

Though  analysis  of  the  water  that  may  be  recovered  from  the  various 
cleaning,  food  preparation,  and  bodily  waste  activities  has  not  been  reported  in  the 
literature,  it  may  be  assumed  that  these  waste  waters  will  resemble  domestic 
sewage  with  respect  to  the  biological  oxygen  demand  [BOD],  Calculations  for 
methane  generation  rates  involving  BOD  will  use  values  from  current  literature  on 
domestic  sewage.  (36:275) 

Cost  estimates  for  this  initiative  include  the  acquisition  cost  of  the 
methane  generation  system,  costs  associated  with  placement  and  operating  the 
system  in  space,  and  conventional  propellant  acquisition  and  shipping  costs  for 
comparison.  Since  a  system  of  this  type  has  not  been  designed  for,  or  utilized  in 
space  before,  costs  are  derived  by  comparison  with  similar  systems  on  Earth.  The 
life  cycle  cost  of  the  system  is  calculated  for  a  useful  lifetime  of  approximately  30 


years.  Conventional  propellant  costs  and  shipping  costs  are  based  on  current  year 
dollars. 


Methodology 

The  first  step  in  handling  a  problem  of  this  size  is  to  partition  the 
problem  into  its  separate  research  components  and  to  identify  the  relationships 
among  those  components.  These  relationships  are  shown  in  Figure  1.  The 
following  discussion  will  cover  each  of  the  four  large  blocks  of  research. 


Figure  1.  Problem  Research  Procedure 


Space  Station  Scenario.  Grew  size,  altitude,  and  size  and  mass  of  the 
space  station  are  all  based  on  previous  studies  done  by  NASA.  Once  location 
and  mass  are  determined,  the  force  necessary  to  overcome  atmospheric  drag  is 
calculated  from  aerodynamic  principles.  This  force  is  then  interpreted  in  terms  of 
fuel  requirements  necessary  to  produce  the  change.  The  use  of  methane  in  a 
cold-gas  system  and  as  a  fuel  in  a  bipropellant  system  is  considered  in 
determining  which  method  best  meets  orbit  maintenance  requirements.  A  cold-gas 
orbit  maintenance  system  is  one  which  uses  inert  gas  jets  incorporating  fast-acting 
valves  and  receiving  its  propellant  supply  from  pressurized  cold  gas  supply  tanks. 
A  bipropellant  system  requires  two  separate  propellants,  an  oxidizer  and  a  fuel.  In 
this  case,  methane  provides  the  fuel.  Various  oxidizing  agents,  such  as  pure 
oxygen,  are  considered.  (30;  34;  45:204,  227) 

Digester  Performance.  The  digester  is  modeled  using  design  criteria  and 
formulas  presented  in  Price  and  Cheremisinoff’s  book,  Btogas  Production  and 
Utilization.  The  model  is  formulated  using  FORTRAN  statements  in  a  SLAM 
continuous  simulation  model.  Input  is  feed  rate  based  on  the  space  station 
population.  The  model  allows  for  variation  of  the  feed  rate,  as  well  as  operating 
temperature  ranges.  Outputs  from  the  model  are  the  volume  rate  of  digester  gas 
produced  and  the  volume  requirement  of  the  digester  design.  Analysis  of  the 
output  data  will  determine  the  proper  combination  of  design  parameters  that  meets 
the  methane  rate  requirement,  and  will  also  determine  the  required  digester  volume. 
If  the  rate  of  methane  production  does  not  meet  the  requirement,  all  further 
costing  comparisons  will  be  based  on  the  portion  of  the  requirement  that  can  be 
met.  (41) 

The  modeled  digester  volume  is  used  to  estimate  the  weight  of  the 


system  components  using  comparable  structural  weights  of  hardware  and  materials 
appropriate  for  use  in  anaerobic  digesters.  Many  of  the  component  parts 
necessary  to  operate  the  digester  have  been  sized  by  comparison  to  a  small  scale 
portable  digester  built  by  Bio-Gas  of  Colorado  and  through  consultation  with 
waste  water  treatment  engineers.  (2.102) 

Cost  of  Digester.  One  of  the  more  difficult  phases  of  this  initiative  is 
estimating  the  cost  of  the  methane  generation  system.  Literature  searches  have 
not  produced  data  on  which  cost  estimates  may  be  based.  Therefore,  the  method 
used  to  base  costs  is  comparison  against  similar  systems  on  Earth.  Factors  used 
in  estimating  cost  are  the  size  of  the  digester  required,  selection  of  construction 
materials,  fabrication,  weight,  temperature,  and  mixing  requirements.  Local  experts, 
such  as  waste  treatment  plant  contractors  and  engineers  were  contacted  to  obtain 
this  information.  The  cost  of  placing  the  unit  into  space  is  based  on  shuttle 
transportation  rates  for  weight  and  volume  of  the  selected  design. 

Operating  costs  are  derived  from  estimated  labor,  maintenance  and 
energy  requirements  based  on  a  comparison  of  the  designed  facility  with  a 
municipal  waste  treatment  plant  having  anaerobic  digesters.  The  energy 
requirement  is  based  on  a  standard  percentage  of  total  facility  cost.  (47) 

In  costing  a  system  of  this  type  an  important  consideration  is  that 
stabilization  of  human  waste  must  be  performed  by  one  means  or  another  to 
prevent  the  spread  of  disease  by  reducing  pathogenic  organism  content.  Therefore, 
the  system  will  have  advantages  associated  with  its  function  as  part  of  a  closed 
cycle  life  support  system.  This  closed  cycle  system  will  be  required  to  produce  a 
fairly  clear  effluent  to  be  used  for  drinking  water  recovery.  There  are  many  waste 
treatment  techniques  and  facility  designs  capable  of  performing  these  functions. 


Therefore,  all  costs  incurred  in  acquiring,  operating  and  placing  an  anaerobic 
digester  in  space  should  be  weighed  against  the  alternative  systems  costs  in 
providing  waste  stabilization  and  against  current  water  recovery  methods.  The 
depth  of  research  and  computational  resource  requirements  involved  in  a  study  of 
this  complexity  is  beyond  the  scope  of  this  project.  This  study  then  represents 
only  a  partial  analysis  which  ran  be  used  as  a  comparison  against  alternative 
systems.  There  are  many  factors  which  may  influence  the  selection  of  a  particular 
system.  One  such  consideration  is  the  advantage  associated  with  handling  a 
reduced  amount  of  waste  solids  from  the  digester.  These  waste  products  could  be 
vacuum  dried  using  current  waste  handling  methods  and  returned  to  Earth  or 
possibly  utilized  in  space  in  such  projects  as  space  produced  fertilizer  for  crop 
growth  on  future  fully  self-contained  space  stations.  (39:923) 

Cost  of  Shuttled  Propellants.  The  final  phase  of  this  initiative  will 
compare  the  methane  generator  cost  to  the  cost  of  supplying  the  required  orbit 
maintenance  propellant  by  the  shuttle. 

Since  nitrogen  is  the  most  common  coid  gas  propellant  used  in  orbit 
maintenance  systems,  it  will  be  used  as  the  comparison  standard  in  costing  the 
delivered  propellant  price  if  the  cold-gas  method  of  orbit  maintenance  is  selected 
as  best  in  meeting  orbit  maintenance  requirements.  Determination  on  how  much 
nitrogen  is  required  will  be  made  by  using  a  mass  density  comparison  of  nitrogen 
to  that  of  digester  quality  methane  based  on  the  volume  of  methane  requirement 
established  in  the  scenario  phase.  (45:228) 

To  determine  if  the  methane  generated  would  best  be  used  in  a 
bipropellant  system  to  meet  orbit  maintenance  requirements,  comparison  against 
conventional  shuttle  orbit  maintenance  bipropellants  {Monomethythydrazine  [MMH] 


and  Nitrogentetroxide  [N204]}  is  made.  Determination  of  the  conventional 

bipropellant  quantity  necessary  to  meet  orbit  maintenance  requirements  is  made  by 
comparing  the  exhaust  velocity  or  the  specific  impulse  of  the  conventional 

bipropellants  with  that  of  the  methane  generated,  when  used  in  a  bipropellant  orbit 
maintenance  mode.  Specific  impulse  is  a  rocket  performance  parameter  which 
relates  thrust  to  the  propellant  mass  flow  rate  [See  Appendix  A  for  a  more 
complete  definition],  (1:13-8;  45:29,  297) 

The  cost  of  the  required  conventional  propellants  is  based  on  current 
market  price.  Shuttle  transport  cost  is  based  on  current  payload  cost  rates  where 
both  propellant  and  container  weight  will  be  considered.  (34:145) 

The  value  of  methane  generation  in  space  for  orbit  maintenance  is 

determined  by  tabulating  the  various  costs  incurred  in  each  system,  calculating  a 

total  life  cycle  cost,  and  comparing  these  costs.  In  using  the  generated  methane 
as  a  bipropellant,  the  additional  cost  of  providing  an  oxidizer  propellant  is 
considered.  Uncertainties  resulting  from  estimations  or  assumptions  are  reflected 
by  stating  a  range  in  which  the  prices  are  valid  and  an  uncertainty  associated  with 
the  price  range.  All  costing  data  acquired  during  this  study  are  standardized  to 
current  dollars  prior  to  making  a  valid  cost  comparison. 

Sequence  of  Presentation 

In  the  chapters  that  remain  a  detailed  examination  is  made  of  the 
methane  generation  process  and  the  SLAM  model  as  well  as  a  cost  comparison  of 
the  anaerobic  digester  propellant  generation  system  versus  the  transportion  of 
conventional  propellants  to  the  space  station.  An  extensive  review  of  current 
literature  related  to  anaerobic  digestion,  space  station  designs,  handling  of  human 
waste  products  in  space,  and  current  methods  of  orbit  maintenance,  fuel  energy, 


and  costing  calculations  are  presented  in  Chapter  II  to  provide  background  for  the 

study.  Chapter  III  provides  a  detailed  description  of  the  development  of  the 

SLAM  model  which  is  used  in  determining  methane  generation  rate  and  digester 

tank  volume.  Data  analysis  is  also  presented  in  Chapter  III.  In  Chapter  IV  the 

calculations  made  in  determining  space  station  size  and  orbit  maintenance 

requirements  are  presented.  The  best  method  for  utilizing  the  generated  methane  is 

selected  and  propellant  requirement  is  weighed  against  available  methane  output. 

Also  in  Chapter  IV,  costing  data  of  the  anaerobic  digester  and  shuttled  propellants 

are  determined  and  compared.  In  Chapter  V  a  discussion  of  the  results  and 

conclusions  are  presented  ba^ed  upon  the  results  of  Chapters  III  and  IV.  Chapter 

V  concludes  with  a  summary  of  the  study  results  and  recommendations  for  further 
study  and  action.  Appendix  A  lists  definitions  of  terms  used  throughout  this 

study.  Appendix  B  through  E  provide  an  example  of  the  SLAM  model  and 

typical  output  data.  Appendix  F  is  an  explanation  of  the  thermochemical 

calculations  made  in  Chapter  IV  to  determine  combustion  temperatures  of  methane 

with  oxygen.  Appendix  G  is  the  program  listing  including  input  and  output  for  the 

cost  analysis  done  in  Chapter  IV. 


II.  LITERATURE  REVIEW 


A  thorough  investigation  of  the  current  literature  was  undertaken  to 
provide  the  knowledge  and  background  in  the  areas  of  anaerobic  methane 
generation,  system  modeling  techniques,  feasible  space  station  concepts, 
astrodynamic  principles  with  regard  to  atmospheric  drag,  space  propulsion,  and 
economics  concerning  anaerobic  digester  components  and  space  transportation. 
This  chapter  provides  a  review  of  the  information  obtained  during  the  literature 
search  which  is  essential  in  establishing  a  foundation  for  the  study  of  space  based 
methane  generation  from  human  waste  and  its  role  as  an  orbit  maintenance 
propellar.  Although  much  research  has  been  done  in  the  area  of  waste 

management  in  space  and  in  using  space  trash  as  a  solid  propellant,  very  little 
attention  has  been  given  to  the  economic  feasibility  of  providing  a  gaseous 
propellant  from  an  environmental  control  life  support  subsystem.  This  review, 
therefore,  starts  by  identifying  approaches  that  may  offer  options  feasible  for  the 
methane  generation  process  with  regard  to  its  capabilities  aboard  a  space  station. 
A  close  look  at  the  anaerobic  process,  space  applicable  anaerobic  digester  systems, 
capacities,  and  through-puts  is  accomplished  to  provide  an  understanding  of  the 
feasibility  of  such  a  concept  for  biogas  production  in  space.  A  suitable  scenario 
specifying  space  station  size  and  weight  must  be  established  in  order  to  size  a 
particular  digester  system  to  its  load  and  determine  its  methane  output  rate. 
Therefore,  following  the  methane  generation  review  is  an  introduction  to  NASA’s 
research  work  in  the  area  of  space  settlements  and  what  can  be  expected  in  a 
near-future  LEO  space  station.  A  coupling  of  the  astronaut's  waste  in  the  space 
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station  with  the  methane  generator  is  accomplished  by  means  of  a  computer 
simulation  model  in  order  to  provide  an  accurate  estimate  of  the  methane  output 
rate  and  the  digester  volume  requirements  in  terms  of  weight  and  size.  Following 
the  methane  generation  review,  a  description  of  the  simulation  language  SLAM 
[simulation  language  for  alternative  modeling]  is  presented  as  background  for 
understanding  the  model  presented  in  Chapter  III.  The  concluding  sections  of  this 
chapter  provide  an  overview  of  information  sources  used  in  determining  digester, 
conventional  propellant  and  shuttle  transportation  costs.  (18;  36) 

The  Anaerobic  Process 

Methane  gas  can  be,  and  has  been,  produced  under  controlled  conditions 
for  many  years.  Sewage  treatment  facilities  use  anaerobic  digestion  extensively. 
To  date  many  research  organizations  both  public  and  private  have  realized  the 
potential  of  producing  heat  energy  from  organic  material  which  would  otherwise  be 
dumped,  unstabilized,  into  our  environment.  Research  efforts  began  as  far  back  as 
1939  with  the  Gobar  Gas  Plant  fabricated  at  the  New  Delhi,  India,  Agricultural 
Research  Institute.  China  currently  has  over  half  a  million  small  scale  digesters. 
India  has  installed  some  100,000  such  plants,  and  Korea  is  building  50,000 
small-scale  anaerobic  operations.  In  1976  the  U.S.  Energy  Research  and 
Development  Administration  awarded  a  research  contract  to  Waste  Management 
Inc.  of  Oak  Brook,  Illinois,  to  build  a  municipal  trash  anaerobic  digestion  facility  in 
Florida  to  further  study  biogas  production  from  municipal  refuse.  (2:95) 

While  looking  toward  the  future,  NASA  has  realized  that  permanent 
stations  in  space  will  present  new  problems,  particularly  in  minimizing  waste 
products  and  the  use  of  consumables  (e.g.  propellants  requiring  resupply  from 
Earth).  A  partial  solution  to  both  of  these  problems  may  be  found  in  the  research 
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work  that  has  been  done  on  producing  energy  from  waste  here  on  Earth.  (15;  36) 

The  primary  goal  in  the  anaerobic  process  is  to  produce  a  stable  sludge 
which  is  not  subject  to  further  biological  decomposition,  which  is  less  odorous  and 
putrescible,  and  to  reduce  the  pathological  organism  content.  The  procedures  which 
are  used  to  accomplish  these  objectives  result  in  a  gaseous  by-product  often 
referred  to  as  "sewer  gas",  "biogas",  or  "digester  gas".  This  gas  is  not  pure  methane 
but  rather  an  approximate  70/30  mix  of  methane/carbon  dioxide,  with  small 
percentages  of  nitrogen  gas,  hydrogen  gas,  hydrogen  sulfide  gas  and  water  vapor. 
Of  all  these  gasses,  only  methane  contributes  any  significant  energy  value.  Pure 
methane,  a  colorless,  odorless  hydrocarbon  is  combustible  at  concentrations  of  5  to 
12%  by  volume  in  air  and  has  an  energy  value  of  approximately  35,800  kJ/m3  of 
methane.  Since  digester  gas  is  about  70%  methane,  it  has  an  energy  value  of 
approximately  25,000  kJ/m3  By  comparison,  natural  gas  has  an  energy  value  of 
approximately  37,300  kJ/m3  Throughout  this  study  the  terms  methane  and  biogas 
will  be  used  interchangably.  In  addition  to  the  gas  produced  during  the 
stabilization  process,  a  large  percentage  of  the  sludge  is  converted  to  liquid  which 
can  be  reclaimed  through  water  recovery  subsystems.  This  conversion  of  the 
digester  sludge  to  gas  and  liquid  can  reduce  the  waste  products  of  a  space  station 
and  can  provide  both  a  gas  as  a  candidate  propellant  and  a  liquid  for  use  in  water 
recovery.  (9:35;  41:104) 

A  detailed  discussion  of  the  microbiology  and  chemistry  of  anaerobic 
digestion  is  beyond  the  scope  of  this  study  and  unnecessary  to  an  understanding 
of  the  conditions  required  for  the  anaerobic  digestion  of  organics  to  gas.  For 
in-depth  information  regarding  the  biochemical  aspects  of  anaerobic  digestion  the 
reader  is  referred  to  references  (41)  and  (49)  in  the  bibliography.  A  general 


knowledge  of  the  requirements  of  the  anaerobic  bacteria  responsible  for  making 
digester  gas  is  necessary.  Therefore,  an  overview  of  the  literature  pertaining  to 
these  bacteria  is  presented  here. 


Anaerobic  digestion  is  a  complex  biochemical  process  in  which  several 
groups  of  anaerobic  organisms  simultaneously  absorb  and  break  down  organic 
material,  human  feces,  industrial  organic  waste,  and  plant  material  in  the  complete 
absence  of  molecular  oxygen.  This  process  can  be  considered  in  two  stages  as 
shown  in  Figure  2.  In  the  first  stage  the  acid-forming  organisms  convert  the 
complex  organic  substrate,  which  is  in  a  particulate  form,  to  volatile  organic  acids 
by  attaching  themselves  to  the  particles  and  secreting  extracellular  enzymes. 
During  the  second  stage  the  methane-forming  bacteria  use  the  acids  to  produce 
carbon  dioxide  and  methane  gas.  It  is  in  the  second  stage  that  waste  stabilization 
occurs.  The  methane  is  insoluble  in  water  and  will  escape  as  a  gas.  A  means  of 
collecting  the  gas  must  be  provided  and  this  also  becomes  a  requirement  for 
pseudogravity  in  providing  a  force  for  separating  the  gas  from  the  liquid  and  the 
liquid  from  the  sludge.  (12:4-6,  4-7) 
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Figure  2.  Diagram  of  Waste  Stabilization  (12:4-6) 


The  methane  bacteria  growth  rate  is  relatively  slow  and  ranges  in  time 
from  2  to  22  days  depending  on  pH,  substrate  composition  and  temperature.  The 
pH  range  tolerance  for  these  bacteria  is  between  6.5  and  7.6.  At  a  given 
temperature,  if  pH  drops  below  6.0,  methane  production  stops  and  there  will  be  no 


decrease  in  the  organic  content  of  the  sludge  and  no  methane  produced. 

Compared  to  the  acid  forming  bacteria  the  methane  bacteria  get  very 
little  energy  from  their  food.  Because  of  this  they  are  fewer  in  numbers  and  are 
more  sensitive  to  changes  in  temperature.  Most  digesters  operate  at  temperatures 
between  26  and  43  °C,  also  known  as  the  mesophilic  range.  Operating 
temperatures  in  the  45  to  55  °C,  or  thermophilic  range,  have  have  been  investigated 
and  show  a  slight  increase  in  digester  efficiency,  but  are  more  sensitive  to 
temperature  fluctuations  and  require  more  energy  to  operate  at  this  higher 
temperature.  (41:18,  88) 

Based  on  the  above  criteria  four  conditions  are  identified  as  being 
essential  for  efficient  digester  operation  and  thus  maximum  methane  production. 
These  are: 

1.  An  environment  free  of  molecular  oxygen. 

2.  A  steady  temperature  in  the  mesophilic  range  [26  to  46  °C]. 

3.  A  proper  pH  between  6.5  and  7.6. 

4.  Sufficient  retention  time  to  allow  methane  bacteria  regeneration. 

When  discussing  the  subject  of  temperature,  the  element  of  time  cannot 
be  ignored  because  solids  stabilization  cannot  be  accomplished  at  low  temperatures 
unless  sufficient  time  is  allowed.  Temperature  and  retention  time  for  efficient 
digester  operation  are  shown  in  Table  I  to  illustrate  this  point.  (12:4-13;  41:95) 


Suggested  Solids  Retention  Time  for  High  Rate  Digesters  (41:95) 

Operating 

Temperature  Suggested  Time 
[°C]  [Days] 

18  28 

24  20 

30  14 

35  10 

40  10 

Methane  Generation 

There  are  two  general  digester  processes  which  have  evolved  over  the 
years.  These  two  processes  are  referred  to  as  standard  rate  digestion  and  high 
rate  digestion.  These  two  terms  are  somewhat  deceptive  in  that  they  do  not  refer 
to  the  rate  at  which  organic  materials  are  converted  to  methane.  However,  the 
high  rate  digester  does  provide  conditions  which  are  more  favorable  to  the 
anaerobic  processes.  From  these  two  processes,  variations  in  staging  and  mixing 
have  resulted  in  several  designs  which  will  be  considered  here  in  selecting  a 
digester  configuration  which  is  appropriate  for  a  space  station.  (48:255) 

In  a  standard  rate  digestion  system  the  contents  of  the  digester  are 
unmixed  and  the  processes  of  digestion,  sludge  thickening,  and  supernatant 
formation  are  carried  out  in  stratified  layers  of  the  tank.  As  a  result,  actively 
digesting  sludge  occurs  only  in  a  portion  of  the  total  digester  volume  reducing  the 
actual  working  volume  of  the  digester  tank.  Because  the  tank  is  unmixed  the 
feeding  of  the  anaerobic  bacteria  is  a  very  slow  process  taking  between  30  to  60 


days  retention  of  the  organic  solids.  The  long  solids  retention  time  requires  a 
large  tank  volume  and  a  corresponding  fresh  sludge  input  rate  usually  in  periodic 
steps  of  two  to  three  times  a  day.  (41:89;  48:255) 

The  high  rate  digestion  process  differs  from  the  standard  rate  process 
in  that  the  sludge  is  mixed  and  active  sludge  digestion  takes  place  throughout  the 
tank.  This  mixing  provides  improved  heating  and  continuous  feeding,  allowing  for  a 
more  complete  interaction  of  the  microorganisms  with  the  organic  sludge  material. 
The  mixing  process  requires  a  nearly  continuous  sludge  input  to  the  digester  and 
therefore  high  rate  refers  to  the  organic  loading  rates  possible  rather  than  the 
methane  generation  rate.  The  incoming  sludge  displaces  the  digested  sludge  in  the 
tank,  therefore  tank  size  is  determined  by  how  long  the  organic  solids  must  remain 
in  the  tank  [solids  retention  time],  which  is  a  function  of  temperature  [Table  I], 
(48:255) 

The  high  rate  digestion  process  has  the  advantages  of  shorter  solids 
retention  time,  continuous  feed,  smaller  tank  volume,  and  uniform  heating 
throughout  the  tank  due  to  mixing.  Also  because  there  is  no  supernatant 
separation,  the  total  solids  are  reduced  by  45  -  50%  and  given  off  as  gas.  These 
factors  make  the  high  rate  digester  process  a  natural  choice  for  the  space  station. 
Because  there  is  no  supernatant  separation  in  the  high  rate  digester  a  second  tank 
must  be  provided  to  allow  settling  of  the  digested  sludge  from  the  liquid.  The 
second  tank  allows  storage  of  the  digested  sludge  until  it  is  concentrated  at  the 
bottom  and  can  be  drawn  off  for  subsequent  conventional  space  vacuum  drying. 
The  supernatant  can  then  be  piped  to  the  on-board  water  recovery  system.  Figure 
3  is  a  simplified  diagram  of  a  high  rate  two  stage  digester.  Both  tanks  are 
covered,  and  the  gas  collection  system  cross  channeled  between  them.  The  second 


tank  will  not  produce  much  gas  since  most  of  the  gas  production  takes  place  in  the 
first  tank.  The  second  tank  does  contain  a  large  volume  of  good  active  methane 
producing  bacteria  which  can  be  used  to  seed  the  first  tank  if  the  digestion  process 
slows  down  due  to  fluctuations  in  temperature  or  feeding  rates.  (12:4-11;  32:82; 
35:922-923;  41:92) 


Figure  3.  High  Rate  Two  Stage  Digester  (12:4-10) 

Mixing  of  the  primary  tank  increases  the  volatile  solids  breakdown  and 
increases  the  amount  of  gas  produced.  Mixing  can  be  accomplished  by  artificial 
means  both  internal  and  external  to  the  tank.  Internal  mixing  can  be  accomplished 
through  impellers  or  turbine  wheels  submerged  within  the  tank.  An  important 
factor  to  consider  in  this  method  of  mixing  is  exposure  of  the  mixing  blades  to  grit 
and  debris.  This  can  cause  wear  of  the  mixing  impellers  which  will  require 
shutdown  of  the  digester  for  drainage  and  repair.  A  more  acceptable  means  of 
mixing  can  be  accomplished  external  to  the  tank  through  the  use  of  a  pump.  The 
total  capacity  of  the  pump  is  generally  less  than  the  circulating  capacity  of  the 
mixers.  The  pump  also  allows  for  the  external  heating  of  the  tank  by  cycling  the 
sludge  through  heat  exchangers.  By  using  external  heating  and  mixing,  equipment 


failures  can  readily  be  handled  without  stopping  the  digestion  process  by  switching 
the  sludge  recirculation  flow  to  redundant  equipment.  (12:4-28) 

The  cover  of  the  digester  has  some  unique  features  which,  though  not 
critical  to  the  results  of  this  study,  should  be  recognized  here.  The  two  types  of 
covers  are  the  fixed  and  floating  cover.  The  fixed  cover  offers  the  advantages  of 
simplicity  in  design.  It  must  be  equipped  with  a  positive  displacement  feed  and 
draw-off  arrangement  so  that  a  negative  gas  pressure  in  the  tank  does  not 
develop.  If  this  were  to  happen,  atmospheric  air  within  the  space  craft  could  be 
drawn  into  the  tank  and  mix  with  the  methane  producing  an  explosive  mixture. 
The  explosive  limits  of  methane  in  air  are  between  5.3  and  14  percent.  This  means 
that  a  methane:air  ratio  between  as  little  as  5.3  to  14  percent  could  present  a 
hazardous  situation.  A  floating  cover  design  provides  the  greatest  degree  of  safety 
since  air  cannot  be  drawn  in  by  a  negative  pressure.  The  problems  with  this  type 
design  are  maintaining  the  cover  guides  in  a  smooth  operating  condition  and 
keeping  the  cover  level.  In  either  case  the  gas  will  be  compressed  and  stored 
external  to  the  digester  to  minimize  any  leakage  and  possible  safety  hazards. 
(12:29;  48:260) 

The  expected  methane  output  from  the  anaerobic  process  varies  from 
author  to  author.  A  popular  figure  based  on  actual  data  is  about  1  cubic  foot  per 
capita  per  day.  Anderson  (2)  cites  a  figure  of  4  -5  ft3  for  every  pound  of  organic 
matter  destroyed.  Price  (41)  suggests  the  use  of  Michaelis-Menton  type  rate 
equations  to  determine  quantity  of  methane  generated  based  on  biological  oxygen 
demand  [BOD]  of  the  sludge,  the  quantity  of  biological  solids  added  per  day,  and 
the  efficiency  of  waste  utilization.  These  equations  are  used  in  this  study  and 
coupled  to  the  expected  BOD  and  input  loads  rate,  based  on  NASA  studies,  using 


simulation  modeling.  The  modeling  and  calculations  of  input  parameters  comprise 
a  majority  of  this  study  and  Chapter  III  describes  the  equations,  the  SLAM  model 
and  the  input  data  used  to  determine  the  methane  generation  rate  and  digester 
volume.  (2:95;  41:37-39,  94,  105-106,  121;  48:111) 

The  Scenario 

In  considering  the  future  utility  of  space  as  a  resource  and  as  a 
laboratory,  it  is  obvious  that  to  fully  exploit  the  potentials  of  space  we  must  have 
a  continuous  manned  presence  in  this  new  frontier.  NASA  has  been  the 
forerunner  in  looking  toward  man’s  future  in  space,  and  has  sponsored  several 
studies  with  the  expert  assistance  of  such  space  pioneers  as  Gerard  O’Neill,  which 
have  put  numbers  to  ideas  in  assessing  the  human  and  economic  implications  as 
well  as  the  technical  feasibility  of  settlements  in  space.  As  man  moves  toward 
establishing  colonies  in  space  one  of  the  first  steps  that  must  be  taken  is  the 
siting  of  small  stations  in  low  earth  orbit  for  research,  development,  demonstration, 
testing  and  evaluation.  Low  earth  orbit  is  necessary  as  a  first  location  because 
materials  must  come  from  the  Earth  during  the  initial  stages  of  space  colonization, 
and  transportation  to  this  orbit  by  means  of  the  Shuttle  is  much  more  technically 
and  economically  feasible  than  going  directly  to  deep  space  orbits  such  as  the 
Lagrangian  liberation  point  L6  or  the  Moon.  The  LEO  space  station  will  be 
capable  of  performing  research  in  materials  fabrication  and  assembly  techniques, 
solar  and  nuclear  power  generation  systems,  and  physiological  effects  of  rotation 
and  reduced  gravity.  (34:v,  1) 

This  study  utilizes  NASA’s  research  as  a  basis  for  analyzing  the 
feasibility  of  space-based  methane  genera  ;on.  AH  figures  relating  to  the  size, 
shape,  capacity  and  orbital  altitude  of  the  LEO  space  station  come  from  research 


work  done  by  NASA  and  published  in  references  (34)  and  (30)  listed  in  the 
Bibliography. 

Requirements  for  pseudograyity  are  also  based  on  these  studies,  and 
past  NASA  research,  which  show  the  physiological  consequences  of  a  gravity  free 
environment  over  a  sustained  length  of  time.  Data  from  past  space  flights  show 
that  decalcification  occurs  at  a  rate  of  1  to  2%  per  month  in  the  absence  of 
gravity,  which  can  result  in  decreased  bone  mass  and  density.  Other  effects  of 
gravityless  environment  include  hormone  and  electrolyte  imbalances  and  unstable 
protein  and  carbohydrate  states.  For  these  reasons  pseudogravity  must  be 
provided  for  the  people  of  the  space  station.  The  most  feasible  way  of  generating 
artificial  gravity  is  to  rotate  the  space  station.  However,  a  rotating  system 
contains  forces  other  than  the  centrifugal  force  which  acts  as  the  gravity.  The 
coriolis  force  is  one  of  the  forces,  caused  by  the  speed  of  motion  and  its  direction 
relative  to  the  axis  of  rotation.  Consequently,  motion  sickness  can  result  even  at 
low  rotation  rates  due  to  the  cross-coupled  angular  accelerations  in  the  ear  when 
the  head  is  turned  out  of  the  rotation  plane.  People  can  adapt  to  rates  below  3 
rpm  and  for  this  reason  the  space  station  is  designed  foT  a  rotation  of  2  rpm. 
Lower  rotation  rates  are  preferable,  but  for  a  platform  of  small  size  this  is  not 
feasible  due  to  the  large  radius  of  rotation  required.  Figure  4  depicts  the  radius  of 
rotation  required  for  various  spin  rates.  (30:36,  40;  34:21,  22;  35:349-408,  37:154) 
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presented  in  Table  II.  Note  that  an  atmosphere  equivalent  to  one-half  of  the 
Earth’s  atmosphere  is  considered  in  uhe  structural  design  to  determine  the  mass  of 
the  station.  By  using  a  lower  atmospheric  pressure  inside  the  space  staion,  less 
structural  strength  is  required  to  hold  the  structure  together  against  the  vacuum  of 
space.  This  is  accomplished  by  increasing  the  amount  of  oxygen  in  the  space 
craft  while  diluting  it  with  an  inert  gas,  such  as  helium  or  nitrogen,  thereby  bringing 
the  total  pressure  to  .5  atmosphere.  (34:41,  144) 


Figure  5.  Dumbell  and  Beaded  Torus  Design 


TABLE  II.  Space  Station  Parameters  (34:46,  87,  147) 


Radius  of  Rotation  [m]  236 

Radius  of  Sphere  [m]  33.3 

Total  Mass  [kg]  1.45  x  106 

Structural  Mass,  1/2  atmosphere  [kg]  3.63  x  10s 

Habitat  Mass,  4.5  xlO*  kg/person  [kg]  5.07  x  10B 

Atmospheric  Mass  [kg]  1.81  x  10B 

Orbital  Altitude  [km]  240 

Population  200 

Area/Person  [mJ]  35 


SLAM  Model 


Simulation  modeling  provides  a  means  of  representing  a  system  in  a 
mathematical  form  so  that  it  can  be  exercised  on  a  digital  computer.  SLAM  was 
chosen  as  the  simulation  language  to  be  used  in  representing  the  anaerobic 
digestion  process  in  this  study  because  of  it’s  capabilities  to  perform  continuous 
simulation  modeling.  In  this  type  of  modeling  the  state  of  the  system,  in  this  case 
the  volatile  solids  loading  rate  and  methane  generation  rate,  are  represented  by 
dependent  variables  which  change  over  time.  This  time-dependent  portrayal  of  the 
variables  provides  an  output  which  can  be  analyzed  to  provide  a  realistic  estimate 
of  system  performance.  Chapter  III  explains  how  SLAM  was  used  to  generate 
time  varying  statistics  on  parameters  such  as  methane  generation,  organic  solids 
loading,  and  digester  capacity.  Detailed  information  on  the  syntax  and  operation 
of  the  SLAM  program  can  be  found  in  Pritsker’s  text,  Introduction  to  Simulation 
and  SLAM  //  (42) 

Cost  Analysis 

The  cost  analysis  used  for  this  study  has  to  be  be  flexible  enough  to 
allow  for  a  technological  uncertainty  in  pricing  the  components  of  a  system  that 
will  be  procured  and  operated  many  years  in  the  future.  It  also  must  allow  an 
input  of  the  expert's  "gut  feeling"  or  past  experience  in  determining  where  the 
actual  price  will  fall  in  the  uncertainty  range.  The  Monte  Carlo  method  is  selected 
to  be  used  as  the  analysis  tool  for  forcasting  the  life  cycle  cost  of  the  digester 
system  and  of  the  shuttled  propellants,  upon  which  the  economic  feasiblility  of  the 
systems  in  question  will  be  based.  This  method  is  recommend  in  the  studies, 
Military  Equipment  Cost  Analysis  (43),  and,  Estimating  Cost  Uncertainty  Using 
Monte  Carlo  Techniques  (7),  done  by  the  Rand  Corporation.  The  Monte  Carlo 
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Technique  requires  that  the  expert  express  his  opinion  of  the  uncertainty  in  the 
price  range  selected  by  picking  a  beta  distribution  curve  from  a  family  of  nine 

which  best  discribes  the  variance  of  the  price  based  on  his  experience.  More 

information  on  the  Monte  Carlo  Techniqe  can  be  found  in  references  (13)  and  (27). 
(7;  43) 

The  pricing  of  the  components  in  the  digester  system  requires  U:at  a 
system  configuration  be  selected  so  that  individual  elements  of  the  system  can 
each  be  analized  to  determine  a  price  range  and  uncertainty  determined  from  its 
application  and  environment.  The  alternative  system,  shuttled  conventional 
propellants,  must  also  be  analized  to  determine  propellant  costs.  Each  system,  in 

addition  to  these  aquisition  costs,  will  have  an  operational  cost  which  will  require 

the  estimation  of  future  space  shuttle  transportation  costs.  A  preview  of  the 
source  of  such  costs  is  described  below. 

Digester  Costs.  Once  system  configuration  is  specified,  the  major 
components  of  the  system  can  be  identified  and  a  cost  value  estimated  for  each. 
Consultation  with  waste  treatment  engineers  with  experience  in  design  of  various 
Environmental  Protection  Agency  projects  and  municipal  waste  treatment  facilities 
is  the  main  source  of  information  in  costing  the  primary  elements  of  the  digester 
system.  Due  to  the  technological  uncertainty  in  the  application  of  equipment  used 
strictly  on  Earth  for  anerobic  digesters  to  a  space  station  environment,  both 
equipment  cost  and  weight  must  be  determined  using  a  Monte  Carlo  Technique. 

Transportation  costs  are  associated  strictly  with  the  initial  transport  of 
equipment  to  the  space  station  and,  in  the  case  of  using  the  biogas  as  a 
bipropellant,  the  transport  of  an  oxidizing  agent  suitable  for  combustion  with  a 
methane/carbon  dioxide  mixture.  These  transportation  costs  are  considered  in  the 


life  cycle  cost  analysis. 


Shuttled  Propellant  Costs.  The  alternative  to  generating  a  biogas 
propellant  from  human  waste  is  to  transport  conventional  propellants  on  a  regular 
basis  to  the  space  station.  The  life  cycle  cost  of  such  an  alternative  must  be 
compared  to  similar  costs  of  the  digester  biogas  generation  system  to  determine  if 
an  economic  advantage  lies  with  either  system.  Pricing  information  on 
conventional  propellants  will  be  obtained  from  personnel  in  logistics  centers 
responsible  for  purchasing  propellants  for  NASA  and  the  DoD.  Transportation 
costs  will  reflect  the  quantity  of  conventional  propellants  required  to  produce  the 
equivalant  capabilities  of  the  digester  system.  These  costs  will  be  based  on 
propellant  weight  based  on  density. 


III.  SLAM  MODEL  FORMULATION 


The  purpose  of  the  model  is  to  generate  data  for  analysis  since  actual 
data  on  small  scale  digesters  applicable  to  a  space  station  scenario  are  not 
available.  Factors  to  be  considered  which  effect  the  anaerobic  process  are  the 
volatile  solids  retention  time  [SRT]  in  the  digester  and  temperature.  Equations 
used  in  formulating  the  model  come  from  Price  and  Cheremisinoff’s  book,  Biogas 
Production  and  UUitzatton  (41)  and  are  presented  below. 

VCH4  =  [0.35  m3/kg][[EQSo[103g/kg]~']-1.42Px]  [1] 

and, 


P,  =  [YQ[ESo][103g/kg]-‘]/[KkdT]  [2] 


where; 

V cH4  =  volume  of  methane  produced,  m3/d 
0.35  =  theoretical  conversion  factor  for  the  amount  of  methane 
produced  from  the  conversion  of  1  kg  of  BODL 
E  =  efficiency  of  waste  utilization  [  8  -  .95] 

Q  =  flow  rate,  m3/d 

S0  =  organic  material  added  [BODL]  g/m3 
1.42  =  conversion  factor  for  volatile  biological  solids  to  BODL 
Px  =  net  mass  of  volatile  biological  solids  [cell  tissue]  produced 
per  day,  kg/d 


Y  =  growth  yield  coefficient 

kd  =  microorganism  decay  coefficient 

T  =  mean  cell  residence  [solids  retention]  time 


Lawrence  and  McCarty  (23)  have  studied  the  growth  kinetics  of  the 
methane  producing  microorganisms  and  have  determined  the  anaerobic  process  to 
be  rate  limited.  Their  mode!,  which  describes  the  net  growth  rate  of  the 
microorganisms  in  a  continuous  flow  completely  mixed  anaerobic  treatment  system, 
is  represented  by  the  expression: 


dX/dt  =  Y[dF/dt]  -  kdX 


where; 


dX/dt  =  microorganism  net  growth  rate  per  unit  volume  of  digester 
dF/dt  =  rate  of  waste  utilization  per  unit  volume  of  digester 
X  =  microorganism  concentration 
Y  =  growth  yield  coefficient 
kd  =  microorganism  decay  coefficient 
Figure  6  contains  the  graph  from  Price  and  Cheremisinoff’s  book  used 
in  determining  the  microorganism  growth  yield  coefficient  [Y],  and  decay 
[endogenous]  coefficient  [kd],  This  is  a  typical  plot  from  laboratory  data  in  which 

the  microorganism  growth  rate  divided  by  the  microorganism  concentration  [dX/dt 
•f  X,  which  is  the  reciprocal  of  the  solids  retention  time  ]  is  plotted  against  the 


rate  of  waste  utilization  divided  by  the  microorganism  concentration  [dF/dt  -r  X], 

The  growth  yield  coefficient  is  then  the  slope  of  this  line  and  the  decay  coefficient 
is  the  intercept. 
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(dX/dt)/X 


(dF/dt)/X 

Figure  6.  Laboratory  Data  Plot  (41:107) 


Values  of  Solids  Retention  Time  [SRT]  which  allow  for  maximum  use 
of  the  waste  solids  by  the  methane  bacteria  are  suggested  by  Price  and 
Cheremisinoff  for  various  temperatures  [Table  I,  page  19]  and  have  been  plotted  in 
Figure  7.  For  the  mesophilic  bacteria  only  temperatures  in  the  30  to  40  C 
temperature  range  must  be  considered. 


ime 

Days) 


^  Temp  (C) 

Figure  7.  Recommended  SRT  for  Various  Temperatures.  (41:95) 


This  graph  provides  two  equations  for  the  SLAM  model  which  will  then 
determine  the  appropriate  solids  retention  time  for  a  given  temperature.  For  the 


temperature  range  from  30  to  35  G,  solids  retention  time,  T,  is: 

T  =  -8[TEMP]  -  38  days 

For  temperatures  between  35  and  40  C,  T  is  a  constant  10  days. 

P'actors  which  also  influence  the  solids  retention  time  are  the  volatile 
solids  loading  on  the  digester,  and  the  percentage  concentration  of  volatile  solids 
suspended  in  the  raw  sludge.  For  this  study  the  percent  solids  in  the  raw  sludge 
is  13%,  based  on  the  quantities  of  waste  expected  in  a  typical  space  station.  A 
curve  showing  the  relationship  among  solids  loading,  solids  retention  time,  and 
sludge  solids  is  presented  in  Figure  8.  From  this  curve  a  constant  of  multiplication 
of  13  can  be  extracted  for  use  in  determining  volatile  solids  loading  from  solids 
retention  time  for  use  in  the  digester  model.  (20:80;  41:97,  99) 


Figure  8.  Volatile  Solids  Loading  vs  SRT  (41:99) 


Model  Development 


A  continuous  model  was  formulated  to  simulate  the  operational 
characteristics  of  the  anaerobic  digester.  The  program  is  written  to  allow 
flexibility  of  input  parameters  and  monitoring  of  various  internal  factors  necessary 
for  model  verification.  Figure  9  depicts  the  interaction  of  the  various  parameters 
and  references  the  equation  or  figure  used  in  determining  the  structure  of  the 


DETERMINATION  OF  METHANE  PRODUCED, 
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A  detailed  description  of  subroutine  STATE  is  now  presented.  It  is 
included  in  the  main  program  [Appendix  B]  and  each  fundamental  line  is  listed 
here  followed  by  a  brief  explanation  of  it’s  contents  and  origin. 
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C  EFFICIENCY: 

SS(7]  =  RNORM  (.875,  .075,  2) 

Allows  for  random  variation  of  system  efficiency  within  .8  -  .95 
range  (41:121) 

C  TEMP  CONSTANT: 

C  SS(ll)  =  XX(1) 

C  TEMP  VARIATION: 

SS(ll)  =  XX(1)-2.8*COS(8*PI*(TNOW-0.1042)) 

Allows  for  test  runs  with  either  a  constant  temperature  or  with  a 
temperature  which  simulates  cooling  of  digester  due  to  cooler  waste 
inputs.  The  cycling  of  temperature  lags  behind  the  loading  cycle  by 
15  minutes  per  hour  to  simulate  heat  transfer  due  to  mixing. 

C  LOAD  CONSTANT: 

C  SS(12)  =  2.6*XX(2)*TNOW 

C  LOAD  VARIATION: 

SS(12)  =  XX(2)*2.6*(1+0.4*COS(8*PI*TNOW)) 

Allows  for  test  runs  with  either  a  constant  load  or  a  loading  which 
simulates  six  hour  input  cycles  which  would  follow  the  normal 
activity  cycles  [prime  waste  facility  use]  of  a  crew  working  shifts 
over  a  24  hour  period.  XX(2)  is  an  input  variable  and  represents 
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people  on  the  space  station.  The  total  waste  expected  per  person  in 
a  space  station  is  2.6  kg/person/day  (20:80). 

FLOW: 

SS(1)  =  SS(12)/100 

Converts  load  input  [kg/day]  to  flow  [m3]. 

SOLIDS  RETENTION  TIME: 

SS(2)  =  10 

IF  (SS(ll)  .LT.  35)  SS(2)  =  -0.8*SS(ll)+38 

Determines  recommended  Solids  Retention  Time  in  digester  by  using 
temperature  and  equations  derived  from  Figure  7.  Figure  7  is  derived 
from  tabled  data  in  Price  and  Cheremisinoff’s  book  (41:95). 

TOTAL  SOLIDS  (1356): 

SS(3)  =  0.13*SS(12) 

Total  of  solids  in  raw  sludge  input  are  determined  to  be  13%  of  the 
weight  of  the  raw  sludge  (20:80). 

VOLATILE  SOLIDS  CALCULATED  (-3%ASH)  (UNITS:KG/D): 

SS(4)  -  SS(3)-0.03*SS(3) 

The  volatile  solids  are  calculated  by  subtracting  3%  ash  which  has 
been  experimentally  determined  to  be  the  non-organic  residue  in  the 
waste.  (35:922) 

METHANE  PRODUCED: 

SS(5)  =  0.35E*SS(4)*SS(7)*(l-(0.0588/(l-0.033*SS(2)))) 

This  equation  for  the  volume  of  methane  produced  is  a  combination 
of  equations  [1)  and  [2)  and  combines  the  constants  Y  and  kd 


derived  from  Figure  6  [slope  and  intercept  respectively)  with  various 


conversion  factors. 

SOLIDS  REMOVED: 

SS(6)  =  SS(3)-1.059*SS(5) 

This  line  of  code  provides  information  on  the  amount  of  waste 
remaining,  after  the  anaerobic  process,  that  will  be  removed  from  the 
digester  and  returned  to  Earth  or  could  possibly  be  utilized  as  a 
fertilizer  in  plant  growth  experiments.  It  is  calculated  by  subtracting 
the  mass  of  the  methane  from  the  input  solids.  The  figure  1.059  is  a 
methane  density  conversion  factor  to  convert  methane  volume  to 
weight. 

VOLATILE  SOLIDS  LOADING  (FIGURE  8)  (UNITS:KG/M3/D): 

SS(8)  =  13/SS(2) 

This  equation  uses  the  Solids  Retention  Time  to  determine  the 
volatile  solids  loading  per  unit  volume  of  digester.  The  equation  was 
derived  from  Figure  8.  The  constant  of  multiplication  [k  =  SRT  x 
VS]  is  13. 

DIGESTER  VOLUME  USED: 

SS(9)  =  SS(4)/SS(8) 

This  calculation  determines  digester  volume  from  the  volatile  solids 
present  and  the  required  volatile  solids  loading  per  unit  volume  of 
digester  calculated  above. 

HYDRAULIC  RETENTION  TIME: 

SS(10)  =  SS(9)/SS(1) 

This  calculation  determines  how  long  the  waste  influent  will  remain  in 
the  digester. 


The  SLAM  input  program  is  written  to  produce  an  output  graph  and 
time  averaged  statistics  on  parameters  which  can  be  compared  to  theoretical 
values  in  verifying  the  model.  The  input  statements  are  listed  in  Appendix  B  and 
a  sample  output  run  is  included  in  Appendix  C. 

Experimentation 

Model  validation  and  verification  requires  a  set  of  input  statements 
which  allows  for  several  runs,  providing  output  of  average  methane  generated  and 
average  digester  volume  used  for  each  run.  This  new  set  of  statements  is  provided 
in  Appendix  D.  Three  input  temperatures  were  specified  and  ten  runs  each 
performed  with  system  efficiency  being  specified  by  antithetic  random  numbers  to 
reduce  variance.  The  output  of  the  thirty  runs  is  summarized  in  Appendix  E. 

The  design  of  the  verification  phase  provides  a  method  of  obtaining  the 
grand  mean  of  methane  generated  for  comparison  against  published  theoretical 
estimates.  Average  digester  volume  used  is  also  observed  to  see  if  it  varies 
between  runs  at  set  input  temperatures.  Since  digester  volume  is  a  function  of  the 
required  volatile  solids  loading  compared  to  the  actual  loading,  which  is  constant 
when  averaged  over  time,  no  change  is  expected  in  the  volume  of  the  digester  when 
operated  at  a  specific  input  temperature.  Changes  would  be  expected  though  when 
the  initial  operating  temperature  is  changed,  which  would  determine  a  new  solids 
retention  time  from  Figure  7,  and  therefore  affect  the  required  volatile  solids 
loading,  Figure  9. 


Validation  provides  a  95%  confidence  interval  in  which  average  methane 
generation  can  be  expected  to  be  valid  for  each  of  the  three  system  operating 
temperatures.  A  95%  confidence  interval  is  also  generated  on  the  mean  difference 
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between  the  three  systems  to  show  that  the  difference  in  the  amount  of  methane 
generated  is  due  solely  to  the  change  in  temperature  and  not  the  randomness  of 
the  system  efficiency. 


Output  Analysis 

The  following  table  summarizes  the  grand  mean  methane  generation 
rate,  variance,  and  95%  confidence  intervals  for  each  of  the  three  initial  operating 
temperatures.  The  data  from  the  thirty  runs  used  to  obtain  the  grand  means  is 
listed  in  Appendix  E. 


SYSTEM  TEMPERATURE  =  32.8G 


Grand  Mean  Methane  Generation  Rate:  18.1402  ms/day 


Variance: 

Confidence  Interval; 


0.0152 

18.0521  to  18.2283  m3/day 


SYSTEM  TEMPERATURE  =  35C 

Grand  Mean  Methane  Generation  Rate:  18.2424  m3/day 


Variance: 

Confidence  Interval: 


0.0213 


18.1381  to  18.3467  m3/day 


SYSTEM  TEMPERATURE  =  37.2C 

Grand  Mean  Methane  Generation  Rate:  18.2830  m3/day 


Variance: 

Confidence  Interval: 


0.0214 

18.1784  to  18.3756  m3/day 


The  methane  generation  rate  of  approximately  18.2  m3/day  agrees  with 
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observed  methane  generation  rates  of  4  -  5  ft3/lb  organic  matter,  given  by 
Anderson  (2:95).  In  this  model  an  average  of  65.6  kg/day  [144.3  Ib/day]  of  organic 
solids  are  available.  This  corresponds  to  a  theoretical  value  of  577.2  -  721.5 
ft3/day  or  16.4  -  20.43  m3/day. 

The  difference  between  methane  generation  rate  means  for  the  three 
operating  temperature  systems  is  tabulated  below  for  each  run  together  with  the 
grand  mean  difference,  variance  and  95%  confidence  interval  for  the  true  mean 
difference.  The  small  variance  in  difference  mean  is  due  to  the  change  in  system 
operating  temperatures  and  not  the  randomness  of  the  modeled  efficiency.  System 
A  operates  at  32.8C,  system  B  at  35C,  and  system  C  at  37.2C. 
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DIFFERENCE  IN  MEAN  METHANE  GENERATION  RATES 
FOR  THREE  OPERATING  TEMEPERATURES 


Run 

Di  FA-B1 

Di  TA-Cl 

DifB-Cl 

1 

-0.102 

-0.143 

-0.041 

2 

-0.101 

-0.141 

-0.040 

3 

-0.104 

-0.145 

-0.041 

4 

-0.102 

-0.142 

-0.040 

5 

-0.100 

-0.139 

-0.039 

6 

-0.100 

-0.144 

-0.041 

7 

-0.102 

-0.143 

-0.041 

8 

-0.103 

-0.144 

-0.041 

9 

-0.104 

-0.145 

-0.041 

10 

-0.103 

-0.144 

-0.041 

ft 
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GRAND  MEAN  OF  DIFFERENCE  OF  MEANS 


Mean 


-0.1021 

0.0015 


-0.1430 

0.0019 


-0.0406 

0.0007 


CONFIDENCE  INTERVAL  (95%)  FOR  TRUE  MEAN  DIFFERNCE 
A-B:  -0.1010  to  -0.1032 
A-C:  -0.1416  to  -0.1444 
B-C:  -0.0401  to  -0.0411 


The  increases  in  methane  generation  due  to  the  increase  in  operating 
temperature  is  extremely  small  and  it  may  not  be  economically  beneficial  in  space 
to  provide  the  extra  energy  required  to  increase  the  methane  generation  by 
operating  in  the  upper  mesophilic  temperature  range.  Digester  volume  used 
remained  constant  for  all  three  operating  temperatures  as  was  expected,  and  did 
decrease  for  increasing  temperatures  due  to  the  increased  metabolism  of  the 
methane  producing  bacteria  at  higher  temperatures.  Due  to  fluctuations  in  input 
loading,  the  max  volume  used  [70.6  ms]  will  be  the  design  parameter  in  sizing  the 
digester  tank. 
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IV.  ANALYSIS  AND  DISCUSSION 

Up  to  this  point  the  anaerobic  process  has  been  described  and  modeled 
as  it  applies  to  a  possible  future  space  station  scenario.  Biogas  output  has  been 
computed  for  a  given  space  station  population.  In  this  chapter  space  station  size 
will  be  determined  based  on  the  population  selection  and  previous  studies  done  by 
NASA  on  the  future  colonization  of  space.  The  mass  and  cross  sectional  area 
will  be  calculated  so  that  an  atmospheric  drag  can  be  determined  for  the  altitude 
at  which  the  space  station  will  be  orbiting.  This  atmospheric  drag,  which  is  an 
acceleration  causing  the  space  craft  to  spiral  towards  Earth,  will  determine  what 
force  must  be  provided  by  the  attitude  control  system  to  maintain  the  orbit. 

The  use  of  the  generated  biogas  as  a  cold  gas  propellant  and  as  a  fuel 
in  a  bipropellant  rocket  system  is  investigated  in  this  chapter  to  determine  which 
system  would  best  provide  the  force  necessary  to  overcome  the  atmospheric  drag. 
The  theoretical  performance  of  each  propulsion  system  is  evaluated  using 
thermochemical  analysis  techniques  to  determine  reaction  temperatures. 

The  remainder  of  the  chapter  compares  the  biogas  generation  system 
performance  with  conventional  propulsion  performance  based  on  the  orbit 
maintenance  requirement  and  what  portion  of  this  requirement  the  biogas  system 
can  meet.  The  common  reference  used  for  comparison  is  the  dollar.  Cost 
estimates  of  the  digester  system  are  determined  by  the  total  acquisition  cost, 
installation  cost  and  operational  cost,  all  amortized  over  an  assumed  lifetime  of 
thirty  years.  This  cost  is  then  compared  to  the  cost  of  providing  the  required 
attitude  control  with  conventional  propellant  delivered  by  the  Shuttle.  All  cost 
estimates  are  determined  using  a  Monte  Carlo  technique. 


Space  Station  Population  and  Size 

As  stated  previously  in  Chapter  II,  the  space  station  population  has 
been  selected  based  on  previous  work  done  by  NASA.  A  population  of  200 
occupants  was  used  in  determining  the  output  of  the  digester  system.  This  figure 
is  also  used  in  determining  the  size  of  the  space  station.  A  brief  survey  of  the 
literature  indicates  that  for  a  small  space  station  with  this  population  a  projected 
living  area  per  person  of  35  m2  is  sufficient  and  is  equivalent  to  that  found  in 
some  small  French  villages.  As  a  comparison  with  an  American  city,  the 
Manhattan  Borough  in  New  York  has  a  per  capita  area  of  38.2  m2/person. 
Projected  area  is  area  of  the  largest  plane  perpendicular  to  the  direction  of  the 
pseudogravity.  Actual  usable  area  can  be  made  larger  than  the  projected  area  by 
constructing  levels  within  the  habitat.  (34:24-25,47) 

With  a  population  of  200  and  a  required  area  of  35  m2/person,  a  space 
station  of  the  dumbell  design  will  have  those  parameters  listed  in  Table  II,  page  26. 
These  parameters  are  used  to  determine  the  acceleration  due  to  drag  on  the  space 
station.  The  formula  for  the  acceleration  due  to  drag  is  (3:423): 

ad  =  CD[l/2]p[A/m]v2  [3] 

where: 

CD  =  Coefficient  of  Drag 
p  =  atmospheric  density 
A  =  average  cross  sectional  area 
m  =  mass  of  space  station 
v  =  velocity  of  space  station 

The  average  cross  sectional  area  of  the  space  station,  A,  is  different  from  the 
total  cross  sectional  area  due  to  the  rotation  necessary  to  provide  a  psuedogravity. 


It  is  calculated  from  the  formula  (21:16); 


A  =  [2/rrfA,  -  At] 

where: 

A,  =  cross  sectional  area  of  a  sphere 

Ak  =  total  cross  sectional  area  [spheres  and  corridor] 

For  the  space  station  the  cross  sectional  area  of  a  sphere  is: 

A,  =  7r[33.3  m]2  =  3.48  x  103  m2 
The  area  of  the  corridor  connecting  the  spheres  is: 

Ac  =  2[2[236m]  -  2[33.3m]][2m]  =  1.62  x  103  m2 
Therefore,  the  total  cross  sectional  area  becomes: 

At  =  2[3.48  x  103  m2]  -  1.62  x  103  m2  =  8.58  x  103  m2 
And  the  average  cross  sectional  area  is  then: 

A  =  [2/ 7T ][3.48  x  10s  m2  +  8.58  x  10s  m2]  =  7.68  x  103  m2 
The  velocity  of  the  space  station  is  dependent  upon  its  orbital  altitude 
and  can  be  calculated  from  the  formula  for  circular  orbit  speed  (19:75): 

v  =  (M/r]'n 

where. 

/i.  =  Earth’s  gravitational  constant  [3.986  x  IQ6  km3/sec2] 
r  =  radius  of  vehicle  from  Earth's  center 

Therefore: 

v  =  [[3.986  x  10B  km3/ sec2]  /  [6378km  +  240km]],/2  =  7.76  km/sec 
The  density  of  the  atmosphere  at  240  km  is  selected  from  NASA 
tables  of  density  for  standard  atmosphere  at  various  altitudes  and  is  1.24  x  10~‘° 
kg/m3.  (29:11-26) 

The  Coefficient  of  Drag  is  a  dimensionless  figure  which  carries  over 
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from  the  early  days  of  aerodynamics  but  remains  applicable  even  for  the 
conditions  in  which  spacecraft  operate.  It  is  dependent  upon  free-molecule  flow  at 
the  orbital  altitude  and  is  calculated  by  assuming  that  the  spacecraft  is  stationary 
and  the  air  molecules  are  flowing  past.  It  is  assumed  that  the  molecules  have  a 
Maxwellian  distribution  of  thermal  velocity  superimposed  on  their  uniform  velocity, 
v,  and  are  temporarily  retained  on  the  surface  after  collision,  then  re-emitted. 
(21:13-14) 

The  drag  coefficients  of  bodies  of  various  shapes,  at  various  angles  to 
the  air  flow,  and  the  average  value  for  rotating  bodies  have  been  evaluated 
experimentally.  For  satellites  with  perigee  heights  between  180  and  500  km  and 
orbital  eccentricities  between  0  and  0.2,  the  value  of  Co  is  between  2.1  and  2.2  for  a 
sphere,  and  approximately  2.15  for  a  cylinder  rotating  like  a  propeller  in  the  air 
flow.  Based  on  this  information  CD  can  be  approximated  by  a  value  of  2.2  with  a 
standard  deviation  which  should  not  exceed  5  percent.  (21:15) 

The  acceleration  due  to  drag,  ad,  is  now  calculated  by  substituting  the 
above  information  into  equation  [3]: 

ad  =  2.2  [1/2]  [1.24  x  lO'10  kg/m3]  [7.68  x  103  m2  /  1.45  x  106  kg]  [7.76  km/sec]2 
=  4.19  x  1CT6  m/sec2 

The  next  section  uses  this  acceleration  to  determine  the  force 
requirement  necessary  to  counteract  the  drag  and  how  much  of  the  force  the 
biogas  from  the  digester  system  can  produce  when  used  in  a  rocket  system. 

Orbit  Maintenance  Requirements 

The  force  of  drag  on  the  space  station  must  be  counteracted  by  an 
opposite  force  provided  by  the  propellant  used  in  the  orbit  maintenance  system. 
This  force  is  equal  to  the  mass  of  the  space  station  times  the  acceleration  due  to 
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atmospheric  drag  acting  on  it. 

F  =  mad  =  [1.45  x  106  kg][4.19  x  10’6  m/sec2] 

=  60.8  kg-m/sec2 

In  calculating  the  required  flow  rate  necessary  to  provide  this  force,  the 
use  of  the  digester  gas  as  a  cold  gas  propellant  is  considered  first.  Assuming  an 
approximate  70/30  mix  of  CH4/C02,  the  biogas  produced,  the  first  step  is  to 
calculate  the  exhaust  velocity  of  the  gas  as  it  exits  from  the  orbit  maintenance 
rocket  nozzle.  Exhaust  velocity  is  a  function  of  the  rocket  chamber 
pressure/exhaust  pressure  ratio  and  the  specific  heat  ratio,  and  is  proportional  to 
the  square  root  of  the  absolute  temperature,  at  the  rocket  nozzle  inlet,  and  the  gas 


constant.  (45:55-56] 


v.  =  {(2gk/(k-l]]RTc[l-(Pe/Pc]|k-1)/k]},/2 


where: 


k  =  specific  heat  ratio  of  the  gas 
R  =  gas  constant 

Tc  =  rocket  nozzle  inlet  temperature 
Pe  =  exit  plane  pressure 
Pc  =  chamber  pressure 
g  =  conversion  factor  [32.2  lbm— ft/lbf— sec2] 

Note,  since  gas  properties  are  tabulated  in  most  reference  books  in  English  units, 
the  following  calculations  will  be  made  in  English  units  and  then  converted  to 


Metric. 


The  pressure  ratio,  Pe/Pc,  is  a  function  of  the  ratio  of  the  nozzle  exit 


area  to  the  throat  area.  In  this  situation  the  nozzle  design  is  an  unknown,  however 
the  pressure  ratio  can  be  assumed  to  be  very  small  since  the  pressure  against  the 
exit  plane  of  the  nozzle  is  the  vacuum  of  space.  In  small  rocket  nozzle  design  the 
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nozzle  exhaust  area  to  throat  area  ratio  can  be  made  very  large,  also  providing  a 
very  small  pressure  ratio.  In  calculation  of  exhaust  velocity,  therefore,  it  is 
assumed  that  the  pressure  ratio  is  negligible,  and  the  value  of  the  quantity  within 
the  square  brackets,  [  j,  in  Eqn  [4]  goes  to  unity.  (8) 

The  gas  constant,  R,  depends  on  the  molecular  weight  of  the  CH4/C02 
mix.  For  the  70/30  mix  combination: 

Mmix  =  £x,mj  =  0.7[16  lbm/lbmole]  G.3[44  lbm/lbmole] 

=  24.4  lbm/lbmole 


where: 


x,  =  percent  of  constituent  by  volume 
m,  =  molecular  weight,  of  constituent 


therefore: 


where: 


Rmix  =  JJ/Mmix  =  [1545  Ft— lbf/lbmole— R]/[24.4  lbm/lbmole] 
=  63.28  ft-lbf/lbm-R 


i l  =  universal  gas  constant  (39:302) 

The  specific  heat  ratio  of  the  gas  is  the  ratio  of  the  specific  heat  at 
constant  pressure  and  the  specific  heat  at  constant  volume.  Each  of  these  specific 
heats  is  dependent  upon  the  CH4/C02  mix  of  the  gas.  The  molar  value  of  specific 
heat  is  first  calculated  as  the  sum  of  the  mole  fraction  and  the  molar  specific  heat 
products  of  each  of  the  constituents  as  was  done  for  molecular  weight.  (8;  46:550) 
CpmiI  =  Ex,Cp,  =  0.7(8.533  Btu/lbmole-R]  0.3(8.934  Btu/lbmole-R] 

=  8,6533  Btu/lbmole-R 


Cvmix  =  Sx/Hv,  =  0.7(6.464  Btu/lbmole-R]  0.3(6.9536  Btu/lbmole-R] 
=  6.6109  Btu/lbmole-R 


Cp  =  specific  heat  at  constant  pressure 

Gv  =  specific  heat  at  constant  volume 
The  ratio  of  specific  heat,  k,  for  the  gas  mixture  is  then: 

k  =  Cp/Cv  =  8.6533/6.6109  =  1.309 

Assuming  a  cold  gas  temperature  of  530R  [ambient  temperature,  approx. 
70F]  and  performing  the  calculations  in  Eqn  [4],  the  velocity  of  the  exhaust  gas,  v„, 
equals  3025  ft/sec  [922  m/sec]. 

The  next  step  is  to  determine  the  mass  flow  rate  of  the  propellant  gas, 
m,  required  to  oppose  the  atmospheric  drag  force,  and  compare  this  value  with  the 
mass  flow  rate  of  the  generated  gas.  Using  m  =  F/v„  where  F  is  the  drag  force, 
a  value  of  mass  flow  rate  of  5898  kg/day  is  calculated.  The  available  mass  flow 
rate  of  biogas  is  a  product  of  the  volume  of  the  gas  generated  times  it’s  density. 
The  density  of  biogas  is  approximately  1.06  kg/m3  [41:123) 

=  pV  =  [1.06  kg/ms][18.2  m3/day]  =  19.26  kg/day 

Clearly,  the  biogas  produced  if  used  in  a  cold  gas  orbit  maintenance 
system  will  not  meet  even  1%  of  the  propulsion  requirement,  and  therefore  is 
considered  no  further. 

Additional  energy  can  be  extracted  from  the  gas  if  burnt  with  an 
oxidizer.  An  oxidizer  commonly  used  in  space  vehicles  today  is  oxygen  which  is 
well  suited  as  a  combustion  reactant  with  methane.  The  only  portion  of  the  biogas 
generated  that  will  react  with  the  oxygen  is  the  70%  methane.  The  calculations 
used  in  determining  the  exhaust  velocity  of  the  hot  biogas/oxygen  combustion 
mixture  are  similar  to  those  made  for  the  cold  gas  system.  The  rocket  nozzle  inlet 
temperature  in  this  case  will  depend  on  the  heat  of  reaction  of  the  methane  with 
the  oxygen: 


CH4  -*•  202  ->  C02  -  2H20 

The  heat  of  reaction  will  also  determine  if  dissociation  of  the  gases  will  take  place. 
Dissociation  is  the  breaking  of  the  gasses  into  various  elements  and  molecules, 
removing  heat  energy  from  the  reaction  and  lowering  the  overall  temperature  of  the 
reaction.  (8) 

To  find  the  temperature  of  the  reaction,  thermochemical  tables  are  used 
to  balance  the  enthalpy  of  the  reaction  with  the  enthalpy  change  necessary  to 
bring  the  products  to  adiabatic  flame  temperature  in  a  trial  and  error  fashion  until 
a  temperature  can  be  located  at  which  these  quantities  are  equal.  The  trial  and 
error  calculations  are  presented  in  Appendix  F.  The  equations  used  in  the 
calculations  are  as  follows  [the  subscripts  p  and  r  refer  to  product  and  reactant 
respectively]: 

AHt  =  2npAHfp  -  En  AHfr  =  Qoul 
AHr  =  Snp/Cp  dT  =  Enp[H-H298]  =  Q,n 

where: 

AHr  =  heat  of  reaction 

n  =  number  of  moles 

Hf  =  heat  of  formation  [kcal/mole] 

H-H29#  =  change  in  enthalpy  in  going  from  reference  temperature 
[298R]  to  adiabatic  flame  temperature 
[from  thermochemical  tables] 

Q  =  energy  [kcal] 

Various  adiabatic  flame  temperatures  are  chosen  until  Qoul  =  Qin.  In  this  case, 
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when  dissociation  is  not  considered,  the  resulting  adiabatic  flame  temperature,  Ta, 
would  be  approximately  4600K  [See  Appendix  F],  At  this  temperature 
dissociation  will  take  place  and  some  energy  will  be  absorbed  in  the  process. 
In-depth  thermochemical  analysis  is  necessary  to  determine  the  exact  composition 
and  adiabatic  flame  temperature  that  will  be  produced  in  such  a  dissociation 
reaction.  The  calculations  involve  a  trial-and-error  simultaneous  solution  of  mass 
balance  and  equilibrium  equations,  which  is  a  problem  of  considerable  complexity. 
For  this  reason  a  computer  program  is  utilized  in  determining  the  adiabatic  flame 
temperature.  An  explanation  of  the  procedure,  the  computer  program  listing  and 
output  are  also  presented  in  Appendix  F.  A  temperature  of  approximately  3380K 
is  determined  using  the  computer  program,  assuming  a  gas  pressure  of  100  psia. 

The  product  gas  used  in  calculating  the  adiabatic  flame  temperature  is 
then  used  to  calculate  the  ratio  of  specific  heats.  The  constituents  of  the  product 
gas  are  H20,  OH,  02,  0,  H2,  and  C02.  The  additional  C02  present  in  the  biogas 
was  considered  in  the  thermochemical  analysis,  as  shown  below.  The  calculations 
are  made  on  a  molar  basis.  The  number  of  moles  of  methane  produced  per  day  is 
equal  to  the  percentage  of  methane  in  the  biogas  times  the  daily  mass  flow  rate  of 
the  biogas  divided  by  the  molecular  weight  of  the  biogas: 

,7[19.26  kg/day]/[24.4  kg/kgmole]  =  .5525  moles/day  [CH4] 

The  number  of  moles  of  C02  present  which  acts  as  a  diluent  is: 

,3[19.26  kg/day]/[24.4  kg/kgmole]  =  .1386  moles/day  [C02] 

The  left  hand  side  of  the  reaction  of  methane  with  oxygen,  including  the  carbon 
dioxide  as  a  diluent  is: 

,5525[CH4  -  202]  hk  ,2368C02 

Using  the  number  of  moles  of  each  constituent  produced  [Appendix  F]  the 
complete  reaction  can  be  written  on  a  per  mole  CH4  basis  as: 


CH4  -  202  0.4286CO2  ->  1.3322Ha  +  0.3173OH  -  0.09680 

-  0.1273O2  *  0.2134H  -  0.4025Ha 

-  1.4286COa 

The  right  hand  side  is  now  used  to  determine  the  molecular  weight  of  the  product 
gas.  First  the  mole  fraction  of  each  constituent  is  determined  using  the  calculated 
total  number  of  product  moles,  3.918,  listed  in  the  program  output  (Appendix  F). 
Mole  Fraction  [x,]: 

H2  :  0.3400 
OH:  0.0810 
O  :  0.0247 
02  ;  0.0325 
H  :  0.0545 
H2  :  0.1027 
C02  0.3646 

The  molecular  weight  of  the  mix  is  then  calculated  as: 

Mmix  =  Sx.m,  =  .3400(18]  -  .0810(17]  +  .0247(16] 

-  .0325(32]  -  .0545(1]  -  .1027(2]  -  .3646(44] 

=  25.23  Ibm/lbmole 

where: 

x,  =  mole  fraction  of  constituent  i 

m,  =  molecular  weight  of  constituent  i  (lbm/lbmole] 

The  gas  constant  of  the  mix  is  next  calculated: 

Rmix  =  Jt/Mmix  =  [1545  ft-lbf/lbmole-R]/[25.231bm/lbmole] 

=  61.24  ft-lbf/lbm-R 


As  the  gas  is  expanded  through  the  nozzle  it  experiences  a  drop  in  temperature 
and  pressure.  To  determine  the  temperature  at  the  nozzle  exit,  the  expansion 
process  is  considered  to  take  place  such  that  the  product  gas  composition  at  the 
nozzle  exit  is  the  same  as  in  the  combustion  chamber.  This  condition  is  known  as 
frozen  flow  or  frozen  equilibrium.  The  temperature  of  the  product  gas  at  the 
nozzle  exit,  Te,  is  then  calculated  from  the  formula:  (45:183;  8) 

Te/Tc  =  [Pe/Pcf ',1/k 

where: 

Pe/Pc  =  exit/chamber  pressure  ratio 

Tc  =  chamber  temp,  [adiabatic  flame  temperature] 

Te  =  nozzle  exit  temperature 
k  =  ratio  of  specific  heats 

In  making  this  calculation  an  estimate  for  the  ratio  of  specific  heats,  k,  and  the 
pressure  ratio,  Pe/Pc,  are  made.  As  an  initial  guess,  let  k  =  1.3.  The  pressure 
ratio  of  chamber  pressure,  Pc,  to  nozzle  exit  pressure,  Pe,  is  assumed  to  be  large, 
since  the  rocket  operates  in  the  vacuum  of  space,  and  a  value  of  1000:1  is  selected. 
The  nozzle  exit  temperature  is  then  estimated  to  be: 

Te  =  Tc[Pe/Pcf-11/k  [5] 

=  3380K  [.001],13-||/,s 
=  686K 

In  reality  the  exit  temperature  will  actually  be  somewhere  in  between 
the  adiabatic  flame  temperature  and  the  exit  temperature  calculated  in  frozen 
equilibrium ;  however,  the  temperature  difference  can  be  used  to  calculate  an 
average  specific  heat  for  each  product  constituent  using  the  formula; 


Gp(avg)  =  AH/AT 


where: 


=  change  in  enthalpy  in  going  from  chamber  temperature  to  exit  temperature 

=  change  in  temperature  between  combustion  chamber  and  exit  nozzle 

A  sample  calculation  is  made  here  for  H20  using  data  from  the 
thermochemical  tables  in  reference  (4],  For  ease  of  extracting  data  from  the  tables 
the  adiabatic  flame  temperature  is  rounded  to  3400K  and  the  exit  temperature  to 
700K. 

Cp(avg)  =  [(35.577  -  3.390)kcal/mole]/[(3400  -  700)K]  =  11.92  cal/mole-K 
It  can  be  shown  that  an  increase  in  the  pressure  ratio  has  very  little  effect  on  the 
resulting  Cp(avg)  value  and  is  not  critical  to  the  analysis.  For  example;  if  the 
pressure  ratio  were  2000:1,  the  resulting  exit  temperature  would  be  approximately 
600K  and  the  resulting  Cp(avg)  value  would  be  11.8.  Therefore,  the  remaining 
Cp(avg)  values  are  calculated  for  a  1000:1  pressure  ratio,  a  resulting  exit 
temperature  of  700K,  and  an  approximate  combustion  temperature  of  3400K. 
Average  specific  heats  [Cp(avg)]  in  cal/mole-K: 

H2  :  11.92 
OH:  8.18 
O  :  4.99 
02  :  9.01 
H  :  4.97 
H2  :  8.14 
C02  14.17 

The  product  gas  specific  heats  are  next  calculated  using  the  above  Cp(avg)  values 
and  the  mole  fractions,  x„  calculated  earlier. 
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Cp  =  Sx.Cp^avg)  =  ,3400[  11.92]  -  ,0810[8.18]  ^  ,0247[4.99] 

-  ,0325[9.01]  +  ,0545[4.97]  ^  ,1027[8.14] 

-  ,3646[14.17] 

=  11.4047  cal/mole-K 

and: 

Cv  =  Cp  -  * 

where: 

S.  =  1.986  cal/mole-K 

therefore: 

Cv  =  11.4074  -  1.986  =  9.4187  cal/mole-K 
The  ratio  of  specific  heats  is  then: 

k  =  Cp/Cv  =  11.4047/9.4187  =  1.211 

When  this  new  value  of  k  is  substituted  into  equation  [5]  the  exit  temperature 
becomes  1000K.  Repeating  the  above  procedures  with  this  new  exit  temperature 
results  in  a  ratio  of  spt  :ific  heats  of  1.206.  An  average  value  of  k  =  1.208  is  now 
used  in  Equation  [4]  with  T  =  3380K  =  6084R,  and  R  =  61.24  ft-lbf/lbm-R 
[assuming  Pe/Pc  negligible]: 

v,  =  {[2[1.208]/.208][61.24][6084][32.2]},/2 
=  11,852ft/ sec 
=  3611  m/sec 

The  required  mass  flow  rate  is  then: 

rii  =  F/v.  =  [60.8  kg-m/secJ]/[3611  m/sec] 

=  16,84  gm/sec 
=  1454.8  kg/day 

This  is  the  mass  flow  rate  of  the  product  gas  of  which  41.7%  [mole 
fraction]  was  initially  biogas.  Therefore  the  required  mass  flow  rate  of  the  biogas 


is  606.6  kg/day.  A  comparison  of  this  figure  with  the  mass  flow  rate  available 
from  the  digester  [19.26  kg/day]  shows  that  the  biogas  will  produce  3.2%  of  the 
force  required  to  overcome  atmospheric  drag  if  reacted  with  approximately  60% 
oxygen  [mole  fraction].  An  optimization  of  this  percentage  may  be  possible  by 
decreasing  the  amount  of  oxygen  in  the  reaction,  thus  providing  a  fuel  rich  mixture 
which  could  lower  the  molecular  weight  of  the  combustion  product.  Considering  the 
uncertainty  of  the  actual  composition  of  the  biogas,  the  3.2%  figure  is  sufficient. 

Cost  Analysis  -  Methane  ys  Shuttled  Propellants 

The  remainder  of  this  chapter  deals  with  estimating  the  cost  of 
providing  the  digester  produced  propellant  and  comparing  that  cost  with  the  cost  of 
supplying  conventional  propellants  to  the  space  station  via  the  space  shuttle.  The 
acquisition  and  operational  costs  are  the  only  costs  considered  in  this  analysis.  It 
is  assumed  that  a  waste  stabilization  system  of  some  type  [biological, 
physiochemically,  electrodialysis,  reverse  osmosis,  etc.]  would  be  necessary  on  a 
space  station  of  this  size,  and  that  the  development  costs  for  all  such  systems 
would  be  similar.  Development  costs  are  therefore  not  considered  to  be  a  key 
factor  in  this  study.  Acquisition  costs,  however,  are  equipment  dependent  and 
would  reflect  the  specific  system  used  to  process  the  waste  materials,  in  this  case 
anaerobic  digestion.  The  cost  of  other  systems  which  might  be  used  to  provide 
waste  stabilization  have  not  been  balanced  against  the  cost  of  the  anaerobic 
digester  because  of  the  complexity  involved  in  evaluating  and  costing  the  multitude 
of  possible  systems.  (13:66-67;  26) 

There  are  many  other  elements  which  may  have  a  bearing  on  the  value 
of  a  methane  generating  system  that  have  not  been  considered  in  the  cost  analysis 
because  of  the  difficulty  in  placing  a  number  on  their  value.  These  include  the 


value  of  the  effluent  returned  from  the  system  after  digestion,  the  value  in 
reducing  the  mass  of  waste  products,  and  the  utility  value  of  stabilized  waste 
material  in  space,  to  mention  a  few.  Since  the  value  of  these  elements  are 
unknown,  this  study  can  only  be  considered  to  be  a  partial  analysis;  however,  these 
factors  must  be  included  in  any  decision  to  use  such  a  system. 

The  Monte  Carlo  cost  estimating  technique  is  used  in  this  study  to 
estimate  the  acquisition  and  operational  costs  of  the  digester  system.  It  is  also 
used  to  provide  a  cost  estimate  of  the  conventional  propellant  requirement 
necessary  to  provide  the  additional  3.2%  thrust  energy  equivalent  of  the  digester 
produced  methane  propellant  for  the  cost  comparison.  The  balance  of  the  thrust 
requirement,  the  additional  96.8%,  must  be  provided  by  conventional  propellants  in 
either  case,  and  is  therefore  not  considered  as  a  cost  element  in  this  study. 

Monte  Carlo  Cost  Estimating  Technique.  The  anaerobic  digester 
proposed  in  this  study  is  based  on  principles  developed  on  Earth  for  large  scale 
municipal  waste  treatment.  Scaling  down  a  system  of  this  type  to  a  size  capable 
of  handling  200  people  and  requiring  that  the  system  operate  in  space  places  a 
good  deal  of  uncertainty  on  the  estimates  of  the  resource  requirements.  There  are 
two  categories  of  uncertainty  involved  in  estimating  the  resource  requirements: 
requirements  uncertainty  and  cost-estimating  uncertainty.  (13:205-207) 

Requirements  uncertainty  refers  to  variations  in  the  cost  estimate  due 
to  various  possible  configurations  of  the  system  under  study.  Possible  causes  of 
requirements  uncertainty  may  stem  from  assumptions  regarding  hardware 
characteristics  or  system  operational  concepts.  Studies  have  shown  that 
requirements  uncertainty  accounts  for  70  to  80  percent  of  the  total  cost  estimate. 
Uncertainties  which  may  affect  requirements  for  the  anaerobic  digester  proposed 


here  include  space  applicable  materials  selection,  for  pumps  and  other  components 
of  the  digester,  configuration  of  the  system  design,  and  estimation  of  flow  rates, 
system  capacities  and  operational  and  maintenance  requirements.  (6:1) 

Cost-estimating  uncertainty  refers  to  the  variation  in  cost  estimates 
due  to  the  unpredictable  nature  of  estimating  relationships.  In  other  words,  there 
is  usually  some  randomness  in  the  variables  that  are  chosen  to  predict  a  particular 
cost  that  must  be  accounted  for  in  the  cost  estimating  relationship.  Sometimes 
this  uncertainty  can  arise  from  the  data  that  is  used  as  a  basis  for  the  cost 
analysis.  The  data  may  be  statistical,  and  therefore  only  as  good  as  the 
observations  from  which  it  was  based,  in  which  case  there  is  some  random 
deviation  of  this  data  from  the  "true"  value.  The  data  may  be  based  on  past  or 
current  experience,  in  which  case  it  must  be  realized  that  relationships  which  held 
in  the  past  may  not  hold  for  future  advanced  systems.  In  estimating  costs  of 
future  systems,  necessity  may  dictate  the  extrapolation  of  relationships  beyond  the 
data  base  causing  expansion  of  prediction  intervals  and  inducing  uncertainty  into 
the  estimates.  (13:206-207;  27:IV-2  -  IV-3) 

The  intent  of  the  Monte  Carlo  method  of  cost  estimating  is  to  describe 
the  above  uncertainties  for  each  of  the  cost  elements  of  the  system  so  that  the 
extent  of  the  uncertainty  of  the  ultimate  system  cost  can  be  anticipated  and 
evaluated.  In  this  way  alternatives  can  be  compared  with  respect  to  their  expected 
cost  and  their  uncertainties.  (6:2,4) 

The  method  requires  the  pricing  expert  to  specify  a  highest  and  lowest 
expected  value  for  each  cost  element  of  the  system  in  addition  to  a  selection  of 
one  of  nine  beta  distributions  which  best  represent  the  cost  uncertainty  of  each 
cost  element  based  on  his  experience.  The  beta  distributions  selected  have  many 
of  the  characteristics  that  would  be  expected  in  input  uncertainties:  upper  and 


lower  bounds,  continuity,  and  unimodal  distribution.  These  distributions  and  their 
associated  CX  and  /?  values  are  depicted  in  Figure  10.  From  each  of  these 
distributions  a  cumulative  distribution  is  plotted  [Figure  11],  A  random  number 
between  0  and  1  is  then  generated  and  located  on  the  vertical  axis  of  the 
cumulative  distribution  plot.  The  corresponding  x  value  from  the  horizontal  axis  is 
then  placed  in  the  formula 

C,  -  CL  [Cl  CH]x 

where: 

C,  =  the  computed  cost  value  for  cost  element  / 

Cl  =  the  lowest  expected  value 
Ch  =  the  highest  expected  value 

to  determine  a  probable  value  of  the  cost  element.  This  procedure  is  repeated  for 
each  cost  element  in  the  system  and  the  computed  value  of  each  cost  element  is 
summed  to  determine  a  probable  total  cost  of  the  system.  The  probable  total  cost 
estimate  of  the  system  is  then  computed  over  several  iterations  to  produce  a  total 
system  cost  distribution.  The  use  of  this  method  requires  that  all  the  input 
parameters  be  mutually  independent  which  is  normally  a  valid  assumption  with 
cost  factor  inputs.  (5;  7:15;  43:154-155) 


Figure  11.  Input  Uncertainty  Distributions  (43:155) 

In  order  to  develop  a  system  life  cycle  cost  using  the  Monte  Carlo 
technique  a  discount  rate  must  be  considered  so  that  future  dollars  can  be 
commensurate  with  todays  dollars.  The  formula  for  computing  the  life  cycle  cost 
is  then: 

LCC  =  S'.fn,  C,/[r  -  1]'] 

where: 

LCC  =  life  cycle  cost 

y  =  number  of  years  from  the  present  that  cost  element  is  purchase 
n,  =  number  of  units  of  cost  element  i  purchased 
r  =  discount  rate 

The  life  cycle  cost  is  computed  over  all  applicable  major  cost  categories  of  the 
system  [Development,  Acquisition,  and  Operations]  to  develop  a  total  cost 
estimate  uncertainty.  (7:16-18;  13:66-69) 


Cost  Catagories.  As  was  mentioned  earlier,  only  two  cost  categories 
are  considered  relevant  to  this  study:  Acquisition  and  Operations.  In  the 


Acquisition  category  of  the  digester  system,  the  initial  investment  costs  include 
the  individual  hardware  components  [pumps,  tanks,  plumbing,  etc.],  shuttle  delivery 
costs,  and  labor  to  assemble  the  digester  as  cost  elements.  Shuttle  delivery  is  a 
variable  which  is  dependent  on  the  weight  of  each  individual  component.  Since  it 
is  unknown  exactly  what  each  component  will  weigh,  the  Monte  Carlo  technique 
was  also  applied  to  capture  the  weight  uncertainty  in  the  life  cycle  cost  estimate. 
A  detailed  cost  breakdown  of  the  cost  elements  under  the  Acquisition  cost 
category  is  presented  in  Table  III  along  with  the  selected  beta  uncertainty  curve 
type  and  other  information  pertinent  to  the  analysis.  All  acquisitions  take  place 
during  the  first  year  of  the  Acquisitions  phase  with  the  exception  of  construction 
and  transport  to  space  which  takes  place  in  the  second  and  final  year  of  the 
Acquisitions  phase.  Table  IV  contains  the  weight  values  of  each  cost  element 
hardware  component  with  the  corresponding  beta  distribution  type  describing  the 
uncertainty  The  hardware  cost  elements  were  based  on  a  system  configuration 
[Figure  12]  developed  with  the  assistance  of  waste  management  engineers  to  be 
compatible  with  the  two  stage  high  rate  digester  concept  initially  proposed.  Cost 
element  and  weight  element  values  were  also  acquired  from  these  experts  through 
personal  interviews  and  correspondence.  (17;  24;  47) 

The  system,  as  configured  in  Figure  12,  is  capable  of  handling  the  waste 
input  of  the  entire  space  station.  Redundancy  is  provided  by  locating  a  system  of 
this  type  in  each  end  of  the  space  station.  A  waste  transfer  pipe  running  the 
length  of  the  corridor  will  transfer  waste  from  one  unit  of  the  space  station  to  the 
other  in  the  event  of  a  system  failure.  Table  III  reflects  the  redundancy  in  the 
equipment  and  additional  spare  components. 


TABLE  III 

Cost  Elements  for  Aquisition  Phase 


Cost[’85$]  Beta 


Cost  Element 

Low 

High  Curve  Quanity 

Remarks 

Tank 

13,000 

20,000 

4 

4 

Includes  delivery 

Tank  Insulation  6,000 

7,600 

2 

2 

3"  Fiberglass  with  service  jacket 

Pump  1 

800 

1,200 

8 

3 

Includes  1  spare 

Pump  2 

2200 

2,600 

8 

3 

Includes  1  spare 

Pump  3 

800 

1,200 

8 

3 

Includes  1  spare 

Pump  4 

800 

1,200 

8 

3 

Includes  1  spare 

Compressor 

800 

1,200 

8 

3 

Includes  1  spare 

Heat  Exchanger  610 

1,000 

1 

3 

Includes  1  spare 

Temp  Control 
Valve 

310 

350 

5 

3 

Includes  1  spare 

Pipe 

15,587 

17,100 

7 

1 

2100'  of  4"  Pipe  and  475’  of  1.5" 

Valve 

270 

300 

9 

16 

4"  Valves,  2  spares 

Coupling 

13 

15 

9 

105 

4"Couplings,  5  spares 

Flame  Trap 

800 

900 

9 

5 

Includes  1  spare 

Pressure  Relief  590 
Valve  <8  Vacuum 
Breaker  w/Flame 

Trap 

700 

9 

5 

Includes  1  spare 

Sediment  <4  Drip  1,050 
Trap  Assembly 

1,150 

9 

2 

Oxygen  Tank 

6,500 

10,000 

1 

1 

2300  ft3  tank  <8  2200psi 

Construction 

30,200 

50,400 

4 

1 

Based  on  30%  total  cost  (47), 
includes  space/grnd  uncertainty 

Shuttle  Transport  656  2,647  6  vt. 


$/lb.  Based  on  STS  or  derived 
vehicle  (25) 


TABLE  IV 
Weight  Elements 

Weight  Beta 


Weight  Element  Low 

High  Curve 

£uanlkZ 

Remarks 

Tank  13,000 

18,000 

4 

4 

Tank  Insulation 

760 

800 

6 

2 

3"  Fiberglass  with 

service  jacket 

Pump  1 

50 

75 

9 

3 

Includes  1  spare 

Pump  2 

180 

200 

8 

3 

Includes  1  spare 

Pump  3 

50 

75 

9 

3 

Includes  1  spare 

Pump  4 

50 

75 

9 

3 

Includes  1  spare 

Compressor 

50 

75 

7 

3 

Includes  1  spare 

Heat  Exchanger 

300 

500 

4 

3 

Includes  1  spare 

Temp  Control 
Valve 

10 

20 

9 

3 

Includes  1  spare 

Pipe  20,000 

23,000 

6 

1 

Includes  Couplings 

Valve 

35 

40 

4 

16 

Includes  2  spares 

Flame  Trap 

40 

50 

8 

5 

Includes  1  spare 

Pressure  Relief 
Valve  i,  Vacuum 
Breaker  w/Flame 
Trap 

40 

50 

8 

5 

Includes  1  spare 

Sediment  i.  Drip 
Trap  Assembly 

140 

155 

8 

2 

Oxygen  Tank 

n/ 

a  39,357 

1 

1 

lbs  oxygen/year  •+■ 

11,000  lb  tank 

63 


r. 


i 


i£S 


Key  to  Major  Items: 


|  PI  -  Will  be  used  to  pump  undigested  sludge  from  one  end  of  the  space  ' 

station  to  the  other  digester  if  a  digester  requires  maintenance.  It  has  a  j 

capacity  to  handle  approximately  140  gals/day  [.1  gpm]  and  can  pump  | 
straight  up  to  a  height  of  780  feet.  At  this  height  the  center  of  rotation  of  j 
the  space  station  is  encountered  and  there  is  zero  gravity.  Therefore  the  I 
average  height  that  the  pump  "sees"  it  must  pump  against  is  about  400  ft. 
At  the  center  of  the  space  station  the  sludge  begins  to  travel  "down-hill" 
toward  the  other  half  of  the  space  station,  assisted  by  artificial  gravity. 

j 

!  P2  -  Will  serve  to  carry  the  digested  sludge  through  the  heat  exchanger  and 

1  also  to  mix  the  contents  of  the  tank  for  proper  feeding  of  the  organisms. 

!  The  tank  is  complete  mixed  once  per  day.  Pump  capacity  is  about  13  gpm. 

|  P3  -  Pumps  contents  of  mix  tank  [Tl]  to  settling  tank  [T2]  at  the  same  rate 
|  as  raw  sludge  is  input  into  Tl. 

'  P4  -  Will  pump  concentrated  settled  sludge  from  T2  to  a  sludge  drying 
j  system  of  some  type.  The  SLAM  model  predicts  an  output  from  T2  of 
about  107  Ibs/day  of  digested  sludge. 

Compressor  -  Will  compress  "biogas",  generated  at  a  rate  of  .45  cfm  [646 
cf/day]  to  100  psia. 

H  E.  -  Heat  exchanger.  Water  type  heat  exchanger  heated  by  electric  coils. 


Figure  12.  System  Configuration 


The  Operations  cost  category  includes  the  purchase  and  shuttle  delivery 
costs  of  oxygen  propellant,  to  be  used  as  an  oxidizer  with  the  digester  generated 
fuel  propellant,  and  the  operations  and  maintenance  costs  of  the  system.  To 
determine  the  amount  of  oxygen  required  per  year,  the  mole  fractions  of  biogas 
and  oxygen  calculated  in  the  last  section  were  used.  For  the  described  reaction 
these  were,  for  oxygen  58.3%,  and  for  biogas  41.7%.  The  mass  flow  rate  of  biogas 
produced  was  calculated  to  be  19.26  kg/day  having  a  molar  weight  of  24.4 
kg/kgmole.  Therefore  the  biogas  molar  portion  of  the  reactants  is: 

[19.26  kg/day]/[24.4  kg/kgmole]  =  0.789  kgmole/day 
and  the  oxygen  molar  portion  of  the  reactants  is: 

[.789  kgmole/day][.583]/[.417]  =  1.104  kgmole/day 
The  molar  weight  of  oxygen  is  32  kg/kgmole,  therefore  the  weight  of  oxygen 
required  for  the  reaction  per  day  is: 

[1.104  kgmole/day][32kg/kgmo!e]  =  35.31  kg/day 
or  77.69  Ibm/day  at  STP  which  represents  a  yearly  requirement  of  28,357  lbs  of 
oxygen.  Operations  and  maintenance  costs  include  labor  and  materials  and  were 
derived  from  various  EPA  charts.  Table  V  lists  the  values  of  each  cost  element 
in  the  operations  phase  with  its  corresponding  beta  distribution  type  describing  the 
uncertainty.  Costs  which  are  based  on  historical  data  were  adjusted  to  current 
year  dollars  by  an  appropriate  inflation  factor  of  9%  which  is  applicable  to  Defense 
Department  studies.  (10:Charts  D-l  -  D— 14;  11:63-71,  Figures  A-17  -  A-40; 


TABLE  V 

Cost  Elements  for  Operational  Phase 
Cost[’85$]  Beta 


Cost  Element  Low 

High  Curve 

Remarks 

Labor 

1.0  M 

1.4  M 

3 

4  Maint.  Men,  1  Lab  Tech.  $120,000/yr/man 
[75$]  (47;  34:160) 

Maintenance 

Supplies 

6,100 

7,100 

5 

Pump  belts,  lubrication,  etc.  (11) 

Energy 

10,200 

16,600 

1 

Based  on  10%  of  total  cost  (47) 

Oxygen 

Purchase 

17,200 

17,400 

8 

Based  on  $.05  per  cubic  foot  at  STP 

Oxygen 

Transport 

25.8  M  104  M 

6 

39,357  lbs  [Oxygen  ■+■  Tank]  mult  by  Shuttle 
cost  per  lb.  (25) 

The  alternative  to  providing  a  portion  of  the  orbit  maintenance 
requirement  with  the  digester  produced  propellant  is  to  use  conventional  fuels  to 
provide  that  portion.  The  life  cycle  cost  of  this  alternative  is  computed  using  the 
same  Monte  Carlo  method  as  before  but  considering  only  the  Operations  cost 
category  since  conventional  propellants  will  be  purchased,  rather  than  produced, 
and  delivered  to  the  space  station  on  a  regular  basis.  The  cost  elements  involved 
with  this  alternative  are  purchasing  and  shipment  of  the  propellants  and  are  listed 
in  Table  VI.  Prcpellant  purchase  costs  are  standard  stock  fund  prices  for  DOD 
and  NASA.  For  MMH  the  price  is  $8  -  $10  per  pound  and  for  for  N204,  $2.50  - 
$2.75  per  pound.  The  variation  in  price  reflects  the  fluctuation  observed  during  the 
1985  calander  year.  (25  4-10;  28) 


TABLE  VI 

Cost  Elements  for  Conventional  Propellants  [Operational  Phase] 

Cost[85$]  Beta 

Cost  Element  Low  High  Curve  Remarks 

Propellant  Tank  3,300  5,100  1  One  Time  Purchase,  Resusable 

N204  67,963  74,758  4  Based  on  27,185  Ibm  (28) 

MMH  131,808  164,760  4  Based  on  16,476  lbm  (28) 

Propellant  31.2M  126.0M  6  Based  on  cost/lb  for  43,661  lbs 

Transport  propellant  and  a  4,000  lb  tank 

The  life  cycle  cost  is  dependent  upon  the  span  of  time  over  which  the 
system  is  procured  and  operated,  and  the  time  frame  of  acquisition  and  operation 
to  which  the  discount  rate  will  be  applied.  In  order  to  calculate  a  life  cycle  cost 
several  assumptions  were  made  in  scheduling  the  acquisition  and  operational 
phases.  After  the  development  stage  is  completed  acquisition  and  on-ground 
preassembly  could  take  place  in  one  year  and  packaging  and  shipping  the  next 
year.  The  operational  phase  would  last  a  lifetime  of  30  years  beginning  in  2005 
and  therefore  the  acquisition  phase  would  begin  in  2003.  Actual  start  dates  for 
these  phases  are  not  important  in  the  analysis,  but  they  do  allow  for  a  common 
reference  point  when  applying  a  discount  rate.  An  average  discount  rate  of  10%, 
appropriate  to  Defense  Department  studies,  is  used  in  this  project.  (13:227) 

Equivalent  Propellant  Determination.  A  comparison  of  the  methane 
generated  propellant  to  shuttled  conventional  propellant  requires  that  the  amount 
and  cost  of  the  conventional  propellant  be  determined.  To  provide  a  convenient 
basis  for  comparison  the  Shuttle  Orbital  Maneuvering  Subsystem  [OMS]  was 
selected  as  a  typical  propulsion  system  that  could  be  employed  on  the  space 


station  to  provide  orbit  maintenance  functions.  This  system  uses  Nitrogen 
Tetroxide  [N204]  as  an  oxidizer  and  Monomethalhydrazine  [MMH]  as  a  fuel  The 
question  is  then,  how  much  of  these  propellants  are  required  to  produce  ‘he 
additional  3.2%  of  the  force  required  to  overcome  atmospheric  drag  if  the  biogas 
system  were  not  used?  (1:13-7) 

The  force  required  to  counteract  the  atmospheric  drag  was  calculated  in 
the  section  before  last  [page  46]  as  being  60.8  kg-m/sec2.  Three  and  two-tenths 
percent  of  this  figure,  or  1.93  kg/sec2  is  then  the  required  force  that  the 
conventional  propellant  system  must  provide.  The  Shuttle  OMS  system  can 
provide  a  specific  impulse  of  313  seconds.  Specific  impulse  is  a  rocket  performance 
parameter  which  relates  thrust  to  the  propellant  mass  flow  rate.  Using  the  formula 
for  specific  impulse  the  mass  flow  rate  required  to  produce  a  force  of  1.93 
kg-m/sec2  can  be  calculated: 

m  =  F/[Igp  ge]  =  1.93  kg-m/sec2  /  313  sec  x  9.8  m/sec2 
=  6.29  x  10""4  kg/sec  =  1.38  x  10“s  lbm/sec 

where: 

gc  =  gravitational  constant 

This  is  equivalent  to  a  propellant  flow  requirement  of  43,661  lbm/yr.  The  mass 
ratio  of  N204  to  MMH  for  this  system  is  1.65  to  1,  therefore  this  flow  requirement 
equates  to  16,476  lbm/yr  of  MMH  and  27,185  lbm/yr  of  N204.  Since  the  shuttle 
can  carry  a  maximum  of  65,000  lbs  to  an  altitude  of  240  km,  this  will  not  pose  a 
weight  restriction  on  how  much  propellant  can  be  delivered  on  a  flight.  To 
determine  if  there  is  a  volume  restraint  on  the  volume  of  propellants  which  can  be 
carried  into  space,  the  densities  [at  STP]  of  N204  [89.1  lbm/ft3]  and  MMH  [48.6 
lbm/ft3]  are  used  to  calculate  the  volume  of  the  above  mass  of  propellants: 


VN204  =  m/p  -  27,185  !bm  /  89.1  lbm/ft3  =  305.1  ft3 


and: 

Vmmh  =  m/p  =  16,476  ibm  /  48.6  lbm/ft3  =  339.0  ft3 
This  represents  a  total  shipping  volume  requirement  of  644.1  ft3.  The  cargo  bay  of 
the  Shuttle  has  a  volume  of  approximately  10,600  ft3,  and  therefore  the  volume 
requirement  of  the  propellants  is  not  a  restraint  and  shuttle  transport  costs  can  be 
calculated  on  solely  a  cost  per  pound  basis.  The  transport  cost  element  is 
reflected  in  Table  VI  together  with  the  purchase  cost.  (1:13-6,13-7;  45:29,  297) 


Analysis.  In  order  to  compute  the  life  cycle  cost  of  the  systems  a  set 
of  equations  is  used  to  determine  the  argument,  x,  of  the  cumulative  distribution 
function,  F(x),  of  the  nine  possible  beta  curves  in  Figure  10.  The  equations  were 
developed  from  tables  of  the  incomplete  beta  function  for  the  given  a  and  /? 
parameters  One  hundred  points  from  the  tables  were  placed  in  a  regression 
program  which  computed  the  polynomials  for  the  curves  that  best  fit  the  data 
points  The  equations  were  then  used  to  create  values  and  these  were  checked 
against  the  tables.  The  equations  are: 

Type  Equation 

1  x  =  4.4519  F(x)  -  11.425  F(x)2  -  14.898  F(x)3  -  7.3393  F(x)4 

2  x  =  1.5828  F(x)  -  1.72666  F(x)2  +  1.1339  F(x)3 

3  x  =  -.25015  F(x)  +2,265  F(x)2  -  3.2916  F(x)3  +  2.2205  F(x)4 

4  x  =  3.5761  (F(x)+.05)  -  9.6707  (F(x)2  +  12.954  F(x)3  -  6.2156  F(x)4 


5  x  =  2.238  F(x)  -  3.6096  F(x)  +  2.3389  F(x)3 

6  x  =  -.37482  F(x)  +  4.8703  F(x)2  -  9.1390  F(x)3  +  5.4841  F(x)4 

7  x  -  2.6547  (F(x)+.l)  -  5.7211  F(x)2  +  6.0833  F(x)3  -  2.3001  F(x)4 


Type 


Equation 


vr?  r*. 
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8  x  =  27511  F(x)  -  5.0861  F(x)2  -3.2677  F(xf 

9  x  —  .014397  F<x)  -  3  3829  F(x)2  -  7,1725  F(x)3  -4.5308  F(x)4 


Constraints:  0  <  F(x)  <  1 


0  <  x  <  1 


These  equations  do  not  fully  describe  the  argument  of  the  cumulative  distribution 
functions  but  are  adequate  except  at  the  extreme  tails  of  the  distribution.  Very 
rarely  the  value  of  x  may  be  computed  outside  of  the  constraint  range,  but  this 
will  have  no  impact  on  the  final  life  cycle  cost  distribution.  These  equations  were 
developed  into  a  FORTRAN  function  called  BETA  and  used  to  determine  random 
variables  for  weight  and  cost  using  the  related  beta  uncertainty  curve  for  the 
particular  cost  element  (5,  40) 

The  flow  chart  describing  the  FORTRAN  program  used  in  determining 
the  system  life  cycle  cost  is  shown  in  Figure  13.  The  range  of  possible  life  cycle 
cost  values  that  the  system  could  have  is  computed  from  the  sum  of  the  high  and 
low  values  for  each  of  the  cost  elements  discounted  back  to  1985.  This  range  is 
then  divided  into  100  intervals  for  later  use  in  testing  each  life  cycle  cost  generated 
to  see  which  interval  it  is  assigned  to.  The  program  then  determines  a  weight  for 
each  cost  element  based  on  the  high  and  low  values  given  and  the  beta  curve 
associated  with  the  weight  uncertainty.  This  weight  is  used  later  in  the  progam  for 
determining  the  shipping  cost  based  on  a  per  pound  cost  estimate  for  the  shuttle. 
The  life  cycle  cost  of  each  cost  element  is  then  generated  from  the  high  and  low 
values  assigned  to  it  and  it's  corresponding  beta  uncertainty  curve.  This  cost  is 
then  multiplied  times  the  number  of  units  needed  and  discounted  back  from  the 
time  purchased  to  1985.  The  program  repeats  this  process  for  each  cost  element 
then  sums  the  life  cycle  costs  to  determine  the  total  life  cycle  cost.  The  life  cycle 
cost  is  then  tested  to  find  the  interval  to  which  it  belongs  on  the  histogram  and  is 
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counted  for  that  interval.  The  above  process  then  repeats  itself  for  a  total  of  500 
iterations  and  provides  a  count  of  each  life  cycle  cost  within  each  interval  of  the 
histogram.  The  same  program  was  used  to  compute  the  life  cycle  cost  of  both  the 
digester  system  and  the  conventional  propellant  alternative.  The  program  listing 
with  inputs  and  outputs  is  presented  in  Appendix  G. 

The  program  output  of  the  life  cycle  costs  have  been  converted  to  a 
graphic  histogram,  from  tabularized  data  in  Appendix  G,  and  are  shown  in  Figures 
14  and  15  for  the  digester  system  and  conventional  propellants,  respectively.  A 
composite  of  these  two  charts  is  shown  in  Figure  16.  It  is  noted  from  the  means 
of  these  charts  that  the  life  cycle  cost  of  the  anaerobic  digester  system  does  have 
a  cost  advantage  over  the  conventional  propellants  alone.  The  additional  cost  of 
some  type  of  substitute  waste  treatement  system  to  replace  the  anerobic  digester 
would  raise  the  cost  of  the  conventional  propellant  system.  Even  without  this 
additional  cost  there  is  a  large  difference  [22.1  M$]  between  the  means  of  the  two 
life  cycle  costs.  The  difference  is  due  primarily  to  the  high  cost  of  space 
transportation  and  the  requirement  for  heavy  density  liquid  propellants.  Where 
digester  fuel  requires  the  shipment  of  light  oxygen  gas  compressed  to  2200  psi, 
orbit  maintenance  requirements  calling  for  heavier  conventional  oxidizer  will 
increase  shiping  costs.  The  posibility  that  shuttle  derived  heavy  lift  vehicles  will 
be  available  to  place  larger  payloads  into  orbit  at  less  cost  is  considered  in  the 
lower  range  of  the  shuttle  transportation  cost  element.  A  comparison  of  the  two 
graphs,  in  Figure  16,  shows  that  within  one  standard  deviation  of  either  mean  the 
gap  is  16.5  M$.  This  comparison  reflects  the  significance  of  the  high  cost  of 
shipping  involved  in  transporting  conventional  propellants  over  the  life  time  of  the 
space  station. 


Figure  13.  Flow  Process  For  Life  Cycle  Cost 
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Figure  14.  Life  Cycle  Cost  for  Digester  System 


The  gap  between  the  two  means  can  be  statistically  tested  to  see  if  the 
difference  can  be  attributed  to  chance.  At  a  S%  level  of  significance  the  test 
favors  the  fact  that  there  is  difference  between  the  two  means.  Though  these 
results  may  seem  to  favor  the  digester  system,  it  must  be  remembered  that  many 
other  factors  must  be  considered  by  the  decision  maker  in  a  selection  of  anerobic 
digestion  as  a  means  of  waste  stabilization  and  propellant  production.  Factors 
which  should  also  be  considered  are  the  costs  involved  with  rocket  design  for 
biogas  combustion  or  alternately,  waste  treatment  facilities  other  than  anerobic 
digestion  if  considering  conventional  propellants.  Weighing  into  the  decission  are 
the  additional  benefits  of  reduced  waste  products,  a  fairly  clear  effluent  from  the 
digestor  for  water  recovery,  and  possible  utility  of  the  stabilized  waste  material. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  previous  chapters  have  brought  together  information  and  analysis 
from  three  widely  separate  disciplines;  1)  simulation  analysis  in  biochemistry,  2) 
thermochemical  analysis,  and  3)  cost  analysis.  This  final  chapter  provides  a 
summary  of  the  major  topics  of  importance  presented  in  the  first  four  chapters 
and  their  relationship  to  one  another.  A  brief  overview  of  the  scenario  chosen  will 
be  presented  with  a  short  discussion  of  the  change  in  output  results  that  could  be 

expected  if  specific  input  parameters  were  altered.  With  regard  to  the  scenario, 

the  output  of  the  SLAM  model  is  discussed  and  its  relationship  to  the 
thermochemical  analysis  and  the  sizing  of  the  anaerobic  digester.  A  review  of  the 
key  elements  considered  in  the  cost  analysis  together  with  their  impact  on  the  final 
life  cycle  costs  is  presented  here.  Also  discussed  are  the  items  that  were  excluded 
from  the  cost  analysis  and  their  impact  on  future  decisions.  Based  on  the  results 
of  the  analysis  in  the  previous  chapters,  recommendations  are  presented  for  further 
study. 

Study  Conclusions 

Scenario.  It  is  clear  from  the  literature  review  that  there  is  a  need  to 

reduce  the  costs  of  space  transportation  and  also  a  need  for  waste  treatment  on 

large  scale  space  stations.  The  literature  provides  very  little  evidence  that  any 
research  has  been  attempted  to  tie  these  two  requirements  together  and  provide  a 
cheap  space  propellant  from  wastes.  Therefore,  to  provide  a  stage  for  the  study 
of  such  a  proposal,  a  scenario  was  selected  based  on  previous  work  done  by 
NASA.  This  scenario  was  then  examined  to  extract  those  elements  that  have  a 


bearing  on  the  design  of  an  anaerobic  digester  and  an  impact  on  it’s  biogas  output 
when  considered  for  use  as  a  prcpeliant. 


The  scenario  selected  provided  three  critical  inputs  to  the  analysis  of 
the  methane  generation  and  it’s  value  as  an  orbit  maintenance  propellant.  These 
inputs  were  psuedogravity,  the  number  of  people  aboard  the  space  station,  and  the 
altitude.  Psuedogravity  was  a  requirement  necessary  for  the  anaerobic  digester  to 
operate.  The  separation  of  the  insoluble  biogas  from  the  waste  sludge  is  only 
possible  if  there  is  gravity  to  pull  the  heavy  sludge  to  the  bottom  of  the  digester 
tank.  The  number  of  people  on  the  space  station  was  used  as  an  input  to  the 
digester  model  in  determining  the  amount  of  methane  that  could  be  expected  to  be 
generated.  It  was  also  a  factor  in  determining  the  size  of  the  space  station.  The 
size  of  the  station  was  then  employed  in  evaluating  it’s  mass  and  shape  using  as 
guidelines  configurations  and  parameters  previously  defined  by  NASA  in  their 
studies  of  future  space  stations.  The  mass  and  shape  of  the  space  station,  along 
with  the  orbital  altitude,  allowed  a  determination  to  be  made  of  the  drag  force  that 
would  be  i.ecessary  to  overcome  if  the  vehicle  were  to  maintain  its  orbit.  The 
altitude  selected  was  compatible  with  the  space  shuttle’s  orbit  since  the  space 
station  is  dependent  on  the  Earth  for  a  majority  of  its  resources.  Any  increase  :n 
the  orbit  altitude  would  decrease  the  drag  on  the  space  station  and  increase  the 
percentage  of  the  force  requirement  that  could  be  produced  by  the  generated 
methane.  Any  change  in  population  requirement  would  not  only  change  the 
methane  generation  rate,  but  also  the  mass,  surface  area,  and  consequently  the 
drag  on  the  space  station. 

SLAM  Model.  The  simulation  language  for  alternative  modeling, 
SLAM,  was  chosen  to  model  the  anaerobic  digester  because  of  SLAM’s  capability 
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to  model  a  continuous  process.  In  order  to  build  the  model  it  was  necessary  to 
have  a  basic  understanding  of  the  anaerobic  process  and  describe  it  in  analytical 
terms.  This  description  was  provide  by  work  done  by  Price  and  Cheremisinoff  in 
their  book,  Btogas  Production  and  Utiluation. 

The  only  input  to  the  model  from  the  scenario  was  the  number  of 
people  aboard  the  space  station.  This  input  was  then  transformed  into  flow  rate, 
percent  solids  present,  organic  solids  present,  and  biological  oxygen  demand,  using 
information  obtained  from  NASA’s  previous  studies  which  looked  at  actual 
astronaut's  waste  product  composition  and  the  expected  waste  output  during  future 
space  activities. 

Other  inputs  to  the  model  which  were  not  dependent  on  the  scenario 
were  digester  temperature  and  efficiency.  A  temperature  range  was  selected,  based 
on  historical  data  of  Earth  based  digesters,  which  would  be  most  favorable  to  the 
mesophilic,  methane  producing,  bacteria.  To  allow  for  a  realistic  modeling  of  the 
environment  within  the  digester  the  temperature  was  varied  to  simulate  the  actual 
fluctuations  expected  in  waste  input  loading.  The  efficiency  was  allowed  to  vary 
randomly  so  the  methane  generation  rate  would  reflect  a  reasonable  estimate  of 
reality.  Shifts  in  the  temperature  range  selection  were  made  in  various  runs  of  the 
model,  providing  output  results  that  showed  very  little  gain  in  methane  production 
rates  for  the  increase  in  thermal  energy. 

The  model  was  constructed  using  many  variables  which  have  their  basis 
in  laboratory  experiments  and  scientific  reasoning.  Variations  in  such  things  as 
suggested  solids  retention  time  in  the  digester,  microorganism  growth  and  decay 
rates,  and  expected  waste  concentrations  on  large  scale  space  stations  can  be 
modeled,  but  in  reality  will  require  human  monitoring  and  intervention  to  provide 


peak  methane  generation. 


Thermochemical  Analysis.  Once  an  approximate  methane  generation 
rate  was  determined  it  was  then  necessary  to  predict  how  to  best  use  the  biogas  in 
providing  a  counteracting  force  to  the  atmospheric  drag.  It  was  determined  by 
looking  at  the  amount  of  biogas  generated  from  the  SLAM  model,  and  its  density, 
that  the  biogas  itself  would  not  suffice  as  a  cold  gas  monopropellant.  An  oxidizer 
suitable  for  combustion  with  methane,  in  this  case  oxygen,  was  selected,  and  first 
cut  thermochemical  calculations  were  made  to  determine  the  temperature  of  the 
combustion  reaction  for  input  into  formulas  that  would  determine  the  exhaust 
velocity  and  the  required  mass  flow  rate  to  produce  the  force  required.  These 
initial  calculations  showed  a  temperature  above  which  dissociation  of  the 
combustion  products  would  take  place.  The  dissociation  would  lower  the 
combustion  temperature  and,  therefore,  it  was  necessary  to  incorporate  in-depth 
thermochemical  analysis  into  the  study. 

Thermochemical  analysis  is  a  science  which  is  totally  dependent  upon 
laboratory  data  tabulated  in  JANAF  (Joint  Army,  Navy,  and  Air  Force) 
thermochemical  tables.  The  results  obtained  provide  the  best  estimate  of 
temperatures  that  can  be  expected  with  given  combustion  reactants  and  expected 
products.  Variables  that  could  change  the  output  temperature  prediction  are 
fluctuations  in  the  expected  70/30  methane/carbond  dioxide  mix  and  deletion  or 
addition  of  expected  output  combustion  products. 

Cost  Analysis.  The  cost  analysis  portion  of  this  study  is  definitely  the 
area  where  uncertainty  has  the  greatest  opportunity  to  alter  the  outcome  of  the 
comparison  of  propellant  systems.  One  of  the  factors  contributing  to  this 


uncertainty  is  the  time  frame  proposed  for  a  space  station  of  this  scale.  Ail  that 
has  been  presented  in  this  study  is  certainly  within  the  limits  of  current  day 
technology,  and  as  such  has  been  analyzed  using  current  day  engineering  principles. 
However,  technology  does  change  and  the  systems  the  future  are  dependent  on 
that  technology.  Therefore,  costs  and  weights  and  the  associated  uncertainty 
factors  have  been  selected  with  their  future  application  in  mind  to  limit  as  much  as 
possible  the  underestimation  of  the  system’s  life  cycle  costs.  Technology 
break-throughs  pending,  the  estimates  made  in  this  study  should  provide  a  realistic 
costing  estimate  of  future  propellant  systems. 

Another  factor  that  can  affect  the  cost  analysis  outcome  is  design 
configuration  of  the  anaerobic  digestion  system.  The  design  selected  for  this  study 
was  one  that  is  standard  for  large  municipal  high  rate  digesters.  The  sizing  of  the 
system  was  scaled  down  to  accommodate  the  loads  predicted  from  the  SLAM 
model.  It  is  probable  that  a  system  that  would  be  placed  in  space  would  have  an 
entirely  different  system  configuration.  One  possibility  is  a  "packaged"  system 
that  would  contain  all  components;  tanks,  pumps,  heat  exchanger,  etc.,  in  a 
completely  self-contained  unit.  Maintenance  acquisition  requirements  to 
component  parts  of  such  a  system  would  require  it  to  be  somewhat  bulky,  however 
design  efforts  could  possibly  overcome  this  difficulty.  Another  possibility  is  a 
single  unit,  rather  than  redundant  units  as  proposed  in  this  study,  with  additional 
equipment  actually  tied  into  the  system  but  isolated  from  the  flow  with  valves. 
This  would  provide  quick  response  to  system  failures  through  the  opening  and 
closing  of  the  proper  valves.  The  first  system  provides  convenience  in  packaging 
and  shipping  but  more  front  end  cost  in  design  which  would  be  hard  to  predict 
since  digesters  of  this  type  and  utility  have  not  been  attempted  before.  The 


second  system  has  the  draw  back  of  having  to  continuously  pump  sewage  across 
the  station. 


For  the  system  proposed  the  configuration  allows  for  one  digester  to 
process  the  waste  of  the  entire  space  station  if  necessary.  Most  of  the  time  it 
will  process  the  waste  on  its  side  of  the  vehicle  while  operating  at  half  capacity. 
The  configuration  proposed  also  allows  the  system  to  be  partially  fabricated  and 
assembled  on  the  ground,  with  the  final  assembly  being  completed  in  space.  This 
ground/space  construction  ratio  is  yet  another  factor  which  is  unknown  and  can 
effect  the  labor  cost  element  in  the  analysis.  The  uncertainty  associated  with  this 
cost  has  been  selected  to  hopefully  capture  the  range  of  possibilities. 

The  outcome  of  the  cost  analysis  shows  a  great  deal  of  spread  in  the 
life  cycle  cost  of  the  two  alternatives  presented.  This  spread  is  due  primarily  to 
the  range  of  uncertainty  in  the  transportation  costs.  The  upper  limit  of  the  cost 
per  pound  transportation  element  is  based  on  current  non-government  subsidized 
costs  for  the  STS  system.  Unfortunately  this  cost  is  much  higher  than  was 
originally  anticipated  for  the  space  shuttle  due  to  cost  overruns  and  correction  of 
many  unforeseen  technical  problems  on  a  new  space  transport  system.  The  lower 
range  is  based  on  an  unmanned  heavy  lift  vehicle  proposed  by  Martin  Marietta. 
The  beta  distribution  selected  favors  the  Martin  concept  for  future  space 
transportation.  If  space  transport  costs  remain  as  high  as  they  are  today,  the 
spread  of  the  life  cycle  cost  can  be  expected  to  narrow  while  the  mean  increases. 

The  spread  between  the  means  of  the  two  systems,  approximately  $22 
M,  is  partially  the  effect  of  the  initial  premise  of  this  study;  the  high  cost  of 
space  transportation.  Another  factor  which  is  reflected  in  this  difference  is  the 


increased  man-hours  in  space  there  is  a  good  possibility  that  the  gap  between  the 
life  cycles  cost  of  these  two  systems  will  narrow.  One  indication  of  this  is  the 
lower  cost  per  pound  prediction  for  the  shuttle  derived  vehicles  which  has  been 
used  as  a  parameter  in  the  cost  analysis  in  the  lower  limits  of  the  cost  range  for 
space  transport. 

Additional  benefits  of  reduction  in  waste  have  not  been  considered  in 
this  study,  but  may  be  a  decision  variable  and  are  readily  available  from  the 
SLAM  model.  The  total  reduction  in  waste  processed  through  the  system  is 
approximately  7%.  This  is  the  percentage  of  waste  converted  to  biogas.  The 
percentage  of  effluent  return  to  the  water  recovery  system  is  approximately  84%. 
The  remaining  9%  is  output  from  the  digester  in  the  form  of  digested  sludge. 

The  conclusion  reached  from  this  analysis  is  that  space  generated 
propellants  from  human  waste  show  a  potential  for  lowering  the  cost  of  routine 
orbit  maintenance  for  future  space  stations.  The  digester  system  can  perform  four 
functions;  waste  reduction,  waste  stabilization,  effluent  recovery,  and  production  of 
a  candidate  propellant.  These  features  can  be  tied  directly  into  a  closed  cycle  life 
support  system  and  provide  a  propellant  for  future  space  missions.  At  higher 
orbital  altitudes  there  is  the  possibility,  with  the  reduced  amount  of  air  density, 
that  the  biogas  propellant  could  provide  100%  of  the  required  force  necessary  to 
overcome  drag.  In  the  very  distant  future  there  is  no  doubt,  in  this  author’s  mind, 
that  long  duration  planetary  missions  will  employ  similar  biological  systems  to 
reduce  wastes  and  provide  propellants  which  can  be  stored  for  use  in  attitude 
control,  orbit  maintenance,  and  landing  control  systems. 

Recommendations 

As  mentioned  early  in  this  study,  a  valid  and  complete  cost  comparison 


between  the  anaerobic  digester  and  the  system  of  conventional  propellant  transfer 
cannot  be  made  without  also  considering  the  cost  of  the  system  that  would  have 
to  replace  the  anaerobic  digester  if  using  conventional  propellants  and  an 
alternative  waste  treatment  system.  It  is  recommended  that  the  results  of  this 
analysis  be  coupled  with  a  cost  analysis  of  other  systems  capable  of  and  feasible 
for  space  based  waste  stabilization.  Such  systems  would  include  reverse  osmosis, 
electrodialysis,  and  physiochemical  processing.  Also  included  in  this  study  would 
be  a  weighted  value  of  the  benefits  of  the  other  systems  being  analyzed,  and  of 
the  anaerobic  digester.  For  the  digester  this  would  include  effluent  recovery,  waste 
reduction,  and  possible  alternative  uses  of  the  generated  biogas. 

As  air  density  decreases  with  increasing  altitude  there  will  be  some 
altitude  at  which  the  methane  generated  would  be  capable  of  providing  all  of  the 
counter-force  necessary  to  overcome  atmospheric  drag.  This  altitude  could  be 
found  by  an  iterative  method  of  guessing  the  altitude,  generating  the  drag  force  at 
that  altitude,  and  comparing  the  drag  force  with  that  which  can  be  produced  using 
the  biogas/oxygen  propellants.  A  study  should  be  conducted  to  determine  this 
altitude  while  considering  the  change  in  drag  coefficient  with  altitude  gain.  This 
altitude  would  then  be  analyzed  to  see  if  it  is  within  shuttle  range  and  below 
harmful  radiation  belts. 

The  energy  required  in  the  operation  of  the  anaerobic  digester  system 
was  taken  from  rule  of  thumb  estimations  for  Earth  based  digesters  and  then  given 
a  spread  and  uncertainty  that  would  attempt  to  capture  the  value  of  space 
generated  energy.  It  is  recommended  that  analysis  be  done  in  the  area  of  future 
space  energy  costs  to  narrow  the  spread  and  uncertainty  of  this  cost  element,  and 
this  new  energy  cost  be  incorporated  in  the  cost  analysis  model.  Other  proposed 


changes  are  the  addition  of  costs  associated  with  space  based  oxygen  generation, 
to  be  used  as  an  oxidizer  agent  in  the  biogas/oxygen  reaction,  and  consideration  of 
the  generated  biogas  for  attitude  control  and  other  uses. 
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APPENDU  A 
DEFINITION  OF  TERMS 

Adiabatic  Flame  Temperature.  The  highest  possible  temperature  which  can  be 
reached  given  certain  reactants.  The  reaction  is  calculated  as  going  to  completion 
[equilibrium  is  reached]  and  taking  place  adiabaticly  [without  heat  gain  or  loss], 
(46:473) 

Anaerobic  Digestion.  A  continuous  microbiological  process  that  decomposes 
organic  material,  by  the  use  of  microorganisms  that  are  capable  of  combining 
hydrogen  with  carbon  in  an  anaerobic  [oxygenless]  environment.  From  a  biological 
viewpoint  this  process  can  be  described  as  taking  place  in  two  stages:  [1]  acid 
production,  and  [2]  methane  production.  In  the  first  stage,  acid-forming  organisms 
convert  complex  organic  compounds  into  primarily  volatile  fatty  acids,  such  as 
acetic,  propionic  and  butyric.  In  the  second  stage,  the  acids  are  converted  to 
methane  and  carbon  dioxide  by  anaerobic  bacteria.  Factors  which  must  be 
controlled  for  successful  anaerobic  digestion  are:  [1]  the  complete  absence  of 
oxygen,  [2]  digestion  time,  [3]  temperature  uniformity,  [4]  proper  balance  of  organic 
carbon,  nitrogen  and  phosphorus  in  the  organic  wastes,  and  [4]  the  absence  of 
elements  that  are  toxic  to  the  bacteria.  The  major  application  of  anaerobic 
digestion  is  the  the  breakdown  of  organic  mater  in  the  concentrated  sludges 
produced  from  the  treatment  of  wastewater.  (2:93;  41:1,88) 
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Antithetic  Random  Number.  Random  numbers  generated  which  have  a  strong 
negative  correlation  to  a  previous  set  of  random  numbers  to  produce  a  small 
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variance  in  the  dependent  variable.  This  negative  correlation  is  created  by  making 
the  random  numbers  dependent  of  the  initial  set  of  random  numbers.  Generation  is 
made  by  subtracting  the  initial  set  of  random  numbers  from  one. 

Biological  Oxygen  Demand  [BOD],  The  oxygen  used  by  the  methane  generating 
bacteria  while  they  utilize  the  sewage  as  their  food  and  energy  source.  BOD  is 
determined  by  a  standard  5  day  test  [BOD5]  which  accounts  for  the  amount  of 
oxygen  utilized  by  the  bacteria.  BODL  accounts  for  the  maximum  or  ultimate 
oxygen  demand  on  the  sewage  by  the  bacteria  and  can  be  determined  from  the 
amount  of  organic  material  added. 

Bipropellant  System.  A  rocket  propellant  system  which  uses  two  separate 
propellants,  an  oxidizer  and  and  a  fuel.  The  propellants  are  stored  separately  and 
not  mixed  until  injected  into  the  combustion  chamber.  Most  liquid  propellant 
rockets  use  bipropellants.  (45:204-205) 

Gold-Gas  System.  A  rocket  propellant  system  which  uses  inert  gas  jets 
incorporating  fast-acting  valves  and  receiving  its  propellant  supply  from 
pressurized  cold  gas  supply  tanks.  (45:227) 

Frozen  Equilibrium  or  Frozen  Flow.  A  term  used  in  describing  the  expansion  of  a 
gas  through  a  nozzle  for  the  purpose  of  analytical  consideration.  The  composition 
of  the  product  gas  in  the  nozzle  is  assumed  to  be  the  same  as  that  in  the 
combustion  chamber.  This  results  in  a  conservative  estimate  of  the  exhaust 
velocity  since  in  the  product  gas,  in  reality,  will  experience  a  drop  in  temperature 


and  pressure  and  a  conversion  of  the  thermal  energy  to  kinetic  energy  causing 
phase  equilibria  to  occur  between  the  gaseous  and  condensed  phases  of  the 
constituents  of  the  gas  (called  shifting  equilibrium  for  a  infinitely  fast  shift  rate). 
Consideration  of  an  appropriate  equilibrium  rate  can  be  a  complex  analytical 
problem  which  is  only  possible  if  the  correct  rate  information  is  known. 
(45:183-184) 

Per  Capita.  The  equivalence  in  people  that  it  would  take  to  produce  a  given 
amount  of  waste.  A  term  used  primarily  in  conjunction  with  municipal  treatment 
facilities  whose  waste  input  comes  from  other  than  domestic  areas  (industrial 
waste,  farm  waste,  etc.) 

Solids  Retention  Time  rSRTl  The  average  time  in  days  that  waste  solids  remain 
in  the  digester. 

Specific  Impulse  QJ.  The  thrust  that  can  be  obtained  from  an  equivalent  rocket 
which  has  a  propellant  weight  flow  rate  of  unity.  If  =  F/w  =  c/g.  Where  F  is 
thrust  in  pounds,  w  is  the  weight  flow  rate  in  pounds  per  second,  c  is  the  effective 
exhaust  velocity,  and  g  is  the  gravitational  constant.  Units  of  specific  impulse  are 
pounds  of  trust  per  pound  per  second  of  propellant  flow,  or  seconds.  (45:30) 

Supernatant.  The  liquid  inside  the  digester.  This  liquid  comes  from  the  waste 
water  input  to  the  digester.  In  a  two  stage  high  rate  digester,  the  primary  mixing 
tank  contains  a  supernatant  which  is  a  mix  of  the  solids  with  the  liquid.  In  the 
secondary  settling  tank,  the  supernatant  is  allowed  to  settle  to  remove  the  solids 


and  withdraw  the  remaining  supernatant  as  a  fairly  clear  effluent. 

Volatile  Solids.  Those  solids  coming  into  the  digester  that  can  be  used  by  the 
bacteria  as  food  mater.  This  is  determined  by  weighing  a  sample  burned  at  550  C 
to  determine  the  ash  content  which  is  subtracted  from  the  total  solid  weight. 


APPENDIX  B 

SLAM  MODEL  MAIN  PROGRAM  AND  INPUT  STATEMENTS 


a 


I 


-SObRAfl 

BIVENS I0H  NatTiiMM) 

C3t«W/SC0«!/ATRIB(liii.DS(ll8i  .DDL  ( 100; .  DTNCH.  1 1 ,  MFA .  1STGP ,  NCLMR 
i .  MCRDR ,  MPPWT  .SiHRuM.  UNSET .  VTAPE .  3S 1 1 M, .  5SL  U  001 .  TNEXT .  TNCM .  »U00) 


CQM1DN  QSEM 


EQUIVALENCE 1 NSET 1 1 1  ,QSET  ( 1  • 

NN5ET=ia000 

NLRDR=5 

NP5NT=6 


NPlOT=2 
Call  Elan 


subroutine  eventii> 

COWON.’SCONI/ATRIB'  100)  ,00(100; . DD_  <108/  .CTNO0. 1 1  .UFA . tISTCP, NCLNR 
1.NCR3P,nPRNT.'4NRUN.HN«ET.Nt  APE, SSU88). SSL  alii. TIEU.TNQtUlilBi) 


ElBRCJTInE  :NTlC 

CCN’ONi SCQHi, ATRIBi 130) , OB! 180) .DDL 1 100) .OTNOi, I !, UFA, 1ST0P, NCLNR 
:  .NCPDR.NPRNT.NNRUN, UNSET, NTAPE.SSaM)  .SSLUMi  .TNEXT, TNON.KX  i 100) 


SUBROUTINE  OTPUT 

CGNNQN/SCOni,  ATRIB(180> ,00(100) , DDL (108) , DTNGM.il, UFA, NSTOP, NCLNR 
1 .NCRDR.NPRNT ,NNRUN,NNSET,NTAPE,SS(100; , SSL (100! ,TNE<T,TNON,< j( 1001 

RETURN 


vo.  /./.v.v.y.  -•  .  /  /  /,  v 


»♦♦*♦******* m***3UE ROUTINE  S' ATE*** *♦*♦**♦♦*♦**♦»*♦* *********** 


SUBROUTINE  STATE 

COMMON;  SC  OMi ;  ATR  IB  1 139)  .3DUM)  .DDL 138)  .DTNOU.Ii  ,MfA  .S1STQP  .NClNR 
1 .  NCR2P .  NCPNT ,  NNRuMNSET ,  PiT  APE ,  5S  ( 1 00 ) ,  SSL  '1801  .TnEXT ,  TNOU.  ( 13fli 
LATh  PI: 3. 1415®/ 

EFFICIENCY: 

S3' -i  -  RnQRM  i.S75,  ,0‘5,  2) 

TEMP  CONSTANT: 

85*11)  -  Jill) 

TEMP  VARIATION: 

S3  ill)  *  n(l)-2.8*COS(8»PInTNOIi-0. 1042)1 
LOAD  CONSTANT: 

SS(12)  =  2.6*»X(2)*TNO0 
LOAD  VARIATION: 

55(12)  =  lXi2)»2.6*(l+0.4»COS(8#Pl*TNC0>i 
FLOX: 

SSil>  =  55(12) /100 
SOLIDS  RETENTION  TIME: 

88*2)  =  10 

IF  SS(li>  ,LT.  35)  SS(2)  =  -0. 8+SS i i  1 ) +38 
TOTAL  SOLIDS  : ll'A): 

55  3)  =  0. 13*55 '12) 

VOLATLE  SOLIDS  CALC  HATES  (-3M5H)  (UNITS^S/Di : 

SS'4t  -  SSi3)-0.03*SS!3) 

METHANE  PRODUCED: 

SS»5:  =  0.35E»SSt4itSS(7)*il-i0.0588/ (1-0.333*55(2)))) 

5QL:DS  REMOVED: 

55(6)  -  5S(3)-l.059*SS(5) 

VOLATILE  SOLIDS  L0ADIN6  (FI6  B)  (UNIT3:).S/H3'D) : 

55(8)  =  13/SSI2) 

DIGESTER  VOLUME  USED: 

SS'R)  =  55(4) / S  5 1 B ) 

HtDRAULIC  RETENTION  TIME: 

S3 1 10)  r  oS:9).SSil) 

RETURN 

END 
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t«»*«*m*+**m*»*******lNFuT  ST ATEPENTS******* **♦♦****♦♦**♦ 


GEn, FALLS i sAD, ANAEROBIC  DIGESTER ,5-  2fl> 1983. 

C3l»T:«UQ’JS.«.22.,..3ii7,N:  3  SZFF  EONS.  22  STATE  E8NS,  CAlC  EVER*  1  HR. 
!NTLC.SSUI=0.'28.SS(2!=iB.8,SS(3)=B.3.5St4)=9i.9.S3i5)=B.3,S5ia;=«.0. 
3S<-3>=0.  fl,S3(9i=B.  0.53(101=0. 0,  ;U1)=35,  :X(2>=200; 

RECORD, TNON, TIME,, B.. 0417;  PLOT  POINT  EVERY  1  HR 
VAR.SSO) ,9, FLOW! 

VAR.SSO, R.SRT;  SOLIDS  RETENTION  TIME 
VAR, SSiJI.S. TOTAL  SOLIDS; 

VAR, SS(4),V, VOLATILE  SOLIDS;  CALCULATED  *4B/D) 

VAR, 55(5 i , M , METHANE  GENERATED; 

VAR, SSibt ,G, SOLIDS  REMOVED; 

VAFi ,35*7)  ,E, EFFICIENCY; 

VAR, SSiBl.L, VOLATILE  SOLIDS;  LOADING  FROM  FJ6  8  UG/M3/D; 

VAR, S5!9),D, DIGESTER  VOLUME  USED; 

VAR.SSi  1B> .H.rtRT;  HYDRAULIC  RETENTION  TIME 
VAR, SSitli.T, TEMP: 

VAR,38'12U,TCTAL  INPUT; 

TIMST.SSO.SRT; 
riMSdSSO, TOTAL  SulIDS; 

T;MST,SS(5i. METHANE; 

TIMET, SS i 9 • , DIGESTER  VOLUME: 

"IMST,'5S(10t  ,HRT; 

INITIALIZE, 0,4.0;  RUN  4  DA*S 
FIN; 
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APPENDIX  C 


SLAM  OUTPUT 


jLii-i  ECHO  REPORT 


SIMULATION  PROJECT  ANAEROBIC  DIGESTER 


DATE  5/22/ 1935 


SLAM  VERSION  JuN  34 


GENERAL  OPTIONS 


PRINT  INPUT  STATEMENTS  (ILIST):  rES 
PRINT  ECHO  REPORT  i IECH0) :  rES 
EXECUTE  SIMULATIONS  lllST):  YES 
NARN  OF  DESTROYED  ENTITIES:  NO 
PRINT  INTERMEDIATE  RESULTS  READING  iIPIRR) :  YES 
PRINT  SUMMARY  REPORT  ; I3MRY ! ;  ;ES 


STATISTICS  FOR  TIME  PERSISTENT  VARIABLES 

TIM3T  VARIABLE  IDENTIFIER  INITIAL 

NUMBER  VALUE 


1 

n 

L 


4 


SSI 

2) 

CRT 

3.133BE+82 

SSI 

3) 

TOTAL  SOLIDS 

0. 8000E+32 

55 1 

5) 

METhAnE 

3.245iE+02 

CC  1 
vu  ■ 

4) 

SOLIDS  REMOVED 

B.b342Et02 

55 1 

?] 

DIGESTER  VOLUME 

3,T862E+02 

SS 1 

13) 

HRT 

3, R708E+81 

BY  FALLSTEAD 

RUN  NUMBER 


HISTOGRAM  SPEC 
NCEL  HLQN 


s 


1 


ts: 

is’ 

C'' 

tr 


INITIALIZATION  options 

3E51HNIN6  TinE  OF  SIHUlATIOH  iTTEE3;:  8.88ME-2B 


ENDlN'i  'I-E  QF  :I!1ULAT;.]N  iT'FIN  :  0,i000E+0l 

STATISTICAL  ARRAyS  CLEARED  i.JJCLR;:  :ES 

VARIABLES  INITIALIZED  iJJVARc  YES 

FILES  INITIALIZED  gjFlL):  *E3 


ySET  QSET  St0RASE  ALLOCATION 


DIMENSION  OF  NSET.-QSET  -.NNSET):  10000 

WORDS  ALLOCATED  TO  FILINS  SYSTEM:  0 

words  allocated  tq  indexed  list  tags:  0 

WORDS  ALLOCATED  TO  NETWORK:  15 

WORDS  AVAILABLE  FOR  PLOTS/ TABLES:  9985 


INPUT  ERRORS  DETECTED:  0 


EXECUTION  WILL  BE  ATTEMPTED 
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*  r,  /; ,  v 
-  V 


M  Z  M  M  j> 

J  _  r  n 


SIWiLATISH  SRCJECT  ANAEROBIC  315ES7ER 


B'<  FAuSiEsD 


DATE  5,-  22/1985 


EiiN  Number 


CURRENT  TIME  0.4000E+01 

STATISTICAL  ARRAUS  CLEARED  AT  TINE  0.3000E+00 


♦♦STATISTICS  FOR  TINE-PERSISTENT  VARIABLES** 


MEAN 

STANDARD 

NiNINUN 

NAXINUH 

TINE 

CURRENT 

VALUE 

DEVIATION 

value 

VALUE 

INTERVAL 

VALUE 

SRI 

B.iBe&E+i: 

0.71 98E+00 

0. 1O00E+02 

0. 1202E+02 

0.4000E+01 

0. 1000E+02 

TOTAL  SOLIDS 

8.a7i0E+02 

0. i S73E+02 

0. 4056E+02 

0. 9464E+02 

0. 4000E+01 

3. 1000E-02 

PEThAnE 

0. 1813E+02 

0.47S7E+01 

0.RS07E+01 

0.2774Et02 

0.4000E+01 

0.271RE+02 

solids  removed 

3.4835E+02 

3.17C5E-02 

0. 3017E+32 

0. *52tE+02 

3.4000Ef01 

0.6584E+02 

DIGESTER  VOLUME 

a.527RE*02 

0. 10L3E-02 

0.35R3E+02 

0.3062E+02 

0.4300E+01 

0.7062E+02 

p:,T 

0. 1031E+02 

0.«R82E+00 

0.R700E+01 

0. 1 16&E+02 

0. 4300E+01 

0.RT30E.01 

VWtajT 


♦♦aTt't  -<NL 

i.-cr‘i .**BLtb** 

Sai  I  i 

E'B  >  1 1 

0.  • LuBE^il 

3.2033E-02 

0. i080Et02 

2.0000Et00 

2. T4c4E+fl2 

8.0228E+00 

2.9130E+82 

0.8080E+00 

2.2719E+02 

0.8008E+00 

2.3581E+02 

2. 0808E+00 

8.9277E+00 

0.0000E+00 

8.  1388E+81 

0. 0300E+00 

0.70t2E+02 

0.0000E+00 

0.9780E+01 

0. 0008E+88 

0.3743E+02 

2.2088E+00 

0.72S8E+03 

0.0000E+00 

♦♦TABLE  NUMBER  !♦♦ 

METHANE 

GENERATE 


RUN  NUMBER  1 


VOLATILE  DIGESTER  TEMP 


TOTAL  IN 
PUT 


MINIMUM  i.  18ME+02  0.98B7E+01  0.1082E+01  8. 3593E+02  0.3248E+B2  0.3128E+0 

NAjIMUM  0. 1202E+06  0.2774E*02  0.1300E+01  8.7062E+02  0.3753E+02  0.7282E+0 
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i*"»  w  ■»  «  w  n  ■'«  »  vi  w  ■-  w  vi  w"  — i—  r~r  ' '  u”  i  v-  ■  v*  v"  -  •^rm^  n.»  r-jr'v* -  » w ,T"irT’  T'nr^  v»  s»v  ^  ■*— w  rw^’Tf'.nn- 

T  VftEEft  i+* 

ftiJN 

ftJ“SEft  1 

llhles  o+  plot 

ft=SRT 

0.  ’.000E+02 

0.  1850E+02 

0. 1101E+02 

3.r.5’.E+0: 

*-i£TKANE  SEN 

E  0.V307E+01 

0. 142'5E*02 

3. 1S77E+02 

2.23 2aE+8: 

L=vQlAtILE  SQL  3.U82E+01 

0. 1 1 3&E+0 1 

0. 1 1  ft  IE-01 

0.  i24;E+£l 

[-DIGESTER  VGL  0. 3593E-02 

0. 44603+02 

0.53Z7E+02 

0.il+4E'0* 

T=TE?1P 

0. :243£+02 

0.3374EX02 

0.3300£t02 

0. 3b2 jt+0i 

I-TGTAL  INpUT 

0.3120E+03 

0.4160E+03 

0. 5200E+03 

0.6240E+03 

8  !  10 

15 

20  25  30  35 

40  45 

J0  JJ  '50  JJ 

70 

TC 

/  w 

i  E0  05 

tihe 

0. 0000£  +  03 

*? 

4 

4 

4- 

ft 

0.4170E-01 

R 

4 

T  ft 

D 

i 

3. 3340E-0 1 

■f  J|_ 

I  (1  D 

4 

4 

3. 1251E+38 

ft&TL 

4 

4 

4- 

a.  istBE+aa 

R 

D 

ft  I 

T 

4- 

a.  28S5E+00 

R 

4 

+  ft 

D 

i 

a.  :532£+aa 

R 

4 

4 

A 

B.  2919E+00 

+R 

4 

T  ft 

D 

T 

0 . 333ot+03 

♦  TL 

,1  D 

4- 

A 

0. 3753E+00 

IDTftL 

4 

4 

4- 

0.4170E+80 

ft  (ID 

I 

T 

4 

3.V587E+M 

R 

4 

+  ft 

D 

T 

0.S004E+00 

n‘ 

4 

4 

F 

1  4 

0. 542IE+00 

+ft 

4 

It 

B 

hi 

3.5333E+00 

*  TL 

Ix  .1  0 

4 

+ 

2.s25!E+00 

0  ft  L 

+ 

4 

4 

0.a672E*00 

R 

D 

11  +! 

+T 

4 

2.  '089E+00 

R 

4 

+  ft 

0 

+1 

0.7506E+00 

? 

4 

+  ft 

4 

0,7r::et08 

4  R 

4 

T+  ft 

B 

!  + 

0.S.:i0E+80 

t  TL 

1  1+  D 

4 

A 

0.8'57E+00 

D  7  L 

ft 

4 

4 

4 

3.5174E+0I 

ft 

0 

+1  ft 

+T 

4 

0. ;5ftlE+00 

r, 

4 

4  P 

D 

+  1 

0.  :00iE+8i 

ft 

4 

4- 

A 

0.1042E+01 

*  ft 

4 

i 

ft 

1  + 

3, 13S4E+01 

+TL 

1+  1  D 

A 

4 

i.H2sE*Bl 

D  TL  H 

4 

4 

•f 

3.  lls3E+01 

? 

D 

+  1  H 

+T 

4 

0.  ;20'3E+01 

u 

+ 

*  ft 

r, 

L 

+  ! 

0. 1 2 j l E+0 1 

R 

4 

+  ft 

4 

3. 1293E+01 

+  r 

+ 

T+  ft 

D 

1+ 

8. 1 3  v4E+®l 

*7L 

n+  j 

4- 

4 

0. 137tE+01 

D  1  L 

4 

4- 

4 

0.iil8E*0l 

ft 

D 

ft  +i 

+  T 

4 

C.  !4f9E+0S 

ft 

+ 

+  ft 

D 

+i 

3.3a78E*0l  R 
3,  ."'1 1  £t0 i  R 
3. 3757E+81  p 
3.:7*'5Et01  » 
3..33at+0l  L 
3.33’8E+0l  D 
0.3«20E+01  R 
0. 39&1E+01  R 
3. 4003E+01  R 


43  45  53  55 


OUTPUT  CONSISTS  OF  99  POINT  SETS  i  594  POINTS) 
STORAGE  ALLOCATED  FOR  1426  POINT  SETS  (  9982  IORDS  i 
STORAGE  NEEDED  FOR  99  POINT  SETS  i  693  WORDS  I 
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APPENDIX  D 


!!»TlC.s*iI<=37.: 

SEEDS. 5ol74i8*5t2;/N0: 
SIHjLAIE; 


SIMULATE 
S I  "ill  AT  £ 
SSflULATE 


SEEDS.' -5si7i4845t2  -NO; 
d  I  flL'LATE ; 


SIMULATE; 

SIMULATE; 

SIMULATE; 


% 

N. 


APPENDIX  E 


OGGQfJ 


■Vr,"TV’  TXW  "  7*  ”V^  r\f 


Ai  INDIX  F 

ADIABATIC  FLAME  TEMPERATURE  CALCULATIONS 


I.  Adiabatic  Flame  Temperature,  no  Dissociation 
Reaction: 


CH4  +  202  C02  -  2H20 

The  actual  reaction  of  biogas  with  02  on  a  molar  basis  is: 


.5525[CH,  -  202]  +  ,2368C02  ,5525[C02  -  2H20]  -  ,2368C02 

On  a  per  mole  CH4  basis: 


CH4  -  202  +  ,4286C02  2H20  «-  1.4286C02 


AH,  =  EnpAHfp  -  En  AHfr 

=noo2AHfc02  -*•  nH20AHfn2o  ~  frcm^Hf  CH4  *  n0iAHf0J 

-  1.4286(— 94.054)  -  2(57.797)  -  [(-17.89o)  -  ,4286(-94.054)] 
=  -134.365  -  115.594  -■  58.207 
=  -191.7548  Kcal/moleCH4 

If  191.7548  Kcal  are  added  to  the  products,  then  by  trial  and  error: 
AHr  =  191.7548  Kcal  =  Enp(H-H298) 


-  nH2o(H-H2gg)H20  ■*"  nc02(H~H298)c02 


®  5500K: 


2(64.949)  ^  1.4286(74.433)  =  236.233  Kcal 


@  5000K: 


1WVW7 


i%r^nIr»,nir 


2(57.829}  *  1.4286(66.753)  =  211.021  Kcal 

®  4500  K: 

2(50.777)  +  1.4286(59.122)  =  186.016  Kcal 
®  4600K: 

2(52.181)  hh  1.4286(60.644)  =  190.998  Kcal 

Therfore,  the  adiabatic  flame  temperature  without  dissociation  is  approximately 
4600K.  At  this  temperature  dissociation  would  be  taking  place.  A  thermochemical 
analysis  involving  all  the  possible  products  must  be  accomplished  to  determine  the 
proper  adiabatic  flame  temperature. 
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II.  Adiabatic  Flame  Temperature.  Dissociation 
Reaction: 

GH4  -  202  -  ,4286C02  -»  aH20  -  bOH  ■+■  cO  +  d02  eH  +  fH2  +  gC02  +  hCO 
where; 

a,b  c,...  =  Number  of  moles  of  particular  product 
Writing  mass  balance  equations  for  H,  0,  and  C: 

H:  2a  ■+■  b  ■+•  e  2f  =  4 
0:  a  +  b^c-2d-2g-h  =  4.8572 
C:  g  -  h  -  1.4286 
The  equilibrium  equations  are: 

1/2  H2  H 

1/2  02  ^  0 

1/2  H2  h-  1/2  02  ^  OH 

H2  -  1/2  02  H20 

C.oi,d  -  1/2  02  ^  CO 

By  using  these  eight  equations  and  the  molar  equilibrium  constants,  Kn  below,  a 
solution  for  the  adiabatic  flame  temperature  can  be  obtained  by  trial  and  error. 

Kn,e  =  nH/nH21/2 
Kn,b  =  nOH/no2l/2nH2l/2 
Kn,a  =  riHzo/nHzfioj^2 
Kn,c  =  n0/n021/2 
Kn,h  =  nC0/n02,/2 


These  molar  equilibrium  constants  are  related  to  pressure  equilibrium  constants, 
Kp,  which  can  be  looked  up  in  thermochemical  tables  f 4],  by  the  following 
relationships: 

Kn,e  =  (Kp,e)(n/p)1/2 
Kn,b  =  Kp,b 
Kn.a  =  (Kp,a)(n/p)'^2 
Kn,c  =  (Kp,c)(n/p]‘/2 
Kn.h  =  (Kp,h)(n/p)"1/2 

where; 

n  =  total  number  of  moles 
p  =  pressure  in  atmospheres 

These  relations  can  be  rewritten  in  terms  of  the  equilibrium  mole  numbers: 

nH  =  (Kn,e]nH21/2 

noH  =  (Kn,b)nH2,/2no21/2 
nH2o  =  (Kn,a)no21/2nH2 
n0  =  (Kn,c)n021/Z 
nco  =  (Kn,hln021/2 

Substituting  these  relations  into  the  mass  balance  equations: 

NH  =  4  =  (Kn,e)nH21/2  -  2nH2  *  (Kn,b)nH21/2n02,/2  -  2(Kn,a)n021/2nH2 
N0  =  4.8572  =  (Kn,c)n02l/2  -  (Kn.bW'V'2  -  2n0a  - 
(Kn,a)n021/2nH2  -  2nC02  +  (Kn,h)n021/2 
Nq  =  1.4286  =  nC02  (Kn,h)no2^2 

where; 


N,  =  Total  number  of  atoms  of  that  species 
Using  Nh,  no2^2  can  be  solved  for  in  terms  of  nu^: 


no2l/2  =  (NH-(Kn1e)nH21/2-2nH2V((Kn,b)nH2'/2no21/2^2(Kn,a)no2l/2nH2)  =  A 


Substituting  into  N0  and  Nc: 

Nq  =  2no2  no21/2((Kn,c)  (Kn,b)nH21^2  -*■  2(Kn,a)nH2) 

2nCQ2  ■+■  (Kn,h)no21/2 


or; 

N0  =  2A2  -  A((Kn,c)  -  (Kn,b)nHa,/2  -  2(Kn,a)nH2)  - 
2nco2  (Kn,h)A 

and; 

Nq  =  nG02  A(Kn,h) 

The  equations  for  N0  and  Nc  are  now  in  terms  of  nH2  and  known  constants  since: 

nG02  =  1.4286  -  (Kn,h)A 


The  computer  program  which  follows  is  a  modification  of  a  program 
written  by  Capt  Robert  Dimmick  (6)  (dealing  with  two  reactants  and  six  products) 
and  iterates  on  a  guess  of  the  one  species,  nH2,  which  the  equations  are  now  in 
terms  of,  at  an  estimated  adiabatic  flame  temperature.  The  actual  adiabatic  flame 
temperature  is  obtained  by  an  analysis  of  the  heat  energies  as  was  done 
previously  when  dissociation  was  not  considered.  For  this  case: 

Qi  - 

The  heat  necessary  to  change  the  reactant  from  a  liquid  to  a  gas  at  the 
boiling  temperature  equals  zero,  since  reactants  are  already  in  a  gaseous 
form. 

Q2  - 


The  heat  necessary  to  raise  the  reactants  from  their  boiling  temperatures 
to  298K  also  equals  zero,  since  reactants  are  stored  at  room  temperature 
(the  reaction  temperature). 


The  heat  of  reaction  which,  as  stated  before,  is: 


Q3  - 

AHr  =  2npAHfp  -  En  AHfr 
which  in  the  program  is  written  in  the  form: 

EnpAHfp  -  58.207  Kcal/mole 

Q4  - 

The  heat  necessary  to  raise  the  products  from  the  reaction  temperature, 
298K,  to  the  combustion  temperature  (Tc),  and  as  stated  before  as  Q,„ 
equals  Enp(H~H294). 

The  program  uses  the  heat  balance  equation; 

Qtot.i  =  Q1  •+•  Q2  +  Q3  +  Q4  =  0 

to  determine  if  the  appropriate  adiabatic  flame  temperature  was  chosen.  If  not 
the  program  calculates  a  new  temperature  and  the  mass  balance  equations  are 
again  iterated.  In  the  program  Q3  has  a  built  in  sign  convention  that  negative  is 
heat  out  so  that  Qi  Q2  Q4  =  -Q3  maintains  the  proper  sign. 

The  computer  program  and  the  last  few  iterations  of  program  output 

follow. 
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1833  sc1*  .-?Q5Rmh  -fiswl* 

:3i3  it1 

.013  "Ei  “M3  croara«  salves  tng  eoui.iMji  comsucr  orcDie*  ter 
.303  RE"  a  rents.  'a:1  at.  !terjti:ns  :s-  ♦:nd  values  :or  tre 
:i 43  -E*  isiacat::  nsae  tefseratu'e  1  Ft  , 

.353  PEN 


izi  FEN 
1870  REN 


.NFuT  F>  OUT  I NE  5 


.383  ?=138»14.®9  : REM  P  is  the  caebustor  oressure  in  ate 

1093  SlIN-0  : REM  81  is  the  <cai  req  d  to  ranee  free  Li  a  to  Sas 

1130  Q2IN=a  : REM  82  is  the  kcai  'SQ  a  to  raise  Jeep  to  298  K 

1110  DIM  PROD l3, 12,3) ,HRE  3» ,MQLE ;3! 

1120  FEU  ;  OANAF  TABLE  DATA  IS  ENTERED  IN  ThE  ARRAY  PROD (CHEM, TEMP . J ) 

1133  REN  ;  THE  FOLLOWING  VALUES  ARE  APPLIED  TO  THE  ARRAY  DIMENSIONS 

1 1 90  REN  CHEN  --  CHEN  1  =  MONATOMIC  MtDRjGEN  ihi 

1150  REN  CHEN  2  *  DIATONIC  H'BRG&EH  iHZi 

1130  REN  CHEN  3  =  MONATOMIC  Q* t GEN  lOl 

1173  FEN  CHEN  4  -  DIATONIC  CXYGEN  -02! 

1133  5EN  C*EH  5  *  h*5SOiYL  jCHi 

1193  RE-  CHEN  s  =  »ATER  >Oi 

11=3  ;EN  .  CHEN  7  =  CARBON  MONO; ICE  iCOi 

11s’  FEN  THEN  8  =  CARBON  2 IOx IDE  C02- 

1  233  RE*1  TENP  —  TENPEPATuRE  OF  PRODUCTS  IN  DEGREES  * 

1 1:0  CEN  J  --  J*1  S'ORASE  Or  TENPERaiuRE  C'ATf* 

1323  FEN  J=2  STORAGE  OF  H  -  hi293!  DATA 

1230  REN  :  J*J  STORAGE  OF  l06  aFi  DATA 

1343  REN  THE  ARRAY  HREiCHEN;  CONTAINS  HEAT  OF  REACTION  DATA  >T=298  *) 

1253  PEN  -  THE  ARRAY  NQitlCHEfli  CONTAINS  THE  EQUILIBRIUM  NOLES  OF  PRODUCTS 

I2a0  FOR  CHEM=  1  TO  8 

1270  PQR  J=1  TO  3 

1280  FOR  TENP=0  TO  9 

1290  READ  PfiODICHEH.TENP.J) 

1730  NEXT  TEMP 


NEXT  Chen 

FOP  lhEN-1  TO  s 
REhD  HRE'ChENi 
NE  1 T  CFEN 


!.■  e  cNiEr*  iNI  IAl  jlicaS  ?QR  THi  ~  •  L  ~wH*E  ’c*£'E:.hT.RE 
ITS*  REN  T-E  5UEES  HOST  BE  IN  THE  BAsgE  :2«?i-T0c3  ZE39EE  t 
IT’D  PEN 

:-j0  input*  ester  initial  mss  pgr  Tc  adiasatici  i-jtc 

uta  T:a=!Tc-:000)«.ai 

1*20  T*SPsi!)T  ;T£0i 

»4’.3  EC4lE-TC0-TENP 

1448  PEN  T„E  .'AR I  ABLE  SCALE  IS  USED  TQ  DO  A  LINEAR  INTERPQLAT[0N  OF  DATA 
1458  REN  CONFUTE  ko  5  FOR  THE  CHErllCALS 

:4i0  K  PH= 10- iPROD 1 1 , TEMP . 3) +SCALE*  tFROO  U . TEflP+1 . 3 ) -PROD 1 1 . TEMP . 3) ) ; 

1470  KPP2= IB" I  PROD ( 2 . TEHP , 3 ) +SCALE* ! PROD ( 2 , TEHP+ 1,3) -PROD ( 2 . TEHP . 3 ) I > 

1480  1 PO=  1 8” ( PROD (3 , TEMP , 3 i +SCALE*  i PROD (3 .TEMP* 1,3) -PROD ( 3 , TEHF , 3 ; i ■ 

1498  k  P02= 1 0* ( PROD ( 4 ,  TEHP ,  3l  4-SCALE*  iPROD  14,  TENF’+l ,  3  >  -PRQD  l4 .  TEMP .  3>  •  > 

1588  kPQH= 1 0* ( PROD ( 5 , TEHP, 3 1 +SCALE* t PROD  1 5 , TEHP+ 1 , 3 ) -PROD ( 5 , TEHP ,3  f  f ) 

1510  KPH20= 13* '.PROD (6 , TEHP , 3 I +SCALE* ( PROD ( 6 , TEHP+1 . 3 ) *f'RQD (6 , TEHP . 3 )  1 1 
1513  kPC0=18"(PROD(7JENP,3)+SCALE*<PRODi7,TEf!P+l,3i-FRGD(7tTEt1P.3)!) 
1517  KPC02* 1IA  <  PROD (8 , TEHP , 3 1 +SCALE* ( PROS  >8 , TEHP* 1 , 3) -PROD • 8 . TEHP , 3 » ) i 

1328  SEN  ENTER  INITIAL  DATA  GUESSES  FOR  THE  EQUILIBRIUM  ITERATIONS 
1570  NE=3.9  '.REN  NE  is  the  guess  For  tne  total  soles  of  products 

1548  NH2=.4  :R£H  NH2  :s  tne  guess  For  the  soles  of  H2  product 

1550  _KS0=4. 3572  :REH  LHSO  is  tne  total  soles  of  0  in  reactants 

1355  -H3C=l. 428ft  :PEN  lHSC  is  the  total  soles  of  C  in  reactants 

1568  REN  CONFUTE  <n's  FOB  THE  PRODUCTS 
15*0  •NH=KPH*3QR(NE/Pi 

losa  kNOsKPO*S8" 'NE- P) 

:5':0  'UGn-‘FQH 

-U08  i'N'H2Q=rPH2j*S3R(P/NE) 

1;05  f  NC2=P  PCO»SfiR iP/ ME) 

till  REN  .•  EQUILIBRIUM  CALCULATION  ITERATIONS  » 


t:i  A*  i  4-KNH»SQR  t»H2» -2*MH2)  /  4kNOH»SQR  !NH2>  s-2*KNH20*NH2> 

a48  N02=A'2 

o50  N0=kN0»SOR(N02) 

660  HH=kNH*SQR(MH2) 

oT0  NC:H=tN0H»SQR;NH2)»SQRiNQ2) 

630  NH20=KNH2C«NH2s30R!N02I 

;85  NCu=*  NC0*SQRi.N02i 

6dT  NCG2=L-SC-NC0 

690  PEN  NE  NON  CHECK  THE  ATONIC  BALANCE  CDR  Oarqen 
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r  v ;  v~.'  T  v  ■ 


1/ 


.  -r. 


*  . 


.730 

^  :  'j  J-rt* ! <  NO*KNun*itfh  iNirZ  **  '«•"  ~  { 

1 J  _ 

17:0 

iF'QK=LHbU“:U  G.\ 

r:a 

IF  ABSlEOl-  .001  7fi£u  in 

1”0 

PRINT  J;ERQR 

i'40  FEN 

HE  HAKE  H  NEW  GUESS  FOR  NH2  AND  1'ERATE 

.'50 

Nh2=NH2-ERjR*. 1  :  3=8"! 

!*;0 

SCTli  1:70 

1*70  :'EM 

*E  NOW  CHECK.  THE  HOLE  BALANCE  FOR  7h£  31, 

£N  NH, 

1730 

NEE=NH+NH2+N0+NQ2+NGH+NH2G+NCG2+NC0 

1790 

ercr=ne-nee 

1-700 

IF  hBSiERORi  ..0001  THEN  1370 

Ibid 

F'RINt:PR!NT:PR!NT 

IS  20 

PRINT'  NE=*;NE:‘  NEE=S;NEE 

1370  REM 

*E  ADJUST  NE  AND  ITERATE  ASA  IN 

1340 

NE=.5*U.5*NEE+.5»NE) 

1850 

60 TO  I5o0 

1340  REN 

OUTPUT  ESUIlIBRIUH  CONDITIONS  ;;; 

1870 

LPPINT"  THE  TOTAL  NUMBER  OF  PRODUCT  HOLES  IS: 

1330 

lPRINT*  THE  NUMBER  OF  MOLES  OF  H  IS:  * 

:NH 

1390 

SPRINT'  H2  IS:  “ 

;NH2 

1900 

LPRINT*  0  |S:  * 

;N0 

;510 

-PRINT"  02  IS:  1 

;NC2 

1920 

-PRINT"  OH  IS:  ' 

:N0h 

ifC3 

.PRINT"  *20  IS:  " 

:NH20 

:  977 

lPRIN'*  CO  IS:  1 

:NC0 

1--77 

-9RI.T”  C02  IE:  " 

:  NC02 

1;  40 

lPRINT 

1:50  9EH 

LOAD  FINAL  HClE  NUMBERS  > 

1"  :0 

MOLE ■ 1 ’ =NH 

.9-0 

MGLEi2/=NH2 

1 :30 

NQlE(3)=NG 

'.990 

HOLE i ♦ ) =N02 

2300 

N0LE<5) =NQH 

7310 

M0LEl6i=NH20 

231.3 

HOLE  1 7  >  =NC0 

73l7 

MOlE t8i :NCQ2 

.320  :'EH 

FINAL  CALCULATIONS  CF  3 

.370  9EB 

iCOuT  REPRESENTS  ThE  HEAT  OF  REACTION  AT  I 

i:8K  F 

7040  R'EM 

mole  numbers  of  products 

»nco:-nc] 


2030  FEN  Q4 1M  RE°FESENTS  THE  -*EAT  NECESSARY  Tj  RAISE  THE  PR0DuC75  FROM  THE 
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-  -  *  '-  *  w  .  -  *  i-N  .  •  ;  •  ■  p  ■%  >  *_ .  f  .  1, ?-WV  A  J  r  jf*.  ■  "_ji  ‘j>  *jt  -A  .  y  .  w* .  ^ -  A  *-«  W  iA«  a  a  *  Jjf „  W  jf ,  kV  g 


I ’ 20  $CAlEH=SCAlE*  PROD  t ChEP* .  TENP+  ROD  anEN .  T£-P , 2, , 

21 13  94  INDOLE  >  CKEH  «  <  £ROD  (CHEfl.  'EPF .  2 1  *5CAlEHi  *94  IN 

:::b  e t t  c.-eh 

21:0  for  ohen=i  to  a 

21 40  j3DUT=- • NGlE  iCKEH) *HRE  iCHEN) ) *Q30UT 

2150  NEK!  CHEM 

2155  53uUT=Q3QuT-5S. 20? 

2160  REM  OUTPUT  9  CALCULATIONS  '» 

2 1 70  LPRINT"  FOR  A  COMBUSTOR  TEMPERATURE  OF  DEGREES  s* 


2180 

lPRINT* 

91 

IN 

IS: 

*:9!IN:" 

Kcal* 

21-50 

lppint* 

92 

IN 

13: 

“:921N;“ 

Kcal* 

2200 

LPRINT" 

33 

OUT 

IS: 

“ : Q30UT ; 

*  teal ' 

2210 

lPRINT* 

64 

IN 

IS: 

*:S4!N; " 

Kcal* 

2220 

LPRINT 

2270  QTQT I N=3 1 1 N+Q2 1 N+Q4 IN 

2235  REN  QTOTO'JT  REPRESENTS  THE  SUN  OF  THE  HEAT  OF  REACTION  FOR  THE 
223s  REN  PRODUCTS  AT  298K  MINUS  THE  HEAT  OF  REACTION  FDR  THE  REACTANTS 
2237  REN  AT  2^3* 

2240  QtOTjUT =Q30UT 

2250  DElTAQ-QTOTIN-OTQTOUT 

2-30  lPRINT"  total  9  in  IS: “jO'CTIN; *  K:a2 * 

22'0  l?RINt'  TOTAL  9  cut  IS: 1  :QTOTOJT : 4  K;ai‘ 

2230  .FftiNT*  the  8  NET  IS: * tDELTAQ; J  Kcal' 

2290  lFPINT:LFRINT:LPRINT 

2300  REN  RESET  VALUES  TO  ZERO,  BESIN  ANOTHER  TEMP  ITERATION  >' 

2310  93003=0  :  Q4IN=0 

2320  REN  <>■;  CHECK  FOR  ADIABATIC  TENP 

2330  IF  ABS(DELTAQk.0B1  THEN  2360 

2340  TC=TC-DELTAQ*10 

2350  SOTO  1410 

i3s8  END 

.ji0  R'tN  Jf*NAF  TABlE  DhTA  ARRAY S  -  • 

1330  REN 

2390  REN  THE  DATA  IS  ARRANGED  8t  CHEMICAL  AND  CONSISTS  OF  'ENPSRATJPE. 
2400  REN  H  -  H298,  and  LOS ^D)  LISTED  SEQUENTIALLY  F0R  THE  PRODUCTS 
2H0  REN 
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2120  f.EM  '  DATr*  -OR  H 

:*:a  data  3000,3100. 3:00. 3300. 340a. 2:00. 3500.3:00. :50a. 2^0 

1140  SATA  iJ. 423. 13.92. 14.417. i4.?l4. 25. 4i. is.W. l3.4i4.1a.T0i..’.:?S.l;.S!:5 

Uli  DATA  -.303. -.671, -.553,-. 439, -.232. -.231. -.135. -.044.. 011.. 123 


2460  REN 
2420  EE* 


jATA  FGR  H2 


2480  DATA  3000. 3100. 3200. 3300. 3400,3508.3600,3700.3800. 3900 

2490  DATA  21. 21. 22.098. 22. 9R2.23. 891, 24. 794. 25.703,26.614. 27. 535. 28. 457. 29. 365 

2300  DATA  0,0. 0,8, 0,0, 0,0, 0.0 

2:10  RE* 

2523  RE*  DATA  FOR  0  > 


2530  DATA 
2540  DATA 
2550  DATA 
2560  RE* 
2570  RE* 
2530  DATA 
2590  DATA 
2600  DATA 
26.0  EE* 


3000.3100,1200,3308,3400.3500.3600,3700,3800.3900 
13.552, 14.023. 14.524, 15.026. 15.529, 16.033. 16.537, 17.043. 17.549, 18.057 
-.949, -.805, -.670, -.543, -.423, -.310. -.204. -.103, -.007,. 084 

'  DATA  FOR  02  > 

3880. 3130, 3200, 3300, 3400, 3500, 3d00, 3700, 3800, 3900 
23.446.24.403, 25.365, 26.331.27.302.28.276,29.254, 31.236.31. 221, 32.209 

0.0. 0,0, 0.0, 0.0. 0,0 


2e20  REN  <.  DATA  FOR  OH  ; 

2  -.30  DATA  3000.3100,3 200 , 3300 , 3400 . 3500 , 3604 . 31 44 . 3844 ,3344 

2640  DATA  21. 404. 22. 283. 23. 166,24. 052. 24. 94. 25.832. 26. 726, 27. e23.28.522.29.i23 

2650  Data  ,  1078 . . 098. .  1  to . . ;33. .  158. .  165. .  1 79, .  192. . 205, . Il6 

2560  CE* 

26'0  3E*  Data  for  hi q  > 

2680  DATA  3000. 3100. 3200. 3300, 3400. 3500. 3600, 3700. 3800,3900 

2690  DATA  30.201. 31. 535.32.876, 34.223,35.577,36. 936,33.3.39.669, 41. 043.42.422 

2*00  DATA  1.343. 1.201. 1.067, .942,. 324. .712.. 607,. 507.. 413,. 323 

2'10  REM 

2711  RE*  DATA  FOR  CO  > 

2712  DATA  3000.3100.3200.3300.3400,3500,3600.3700,3800.3900 

2713  DATA  22. 357, 23.248. 24. 139. 25.032. 25.927,26.822. 27.719, 28. 617, 29. 516, 30. 416 

2714  DATA  6.407,6.336.6.269.6.206,6.145.6.088.6.034.5,982,5.933.5.886 


2715  *EN 
2?ib  6E* 
2117  DATA 
.718  DATA 
2714  DATA 
2 7 20  f'E* 


DATA  FOR  C02  ; 

3000, 3 100, 3100, 3300. 3400, 3500, 3600, 3700, 3800, 3900 
36, 535, 38. 024, 39. 515, 41. 018, 42. 507, 44. 006, 45. 508, 47. 012, 48. 513, 50. 027 
5.842,6.668,6.458,6.260.6.074.5.898,5.732,5.574,5.425.5.283 


■ai&i8S3iii&fliiii^^ 


'V 


i  GUt  *  -■  -  «  .•  I  j  ?:£di 

h£  5  NET  15:  1.37a83E-i3  Heal 


T«E  '3TAL  NUMBER  uF  PRODUCT  HOLES  15: 
TnE  NUMBER  C-F  HOLES  0C  H  IS:  . 2 i 3394 


«2 

IS: 

.40249 

a 

IS: 

.09.57504 

02 

re. 

.  j  ■ 

.127308 

OH 

IS: 

.317272 

120 

T  C 

1.33218 

■20 

IS: 

0 

C02 

IS: 

1.428a 

H2Q 

IS: 

!.33ilS 

CO 

IS: 

0 

C02 

IS: 

1, 4i.3o 

FOP  A  OOMBuSTOR  TEMPERATURE  OF  3383. la 
31  In  IS:  a  Real 

52  IN  IS:  0  Heal 

93  OUT  IS:  133.315  Kcal 

jl  IN  IS:  133,314  kcal 

•07AL  3  m  IS:  133.314  Heal 
TCTAc  3  out  IS:  133.315  <cal 
'he  9  NET  1 5: -1 . 39363E-03  Heal 


Th£  TOTAL  NUMBER  OF  PRODUCT  MOLES  IS:  3 
THE  NUMBER  OF  MOLES  OF  H  IS:  ,213483 
92  IS:  .40249T 
0  IS:  .0^7552 
02  IS:  . 12T31 1 
Oh  15:  .31728 

HtO  .5:  l.jjklo 

30  IS:  0 


5E6REES 


91821 


OR  A  COflBUbTGK  T;i1PzRh  ,  ijn'E 


In 

i- ! 

a  * 

-■  -  a  i 

L 

.n 

:5: 

a  * 

cal 

OUT 

Is: 

133. 

313 

*cal 

4 

IN 

IS: 

:i4 

fcai 

0* 

4i_  ; 

j  :n 

*5 

\  T" 

,314 

seal 

T 

< J  ' 

AL 

3  cu 

*  7  C 

w  i  J 

«  73 

1  U  1  j 

seal 

ME 

3  1 

MET 

IS 

9.4s#i5E- 

■24 

APPENDIX  G 


Cv. 


I  ** 


MONTE  CARLO  PROGRAM  FOR  LIFE  CYCLE  COST 


Program: 

C  ASui SIT iOri  ■_ _ £T I *«S  Rjj  InE  Tj  DETERMINE  LIFE  I t CLE  CG5T 
C  OF  ;PALL  3 h£ ED  mNAE-OBIC  iioEaTtn'. 

C 

C  .ARIriBLE  GEFF INI T IONS: 

C  ♦*«RRAt  VARIABLES** 

C  Zrtil)  -  AIStiEST  POSSIBLE  COST  GF  ACQUISITION  COST  ELEMENT 
C  CLil)  -  luNEST  POSSIBLE  lOlT  OF  ACQUISITION  COST  ELEMENT 
C  ISiI)  -  BETA  UNCERTAINTY  CURVE  TYPE  FOR  COST  CF  ACQUISITION  E-EMEN' 
C  IfR‘1)  -  (EAR  IN  ACQUISITION  PH-,SE  !H«T  COST  ELEMENT  IS  PURCHASED 
C  Nil;  -  NUMBER  OF  ACQUISITION  COST  ELEMENT  PURCHASED 
C  Hhil)  -  HiSHEST  POSSIBLE  nEIShT  OF  ACQUISITION  COST  ELEMENT 
0  *L i I  -  uG»EST  POSSIBLE  NEIGH!  GF  uCQuISITISn  cost  ELEMENT 
L  IWB  ( 1 1  ■  BE  *  A  UNCERTAINTY  LUFvE  cCR  <EI‘jHT  Q-  ACQUISITION  ELEMENT 

u  LHlili  -  HIGHEST  POSSIBLE  COST  GF  OPERATIONS  COST  ELEMENT 

C  CLU ' i  ■  -  iO»EST  POSSIBLE  COST  OF  OPERATIONS  cost  element 

I  ;§£■•!»  -  BETA  UNCERTAINTY  CuRvE  7fPE  PGR  COST  OF  OPERATIONS  ELEMENT 

:  CI.l!  -  COST  VALUE  GF  UPPER  RANGE  CF  COST  INTERVAL  l  ON  HISTOGRAM 
C  ICO JN" ■  I j  -  NUMBER  CF  COST  ELEMENTS  IN  CCST  Interval  I  ON  HiSTCEM 


i 


i 


i 


El.E  -  w-HEal  POSSIBLE  7E'-L  LlC  E2:r 

*  -  ><;.ShT  jF  A ■_ Q u i : I  . 0 *< ;  COST  .  _E.';;  EE7-*  ■-‘•EEF"-!l'iT 

LL;*  ■  Jl-L  “LJLi.EITDN  _ll  -.3  _;.Ns  :t  -  „’.i,E' 1  - ."i ■< 

..;'C  ■  -.'-L  GPERAT!’3(lS  .El  ECS'  _i>c  Bi  -  .NEE-'-!  *T r 
■L'.-  *  jT-L  LL’C  Ui-ING  svh  JN.FF  I  NT'" 

*•:  -  t;-hl  of  au  leg  s  over  ter  iterations 

S'lE E  -  CF  Sa-SES  E;  'LCC 
*£-*  -  ’*£-*'»  Ec  lEE  S 

~  :  -*\wm-.£  Jt  • •  I-jN  jF  L1L  3 

si-ensign  :■>  eo  ,  el  :t  ,  ibi :i‘ .  h  :e,- .  j  chi •  g ai. 

liXZ'ji'.  i .0-!,  .  „l  100  .  4H',j0),  l-_l_0', 

\:ti  E0.; 

;E*l  Mt«N 
E=0.a 
:T.EE=0.0 
.«E-N=0 

.  Lj=0 

-  i‘iE=0,0 
:  .!iNC=0 

EL--3.0 
ELEB=0.0 
E^=0.3 
Eru:0. 0 
7»L  =  0.0 

r4h  =  0.0 

hLL  ..REN  •il*0ri '  4,  I  N 
CALL  E'FEN  '1.  3DAFA.  in 
EAlL  EPEN  5,  LiT : 

READ  'C.Ej  Nl*.  1tE-.  /•!.  I'E^R,  li-ST  *.  EP.  N.MC. 

'.I <E-R0.  SEES 

2  -O'- Hu"  M.  4 ,  14,  [‘,  IP,  - E ,  I A .  .4,  F* .  E 
«'.Tl  ’...  ‘<L'‘i  I“.C,  .  lE-R,  --liT'A,  jR,  \jN0, 

i  I -E-PE.  SEE j 
LE  E  I  -  '..IK 

e  ;  e  juM 7  ■ : 1  =0 


1 1  *  r  ■  c  ” 

j  -  -  j 

E..EMEN7  BAT 

A  FOR 

ACSL.a'i  ■  ;lN 

PEhS  e, 

,4/  t EL ■ I ■ . 

CH(  I) , 

,  i  B  ’  I  ‘ .  N 1 1 

*  •*  W  •  \  t 

•h '  •  ,  j 

iWB- 1  ■ .  1=1. 

NUM) 

FORMAT 

iFa.0.  Fs.0, 

IE,  i 

13.  12.  f-,0. 

ci . 0 ,  12 

»R!TE  ; 

[ .4!  (CliIi. 

Ch  :  I 

.  ISui ,  n  •  I 

■ ,  .  r  R  • . 

*rt  i  i :  ,  ] 

I  <4  B  <  I  .  :=1 . 

NUM' 

:.r 

element  d«' 

A  FOR 

ULtMi  I'ONS 

- 1  *D  4 . 

,j,  LL'Jlif. 

uhl  1 1 

[..  i ee 1 1 ■ ,  ; 

d.NLMQl 

FORMAT 

iF5.0,  -'9.0. 

12 

rlrUi  |! 

1,5)  iLLii'.Ii 

,  lmI 

i  Ei .  ibElI, , 

i-i  ,N'jMU) 

i^RAND1; 

iEEDl 

'>Er>E 


;-'jH  T:„T«L  LJi  -*ND  ^E.3«T5 


S: 

»*. 

W\ 


00  3  i  1  1.  nuM 
'Nl  =  7NL  +  KLvil 
'<H  =  iHH  -  Nrtli) 

5  CONTINUE 

l  -.OMRijT;.  -* I  ->‘*l  {.G  V$u]"t  OF  LIFE  CyClE  IQET  AND 
-  t  vSljE  J1"  t-Ch  INCREMENT  OF  HlSTOfiRAM 
C 

DO  a  I  =  1.  NUl* 

N' 18!  -  TNi 

CAL  •  M ( I J *CL I! » .  it DR+ 1  * ** 1 1 1* E6R+ 1 V R C I )  >  i 
CL->  =  CAL+CLA 
N « 1 8 !  =  TKH 

CAM  =  'Nil! *CH » I ) > / i <  DR + 1 ) *  * { I YE AR+ 1 1 R 1 1 ) ) ) 
CHA  =  CAH+CrtA 

6  CONTINUE 
L=LASTYR-I¥EARO 
DO  7  J  1  l,  L 

DO  7  1  =  1 ,  NUMO 

COL  =  CLO ( I ) / l (DR+1 i *♦ t I YEARC+J ) i 
ClQB  =  CCl+CLQB 

COh  5  CHI i I ) /  i(DR+li»*iIrEARO+J!) 

CHO  s  COH+CHQ 

7  CONTINUE 

CALC  *  CHA+CHO 
LL'.C  -  CL'J6+LLA 
vINC  s  iCHlC-ClLO/INC 
c 

C  BE5IN  ITERATIONS  OF  LIFE  CyCLE  COS7 


DO  40  M  =  1.  ITER 


“cost =a. 0 
L jSTU-8. D 
TlDD=0.0 
N'  18/  =  0.0 
DO  Cl  I  -  ,  Null 

place  total  nt  of  all  cost  elements  in  array  n«i8;  -  shuttle 
cost  ter  pound  element 

IF  1 1  ,E3.  18)  THEN 
DO  20  J  -  1,  SUM 

*T  =  NLIJi  +  BETA! IHBCJ) )*(MH< J>  -  NuJl) 

N «.  1 3 1  =  NT  *  NUB) 

C0  CONTINUE 

cNuiF 

C05T=  N !  I *  i  Cl  i  I +BETA  t  IB  ( I  /  >  ♦  (CH ( I )  -CL  1 1  i  /  >  i .  t 1 1 +DR !  *♦  1 1  yEAR+ 
1 L  ffi  '  1  i » 

C„bT  -  C Jb / +TCCST 
Cl  CONTINUE 
DO  05  J  -  I.l 
DO  C5  I  =  l.NUMO 


127 


1 


L-S I  J:  ■  Llu 1 1 '  +i?s !  n  ;St  . :  •  * !(.<•:  ■  i  .■  -  1  ** 

l  i  ?  E ARu  J  * ; 
rLU5:0  -  LllSTG+TCGSTG 
35  C'JNT  IVUE 

*LCt  =  TCQSUTCQSTG 
WRITE  U. *'  ITERATION  ♦:  ,N 

CuCuIaTE  "EAN.  5'AHE-ARD  DEVIATION,  AND  HI5T08RAN 

3D  4 1  i  =  1 ,  i  NC 
Cliii  =  Ci-LC+vINC+I 

41  CCNTiNUt 

DO  42  l  *  1,  INC 
»F  i T s_wC  ,iE.  Cl  ill)  ThEN 
’.COUNTS!)  =  ICQUNT 1 1)  +  1 
GO  TO  45 
ENDIF 

42  CONTINUE 

45  TC  -  TC  +  TLCC 
STLCC  =  TLCC*»2*5TLCC 
60  lGNT.NUE 
MEAN  =  TO/ [TER 

GO  =  1 1  ITER*STLCC-TC-**2)  / '  ITER* .  I  TER-i  J  > .  **8. 5 
WRITE  (1,18)  THE  HEAN  LIFE  CYCLE  COST  IS:  t  , MEAN 
WRITE  .1,18)  THE  5TANDARD  DEVIATION  IS:  S  .SO 
WRITE  11. ♦>  HISTObftAN: 

WRITE  ,1.10)  THE  LO  VALUE  ON  ThE  HISTOGRAM  IS:  i  ,CLLC 
WRITE  .1,18)  THE  HI  'VALUE  ]N  THE  HISTOGRAM  Is:  t  .OHIO 
bR I TE  1.1,18)  EACH  INCREMENT  ON  ThE  HISTOGRAM  IS  =  i  ,vlNC 
•RITE  ‘i,l0i  INCREMENT:  UPFER  LIMIT  vSlUE:  COUNT  WITHIN  INCREMENT: 
12  :GR"AT  .-12.0? 

WRI’E  it.'5)  1 1 , Cl 1 1 ; , IC ju’wT .  [ i ,  1  =  1.  INC) 

4  c,3rnat  'Jt.n.m.  i  .= i2.8.:ax. 

END 


'CjTInE  *3  CALCULATE  A  RANDOM  VAljE  F53P  BETA  RaNCT 3 LINS 
&“ siD  jN  sElEL'ICN  OF  SETA  CURVES  1  ""Ru  a 

FUNCT iCN  SETA  IE > 

X  =  RAND'01 

53  Tfl  ill, 12. 13, 14, 15. 16, 17, 18. 19). IB 

11  BETAM.45l9»Ml.  245*1**2*14. 898»l»*3-\3393*l«*4 
60  TQ  20 

12  BETAM. 5828*1-1. 7266*X-»*I+1 . 1 33'5*»**3 
GO  TO  20 

13  BETA=-0. 25015*1*2. 2a5*l**2-3. 2;la*i **3*2. 2205*i**4 

50  TO  28 

14  B£Ti=3.:7al*! i  +  .05;-9.a707*m2*12.954*l**3-a.2l56*i**4 
60  TO  20 

15  BET As2. 233*1- 3. fc09o*i** 2*2. 3339*1**3 
50  TO  20 

la  3ETA=-0. 37482*1*4. S703*»**:-9.  1390*i**3*5. 4641*1**4 
50  TO  28 

1T  BETA=2. a547# i «*.l i -5. 721!*l**2*a. 0833*1**3-2. 3001*1**1 


INPUT  FOR  ANAEROBIC  DIGESTOR  SYSTEM 


(Data  not  in  proper  field  locations) 

CARD  1  (PARAMETERS) 

NUM  ITER  INC  IYEAR  LASTYR  DR  NUMO  IYEARO  SEED 

18  10  100  18  50  0.10  5  20  0.327 


CARD  2-20  (ACQUISITION  DATA) 


ITEM 

CL(I) 

CH(I) 

IB(I) 

NCI) 

IYR(I) 

WL(I) 

WH(I) 

1WB(I) 

Tank 

13000 

20000 

1 

4 

1 

13000 

20000 

4 

Insulation 

6000 

7600 

2 

2 

1 

760 

800 

6 

Pump  1 

800 

1200 

8 

3 

1 

50 

75 

9 

Pump  2 

800 

1200 

8 

3 

1 

180 

200 

8 

Pump  3 

800 

1200 

8 

3 

1 

50 

75 

9 

Pump  4 

2200 

2600 

8 

3 

1 

50 

75 

9 

Compressor 

800 

1200 

8 

3 

1 

50 

75 

7 

Heat  Exchngr 

610 

1000 

1 

3 

1 

300 

500 

4 

Temp  Cntrl 

310 

350 

5 

3 

1 

10 

20 

9 

Pipe 

15587 

17100 

7 

1 

1 

20000 

23000 

6 

Valve 

270 

300 

9 

16 

1 

35 

45 

4 

Couplings 

13 

15 

9 

105 

1 

0 

0 

1 

Flame  Trap 

100 

200 

1 

5 

1 

75 

100 

7 

PRVAVB  w/FT  350 

500 

4 

5 

1 

150 

200 

7 

Sedmnt  Tank 

75 

150 

4 

3 

1 

40 

80 

8 

02  Tank 

6500 

10000 

1 

1 

1 

39357 

39357 

1 

Cnstruct’n 

30200 

50400 

4 

1 

2 

0 

0 

1 

Transport 

656 

2647 

6 

1 

2 

0 

0 

1 

OPERATIONS  DATA: 


CARD  1  -  5 
ITEM 


CLO(I)  CHI(I) 


IBE(I) 


Ops  Labor  1000000  1400000 


Maint. 

Energy 


10200 


02  Purchase  17200 


16600 

17400 


Transport  25818192  104177979  6 


INPUT  FOR  CONVENTIONAL  PROPELLANTS 
(Data  not  in  proper  field  locations) 

CARD  1  (PARAMETERS) 

NUM  ITER  INC  IYEAR  LASTYR  DR  NUMO  IYEARO  SEED 

1  500  100  18  50  0.10  3  20  0.327 

CARD  2  (AQUISITION  DATA) 

ITEM  CL(I)  CH(I)  IB(I)  N(I)  IYR(I)  WL(I)  WH(I)  IWB(I) 

Tank  3300  5100  1  1  2  0.0  0.0  1 

OPERATIONS  DATA: 

CARD  1  -  3 

ITEM  CLO(I)  CHI(I)  IBE(I) 

N204  67963  74758  4 

MMH  131808  164760  4 

Transport 


31232816 


126026317  6 


OUTPUT  FOR  ANEROBIG  DIGESTOR 
(All  data  within  first  49  intervals) 

THE  MEAN  LIFE  CYCLE  COST  13:  X  54771722. 

THE  STANDARD  DEVIATION  IS:  X  5309317. 

HISTOGRAM: 

THE  LO  VALUE  ON  THE  HISTOGRAM  IS:  X  31359310. 

THE  HI  VALUE  ON  THE  HISTOGRAM  IS:  $  125958890. 

EACH  INCREMENT  ON  THE  HISTOGRAM  IS  =  X  945996. 

INCREMENT:  UPPER  LIMIT  VALUE:  COUNT  WITHIN  INCREMENT: 

1  X  32305306. 

2  *•  33251302. 

3  X  34197297. 

4  *  35143293. 

5  X  36089289. 

6  X  37035285. 

7  *  37981281. 

8  X  38927276. 

9  X  39873272. 

10  X  40819268. 

11  X  41765264. 

12  X  42711260. 

13  X  43657255. 

14  X  44603251. 

15  X  45549247. 

16  X  46495243. 

17  X  47441239. 

18  X  48387234. 

19  X  49333230. 

20  X  50279226. 

22  X  52171218. 

23  X  53117213. 

24  X  54063209. 

25  X  55009205. 

26  X  55955201. 

27  X  56901197. 

28  X  57847192. 

29  X  58793188. 

30  X  59739184. 

31  X  60685180. 

32  *  61631176. 

33  X  62577171. 

34  X  63523167. 

35  X  64469163. 

36  X  65415159. 

37  X  66361155. 

38  X  67307150. 

39  X  68253146. 

40  X  69199142. 

41  X  70145138. 

42  X  71091134. 

43  X  72037129. 

44  X  72983125. 

45  *  73929121. 

46  X  74875117. 

47  X  75821113. 

48  X  76767108. 

49  X  77713104. 


►*  ^  i-*  w  is)  r-j  w  -b  C--I  ot  w  M  t  )  *-‘►^1-* 

SSQS^SQ^hHI.  Vj  CD  <l  CD  vl  10  OP  l-.J  O  W  A  LO  f-J  W  (J)  si  nJ  h  |j  (0  U  S  h  S  S  S  S  S  S  S  S  S  S 


OUTPUT  FOR  CONVENTIONAL  PROPELLANTS 
(Al!  data  within  first  49  intervals) 

THE  MEAN  LIFE  CYCLE  COST  IS:  X  7<bSt>17Q0. 

THE  STANDARD  DEVIATION  IS:  X  6025964. 

HISTOGRAM: 


THE  LO  VALUE  ON  THE  HISTOGRAM  IS: 
THE  HI  VALUE  ON  THE  HISTOGRAM  IS: 
EACH  INCREMENT  ON  THE  HISTOGRAM  IS 


INCREMENT: 

UPPER  LIMIT  VALUE: 

1 

$ 

45337116 

$ 

46666344 

X 

47995573 

4 

X 

49324801 

5 

X 

50654030 

6 

X 

51983258 

7 

X 

53312487 

a 

X 

54641715 

9 

X 

55970944 

10 

X 

57300172 

11 

X 

58629401 

12 

X 

59958629 

13 

X 

61287858 

14 

X 

62617086 

15 

X 

63946315 

16 

X 

65275543 

17 

X 

66604772 

IB 

X 

67934000 

19 

X 

69263229 

20 

X 

70592457 

21 

X 

71921686 

MO 

X 

73250914 

X 

74580143 

24 

X 

75909371 

25 

X 

77238600 

26 

X 

78567828 

27 

X 

79897057 

28 

X 

81226235 

29 

X 

82555514 

30 

X 

83884742 

31 

X 

85213971 

—  *-> 

X 

86543199 

*T  ~T 

X 

87872428 

34 

X 

89201656 

35 

X 

90530885 

36 

X 

918601 13 

37 

X 

93189342 

38 

X 

94518570 

39 

X 

95847799 

40 

X 

97177027 

41 

X 

98506256 

42 

X 

99835484 

43 

X 

101164710 

44 

X 

102493940 

45 

X 

103823170 

46 

X 

105152400 

47 

X 

106481630 

48 

X 

107810860 

49 

X 

109140080 

X  44007S87. 

X  176930740. 

=  X  1 329229 . 

COUNT  WITHIN  INCREMENT 


►*  >-*  M  r-J  M  LI  U1  -L  C-J  LI  ^  W  ^  h 

SS88S88SShSH*WU'J  W  vl  vj  M  CD  Lfl  t-  CD  Cl  <1  LI  W  Q  IJ  O'  C-4  -0  Ul  H  B  8  SI  S  S  S  Gl  S  8  Q  S  S  S 


Captain  Coral  C.  Fallstead  was  born  on  24  September  1949  in  Palo 
Alto,  California.  He  graduated  from  high  school  in  Redwood  City,  California,  in 
1967.  He  enlisted  in  the  Air  Force  in  April  1971,  and  was  stationed  as  a  Weather 
Equipment  Repairman  in  the  6th  Weather  Squadron  at  Tinker  AFB,  Oklahoma. 
He  was  accepted  to  the  Airman’s  Education  and  Commissioning  Program  in  1972, 
and  graduated  from  Oklahoma  State  University  in  May  1974,  receiving  the  degree 
of  Bachelor  of  Science  in  General  Engineering.  Upon  graduation  he  attended 
Officers  Training  School  at  Randolf  AFB,  Texas,  and  received  a  the  commission 
of  2nd  Lieutenant  in  the  United  States  Air  Force,  16  August  1974.  He  served  as 
Deputy  Commander,  and  Commander,  Missile  Combat  Crew,  571  Strategic  Missile 
Squadron,  390  Missile  Wing,  Davis  Monthan  AFB,  Arizona,  from  August  1974  till 
Feberuary  1979.  In  Feberuary  1979  he  was  stationed  at  NORAD/Cheyenne 
Mountain  Complex,  Colorado,  as  Space  Systems  Senior  Director  and  in  August 
1980  became  Chief,  Deep  Space  Analyst.  From  June  1981  through  June  1982  he 
was  Space  Systems  Director  at  the  13th  Missile  Warning  Squadron,  Clear,  Alaska. 
In  June  1982  he  became  the  Space  Surveillance  Program  Director  for  Ground  Based 
Radars  at  Space  Command,  Peterson  AFB,  Colorado.  He  entered  the  Space 
Operations  Management  Course  in  the  School  of  Engineering,  Air  Force  Institute 
of  Technology,  in  June  1984. 

Permanent  address:  5545  Galena  Dr. 


Colorado  Springs,  Colorado  80918 
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